Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,823 +A pseudogene (1174/1275 nt) PA14_RS16215 → fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,566,8221.A99.6% 1859.4 / ‑0.8 453pseudogene (1173/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (210/240);  total (211/242)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 478 total aligned reads displayed.

CGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCACCCA  >  NC_008463/3566672‑3566969
                                                                                                                                                       |                                                                                                                                                   
cGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCgg                                                                                                                                                      <  1:5052430/151‑1 (MQ=255)
   gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaat                                                                                                                                                  >  2:869773/1‑149 (MQ=255)
   gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaat                                                                                                                                                  <  1:8104576/151‑3 (MQ=255)
     gTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatcg                                                                                                                                                <  2:11157934/151‑5 (MQ=255)
      tGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatcgt                                                                                                                                               <  1:10226235/151‑6 (MQ=255)
      tGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatcgt                                                                                                                                               >  1:9245873/1‑146 (MQ=255)
       gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                                                               >  2:8989434/1‑151 (MQ=255)
       gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                                                               <  1:3041591/151‑1 (MQ=255)
       gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                                                               <  1:2396227/151‑1 (MQ=255)
         aCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGa                                                                                                                                             >  1:4267636/1‑151 (MQ=255)
         aCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGa                                                                                                                                             >  2:616208/1‑151 (MQ=255)
           ttCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACt                                                                                                                                           >  1:10685152/1‑151 (MQ=255)
            tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                                                          <  1:6210316/151‑1 (MQ=255)
              tGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGt                                                                                                                                        <  2:6505230/151‑1 (MQ=255)
                   cATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgc                                                                                                                                   <  1:5205432/151‑1 (MQ=255)
                    aTCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgcg                                                                                                                                  >  2:3666950/1‑151 (MQ=255)
                       gcgcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgaga                                                                                                                               <  1:7094901/151‑1 (MQ=255)
                         gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                                                             >  1:4879744/1‑151 (MQ=255)
                           cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                                                                                           >  2:451157/1‑151 (MQ=255)
                            ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGCAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAt                                                                                                                          <  1:9482627/151‑1 (MQ=255)
                                cgcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGccc                                                                                                                      >  2:10642203/1‑151 (MQ=255)
                                 gcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCg                                                                                                                     <  2:2288342/151‑1 (MQ=255)
                                 gcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCg                                                                                                                     <  2:8005040/151‑1 (MQ=255)
                                 gcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCg                                                                                                                     <  1:2958991/151‑1 (MQ=255)
                                        caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                                                              <  1:7891421/151‑1 (MQ=255)
                                        caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                                                              <  1:5890796/151‑1 (MQ=255)
                                        caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                                                              <  2:3199485/151‑1 (MQ=255)
                                          acgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTcc                                                                                                            <  2:6240744/148‑1 (MQ=255)
                                           cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGGGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGGCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCg                                                                                                           >  1:7669642/3‑151 (MQ=255)
                                           cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCg                                                                                                           <  1:9571851/149‑1 (MQ=255)
                                            gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg                                                                                                          <  1:3623893/150‑1 (MQ=255)
                                             aGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGGTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCGCGCTCCGga                                                                                                         >  2:4888503/1‑151 (MQ=255)
                                              ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgag                                                                                                        <  2:1626377/151‑1 (MQ=255)
                                                cgcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAc                                                                                                      <  1:6581361/151‑1 (MQ=255)
                                                  cgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCg                                                                                                    <  2:513855/151‑1 (MQ=255)
                                                  cgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCg                                                                                                    >  2:10991323/1‑151 (MQ=255)
                                                  cgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCg                                                                                                    <  2:1870793/151‑1 (MQ=255)
                                                    ccccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt                                                                                                  <  2:5156837/151‑1 (MQ=255)
                                                     cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATg                                                                                                 >  2:5387251/1‑151 (MQ=255)
                                                     cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATg                                                                                                 <  2:1319377/151‑1 (MQ=255)
                                                     cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATg                                                                                                 <  1:6661816/151‑1 (MQ=255)
                                                     cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATg                                                                                                 >  2:3552084/1‑151 (MQ=255)
                                                          gaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                                                            >  1:1208697/2‑151 (MQ=255)
                                                          gaaGTTCGGCCCGATGACCGAGAACTGGGGCATGTGGGATAACTACCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGATCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                                                            >  2:890492/2‑151 (MQ=255)
                                                              ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                        <  1:1946192/151‑1 (MQ=255)
                                                              ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                        <  1:8722628/151‑1 (MQ=255)
                                                               acGGCCCGATGACCGAGACATGGGGCATGAGGGAGAACTACCTGATCCTGCGCGACCACCTCAGCCGCGGCCACACGGACCCGTACGGAAACGTCCACTGGTACCGAGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCa                                                                                       >  2:8325393/2‑151 (MQ=255)
                                                                cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                                                      <  2:4668053/151‑1 (MQ=255)
                                                                    ccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcccc                                                                                  >  2:3744889/1‑151 (MQ=255)
                                                                    ccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcccc                                                                                  >  2:9856026/1‑151 (MQ=255)
                                                                     cGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCca                                                                                 >  2:10764721/1‑151 (MQ=255)
                                                                         gACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACacc                                                                             <  1:10461886/151‑1 (MQ=255)
                                                                            cgagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaac                                                                          >  2:7348065/5‑151 (MQ=255)
                                                                             gagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacc                                                                         <  2:6943151/148‑1 (MQ=255)
                                                                             gagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacc                                                                         <  2:9195535/148‑1 (MQ=255)
                                                                              agaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacca                                                                        <  1:4671612/149‑1 (MQ=255)
                                                                                        cATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATcc                                                              <  1:1477055/151‑1 (MQ=255)
                                                                                             gggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGc                                                         >  1:10065369/1‑151 (MQ=255)
                                                                                             gggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGc                                                         >  1:350912/1‑151 (MQ=255)
                                                                                              ggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCt                                                        <  2:9740945/151‑1 (MQ=255)
                                                                                                  aaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCg                                                    <  2:8161153/151‑1 (MQ=255)
                                                                                                   aCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGc                                                   <  1:8310445/151‑1 (MQ=255)
                                                                                                        cTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGAt                                              <  2:6746066/151‑1 (MQ=255)
                                                                                                         tgctCCTGCGCGCCCCCCTCAGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCACCCCGCATCCCAGGCTACCGCCTGATc                                             <  2:4340263/148‑1 (MQ=255)
                                                                                                          gATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATca                                            <  2:7133163/151‑1 (MQ=255)
                                                                                                          gATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATca                                            >  1:3451142/1‑151 (MQ=255)
                                                                                                          gATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATca                                            >  1:4170706/1‑151 (MQ=255)
                                                                                                           aTCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcac                                           <  1:8963035/151‑1 (MQ=255)
                                                                                                              cTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcacaca                                        >  2:9357661/1‑151 (MQ=255)
                                                                                                               tGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAg                                       <  2:2601254/151‑1 (MQ=255)
                                                                                                                  gcgaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCTCACAGGaa                                    <  1:326626/148‑1 (MQ=255)
                                                                                                                    gaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                                  <  1:6414888/150‑1 (MQ=255)
                                                                                                                    gaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                                  <  2:2073444/150‑1 (MQ=255)
                                                                                                                     accaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAcc                                 >  1:11024937/1‑151 (MQ=255)
                                                                                                                       caccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAcccc                               <  2:5911641/151‑1 (MQ=255)
                                                                                                                         ccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCAt                             >  1:3379351/1‑151 (MQ=255)
                                                                                                                          ctctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATg                            <  1:4500024/151‑1 (MQ=255)
                                                                                                                             tGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAAt                         <  2:8492863/151‑1 (MQ=255)
                                                                                                                               ccGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATAc                       <  2:2972346/151‑1 (MQ=255)
                                                                                                                                cgcgGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACt                      <  1:2120822/151‑1 (MQ=255)
                                                                                                                                   ggCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAg                   <  1:377646/151‑1 (MQ=255)
                                                                                                                                   ggCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAg                   >  1:7567467/1‑151 (MQ=255)
                                                                                                                                   ggCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAg                   <  1:9827860/151‑1 (MQ=255)
                                                                                                                                      cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAtt                >  1:5011014/1‑151 (MQ=255)
                                                                                                                                      cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAtt                >  1:8615056/1‑151 (MQ=255)
                                                                                                                                      cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAtt                >  1:5604448/1‑151 (MQ=255)
                                                                                                                                       aGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACTCCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTg               >  1:2737424/1‑151 (MQ=255)
                                                                                                                                       aGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTg               >  1:5432324/1‑151 (MQ=255)
                                                                                                                                        gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGc              >  1:1685708/1‑151 (MQ=255)
                                                                                                                                        gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGc              <  2:5274305/151‑1 (MQ=255)
                                                                                                                                        gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGc              <  2:7592624/151‑1 (MQ=255)
                                                                                                                                         tCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGcc             >  2:4132860/1‑151 (MQ=255)
                                                                                                                                         tCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGcc             >  2:10538844/1‑151 (MQ=255)
                                                                                                                                          cGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCAGGAATACTCAGATAGccc            >  2:9544918/1‑151 (MQ=255)
                                                                                                                                           ggACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCa           >  2:1441934/1‑151 (MQ=255)
                                                                                                                                           ggACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCa           <  1:998368/151‑1 (MQ=255)
                                                                                                                                            gACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAg          <  1:2190584/151‑1 (MQ=255)
                                                                                                                                             aCCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGa         >  2:5639218/1‑151 (MQ=255)
                                                                                                                                             aCCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGa         >  2:10282676/1‑151 (MQ=255)
                                                                                                                                               ccGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATc       <  2:7806496/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCACCCA  >  NC_008463/3566672‑3566969

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: