Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,252,630 +G 98.7% coding (223/1050 nt) ltaE → threonine aldolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,252,6301.→G98.7% 3.2 ‑1.2 ‑0.4313coding (223/1050 nt)ltaEthreonine aldolase
Reads supporting (aligned to +/- strand):  new base (100/209):  ref base (0/4):  total (100/213)
Fisher's exact test strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.23e-01

Only 100 of 360 total aligned reads displayed.

GCATCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.G.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAGCTGGTGCTGGT  >  NC_014034/2252530‑2252724
                                                                                                     |                                                                                                
gCATCTGCGCGAGGTTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCt                                                                                                   <  2:5594988/101‑1 (MQ=255)
 cATCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g                                                                                                 <  2:8689718/101‑1 (MQ=255)
 cATCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g                                                                                                 >  1:6700076/1‑101 (MQ=255)
 cATCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g                                                                                                 <  1:5330101/101‑1 (MQ=255)
  aTCTGCGCGAGATTTTCGAGGCTCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.g                                                                                               <  2:3253691/101‑1 (MQ=255)
   tCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.gg                                                                                              <  1:158513/101‑1 (MQ=255)
   tCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.gg                                                                                              >  1:7737827/1‑101 (MQ=255)
    cTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.ggg                                                                                             <  1:8361106/101‑1 (MQ=255)
    cTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGAGGCGAATGCGCTCTCCTCTGCGGTGCTGACCCAGCCCT.g.ggg                                                                                             <  2:4369962/101‑1 (MQ=255)
    cTGCGCGAGATTTTCGAGGCGCCACAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGTTCTCGGGTGCTGTGCTGAACCAGCCCT.g.ggg                                                                                             >  2:1245629/1‑101 (MQ=255)
     tGCGTGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAGTGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.gggg                                                                                            <  1:1348247/101‑1 (MQ=255)
       cgcgAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.gggggg                                                                                          <  2:6324817/101‑2 (MQ=255)
       cgcgAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.gggggg                                                                                          <  1:438707/101‑2 (MQ=255)
        gcgAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.ggggggc                                                                                         <  2:709747/101‑3 (MQ=255)
         cgAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.g.ggggggcg                                                                                        <  1:6476429/101‑4 (MQ=255)
            gaTTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGAtt                                                                                     <  1:577323/101‑1 (MQ=255)
             aTTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGAttt                                                                                    <  2:7532617/101‑1 (MQ=255)
              gtttCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.G.GGGGGTCGAtttt                                                                                   <  2:2112961/100‑1 (MQ=255)
                ttCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCt                                                                                 <  2:4146865/101‑1 (MQ=255)
                   gAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTc                                                                              <  1:3680163/101‑1 (MQ=255)
                   gAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTc                                                                              <  1:7692743/101‑1 (MQ=255)
                   gAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTc                                                                              <  2:6880537/101‑1 (MQ=255)
                    aGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCa                                                                             <  1:7087324/101‑1 (MQ=255)
                    aGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCa                                                                             <  1:4058320/101‑1 (MQ=255)
                      gcgcCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCAcc                                                                           <  2:4828580/101‑1 (MQ=255)
                      gcgcCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCAcc                                                                           <  2:5836253/101‑1 (MQ=255)
                        gccgcAGTCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.GCGGGGG.CGATTTTCTGTCACCGc                                                                         <  2:2978923/101‑1 (MQ=255)
                        gccgcAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGc                                                                         >  1:7406021/1‑101 (MQ=255)
                        gccgcAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGc                                                                         >  2:7832762/1‑101 (MQ=255)
                        gccgcAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGc                                                                         >  2:9026343/1‑101 (MQ=255)
                         ccgcAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCa                                                                        <  2:3442771/101‑1 (MQ=255)
                         ccgcAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCa                                                                        <  2:6024152/101‑1 (MQ=255)
                            cAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACg                                                                     >  2:1877721/1‑101 (MQ=255)
                               gCCGTGGTCTATCTGGTGGCTACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGccc                                                                  <  1:7241141/101‑1 (MQ=255)
                               gCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGccc                                                                  <  1:2894979/101‑1 (MQ=255)
                                 cGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCat                                                                <  1:8038254/101‑1 (MQ=255)
                                 cGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCat                                                                <  1:510241/101‑1 (MQ=255)
                                 cGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCat                                                                <  1:1245661/101‑1 (MQ=255)
                                 cGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCat                                                                <  1:102517/101‑1 (MQ=255)
                                  gTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCata                                                               >  1:6284695/1‑101 (MQ=255)
                                  gTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCata                                                               <  1:7841784/101‑1 (MQ=255)
                                  gTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCata                                                               <  2:3299766/101‑1 (MQ=255)
                                  gTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCata                                                               >  2:7659710/1‑101 (MQ=255)
                                   tGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCatat                                                              >  2:1762858/1‑101 (MQ=255)
                                    ggTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATc                                                             >  2:8805368/1‑101 (MQ=255)
                                    ggTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATc                                                             <  1:6125408/101‑1 (MQ=255)
                                    ggTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATc                                                             >  1:1809536/1‑101 (MQ=255)
                                     gTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCg                                                            <  1:1439741/101‑1 (MQ=255)
                                         aTCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGa                                                        <  2:6736807/101‑1 (MQ=255)
                                          tCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGaa                                                       <  2:8496283/101‑1 (MQ=255)
                                           cTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAg                                                      <  1:1469200/101‑1 (MQ=255)
                                            tggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGa                                                     <  1:3321862/101‑1 (MQ=255)
                                            tggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGa                                                     <  2:6129608/101‑1 (MQ=255)
                                            tggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGa                                                     <  1:974544/101‑1 (MQ=255)
                                             ggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGAc                                                    >  1:5909268/1‑101 (MQ=255)
                                             ggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGAc                                                    <  1:5646869/101‑1 (MQ=255)
                                             ggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGAc                                                    <  2:5549569/101‑1 (MQ=255)
                                             ggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGAc                                                    <  1:4584524/101‑1 (MQ=255)
                                             ggtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGAc                                                    >  1:397782/1‑101 (MQ=255)
                                              gtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACg                                                   <  2:7464334/101‑1 (MQ=255)
                                              gtggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACg                                                   >  2:2539980/1‑101 (MQ=255)
                                               tggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGa                                                  >  1:2998712/1‑101 (MQ=255)
                                               tggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGa                                                  >  1:9126111/1‑101 (MQ=255)
                                               tggCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGa                                                  <  1:257038/101‑1 (MQ=255)
                                                 gCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAAt                                                <  2:7957373/101‑1 (MQ=255)
                                                  cGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATg                                               <  1:884616/101‑1 (MQ=255)
                                                   gACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATgc                                              >  2:1655000/1‑101 (MQ=255)
                                                    aCGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATgcg                                             <  2:3012502/101‑1 (MQ=255)
                                                     cGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATgcgg                                            <  2:8153771/101‑1 (MQ=255)
                                                        gCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGgcgc                                         >  1:4823002/1‑101 (MQ=255)
                                                          aCGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGccc                                       >  2:4062970/1‑101 (MQ=255)
                                                           cggcggCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGcccc                                      <  1:4734156/101‑1 (MQ=255)
                                                            ggcggcGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGccccc                                     <  2:7712799/101‑1 (MQ=255)
                                                             gcggcgAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCg                                    <  2:5230233/101‑1 (MQ=255)
                                                             gcggcgAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCg                                    <  2:1307197/101‑1 (MQ=255)
                                                                   aaTGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCt                              <  1:5407847/101‑1 (MQ=255)
                                                                   aaTGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCt                              >  2:362251/1‑101 (MQ=255)
                                                                   aaTGCGCTCTCCTGTGCGGTACTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAACTCt                              <  1:7200506/101‑1 (MQ=255)
                                                                    aTGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTa                             <  1:6414050/101‑1 (MQ=255)
                                                                     tGCGCTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTac                            <  1:3169319/101‑1 (MQ=255)
                                                                        gcTCTGCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACAcc                         <  1:106807/101‑1 (MQ=255)
                                                                        gcTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACAcc                         <  1:6758317/101‑1 (MQ=255)
                                                                        gcTCTCCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACAcc                         >  1:7563812/1‑101 (MQ=255)
                                                                         cgcgccTGTGCGGTGCTGACCCAGCTCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACcg                        <  2:2717571/97‑1 (MQ=255)
                                                                           ctcCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAGGACGAATGCGGGGCCCCCGAATTCTACACtggc                      >  1:2498659/1‑101 (MQ=255)
                                                                           ctcCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACcggc                      <  2:397774/101‑1 (MQ=255)
                                                                           ctcCTGTGCGGTGCTGACCCAGCCCTGG.GGGGG.CGAGTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGATTTCTACACcggc                      <  1:8168815/101‑1 (MQ=255)
                                                                               tgtgCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCt                  <  1:276216/101‑1 (MQ=255)
                                                                                 tgCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTcc                <  2:445793/101‑1 (MQ=255)
                                                                                 tgCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTcc                <  2:8317075/101‑1 (MQ=255)
                                                                                 tgCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTcc                <  1:214129/101‑1 (MQ=255)
                                                                                  gCGGTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCa               <  1:8658340/101‑1 (MQ=255)
                                                                                   cGGTGCTGGCCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGTAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCaa              <  2:2264910/101‑1 (MQ=255)
                                                                                     gTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgc            <  2:4700731/101‑1 (MQ=255)
                                                                                     gTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgc            <  2:7226651/101‑1 (MQ=255)
                                                                                     gTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgc            <  2:7405697/101‑1 (MQ=255)
                                                                                     gTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgc            <  2:3897984/101‑1 (MQ=255)
                                                                                     gTGCTGACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgc            >  1:820494/1‑101 (MQ=255)
                                                                                          gACCCAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgctggtg       >  2:1492088/1‑101 (MQ=255)
                                                                                             ccAGCCCTGG.GGGGG.CGATTTTCTGTCACCGCCACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAgctggtgctg    <  1:4562757/101‑1 (MQ=255)
                                                                                                     |                                                                                                
GCATCTGCGCGAGATTTTCGAGGCGCCGCAGGCCGTGGTCTATCTGGTGGCGACGGGCACGGCGGCGAATGCGCTCTCCTGTGCGGTGCTGACCCAGCCCT.G.GGGGG.CGATTTTCTGTCACCGCAACGCCCATATCGAGGAAGACGAATGCGGCGCCCCCGAATTCTACACCGGCGGCTCCAAGCTGGTGCTGGT  >  NC_014034/2252530‑2252724

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG