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.2% 1566.9 / 2.6 387pseudogene (1173/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (176/208);  total (179/208)
Fisher's exact test for biased strand distribution p-value = 9.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 409 total aligned reads displayed.

GGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑ATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAC  >  NC_008463/3566673‑3566966
                                                                                                                                                      |                                                                                                                                                
gggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCgg                                                                                                                                                   >  2:6718614/1‑150 (MQ=255)
 ggggTGACCTTCGGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCACTCGGACCCGTTCGG‑a                                                                                                                                                 <  2:5672721/150‑1 (MQ=255)
 ggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑a                                                                                                                                                 >  1:10306892/1‑150 (MQ=255)
 ggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑a                                                                                                                                                 >  1:11090399/1‑150 (MQ=255)
  gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGAGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑aa                                                                                                                                                >  1:3749931/1‑149 (MQ=255)
    gTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGTTGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑aatc                                                                                                                                              >  2:5201931/1‑147 (MQ=255)
     tGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑aatcg                                                                                                                                             <  1:1441889/150‑5 (MQ=255)
      gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑aatcgt                                                                                                                                            <  1:10053705/150‑6 (MQ=255)
       aaCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                                                            >  1:10565078/1‑150 (MQ=255)
              ggCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGt                                                                                                                                     <  1:10410949/150‑1 (MQ=255)
                   aTCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgc                                                                                                                                <  2:3130322/150‑1 (MQ=255)
                      gcgcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgag                                                                                                                             >  2:1711231/1‑150 (MQ=255)
                       cgcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgaga                                                                                                                            <  1:4973169/150‑1 (MQ=255)
                       cgcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAGGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgaga                                                                                                                            <  1:2171496/150‑1 (MQ=255)
                         ccGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                                                          >  2:1211505/1‑150 (MQ=255)
                         ccGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                                                          >  1:8680668/1‑150 (MQ=255)
                         ccGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                                                          <  2:7067932/150‑1 (MQ=255)
                         ccGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                                                          >  2:2426630/1‑150 (MQ=255)
                          cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAcc                                                                                                                         <  1:9401360/150‑1 (MQ=255)
                            gAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAt                                                                                                                       <  2:1379426/150‑1 (MQ=255)
                                gcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGccc                                                                                                                   <  2:10974562/150‑1 (MQ=255)
                                    ggCCAACGAGGCGCGCCCCTCGAAGTTCGCCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgct                                                                                                               >  2:9622059/1‑150 (MQ=255)
                                        aacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                                                           >  1:9089106/1‑150 (MQ=255)
                                         acgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgctc                                                                                                          >  2:804873/4‑150 (MQ=255)
                                          cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGCGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGGTCGGAATCGTCGACTGGTACCGCGAGCACCATGCCCGGCTCGCTcc                                                                                                         >  2:10366506/3‑150 (MQ=255)
                                           gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCg                                                                                                        >  2:5409387/2‑150 (MQ=255)
                                             ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGga                                                                                                      <  2:5091232/150‑1 (MQ=255)
                                             ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGga                                                                                                      <  2:680338/150‑1 (MQ=255)
                                                    cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCTACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt                                                                                               >  1:7484065/1‑150 (MQ=255)
                                                          aaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                                                         <  1:10845150/150‑1 (MQ=255)
                                                          aaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                                                         >  2:2994345/1‑150 (MQ=255)
                                                            gTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCaa                                                                                       <  1:5652338/150‑1 (MQ=255)
                                                            gTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCaa                                                                                       <  1:7566437/150‑1 (MQ=255)
                                                             ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAc                                                                                      <  1:10736959/150‑1 (MQ=255)
                                                             ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCAGCCGCGGCCAGTCGGACCCGTTCGGAATCGGCGACTGGTACCGCGAGAACCATGCCCGGCTCGCACCGGATACCGATGCCGCCAAc                                                                                      <  2:604650/150‑1 (MQ=255)
                                                              tCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                     <  1:4699151/150‑1 (MQ=255)
                                                              tCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                     <  2:944134/150‑1 (MQ=255)
                                                              tCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACGTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                                                     <  2:559394/150‑1 (MQ=255)
                                                                 gCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGc                                                                                  >  2:4992278/1‑150 (MQ=255)
                                                                 gCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGc                                                                                  <  2:918451/150‑1 (MQ=255)
                                                                   ccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGccc                                                                                <  1:1264184/150‑1 (MQ=255)
                                                                     gATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCca                                                                              <  2:3598501/150‑1 (MQ=255)
                                                                       tGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCcaca                                                                            <  1:3519751/150‑1 (MQ=255)
                                                                       tGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCcaca                                                                            <  2:8739574/150‑1 (MQ=255)
                                                                       tGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCcaca                                                                            <  1:5211395/150‑1 (MQ=255)
                                                                        gACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCacac                                                                           <  1:8875942/150‑1 (MQ=255)
                                                                           cgagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaa                                                                        >  2:58433/5‑150 (MQ=255)
                                                                            gagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaac                                                                       <  1:8385498/147‑1 (MQ=255)
                                                                             agaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacc                                                                      >  1:2355978/3‑150 (MQ=255)
                                                                              gaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacca                                                                     <  2:9821081/149‑1 (MQ=255)
                                                                                ccTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACcaaccaac                                                                   <  1:5986505/150‑1 (MQ=255)
                                                                                 cTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCaaccaacc                                                                  <  1:181478/150‑1 (MQ=255)
                                                                                  tGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAAccc                                                                 <  1:4271151/150‑1 (MQ=255)
                                                                                      gCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCAt                                                             <  1:2480126/150‑1 (MQ=255)
                                                                                        aTGTGGGAGAACGTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATcc                                                           <  2:1899766/150‑1 (MQ=255)
                                                                                         tgtgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATccc                                                          >  2:9496614/1‑150 (MQ=255)
                                                                                           tgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAg                                                        >  1:2647765/1‑150 (MQ=255)
                                                                                                     tCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCt                                              >  2:659851/1‑150 (MQ=255)
                                                                                                       cTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGa                                            <  2:7709705/150‑1 (MQ=255)
                                                                                                       cTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGa                                            <  1:281121/150‑1 (MQ=255)
                                                                                                        tGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGAt                                           <  1:8199501/150‑1 (MQ=255)
                                                                                                         gATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATc                                          <  1:9900304/150‑1 (MQ=255)
                                                                                                          ctCCCGCGCCACCCCCACAGCCGCGGCCAGGCGGCCCCGCCCGGAACCGTCGCCTGGTACCGCGAGAACCCCGCCCGGCTCGCTCCGGAGACCGAAGCCGCCAACCAGCCCCACACCACCCCACCCGCACCCCAGGCTACCGCCTGAAca                                         <  2:7916401/149‑1 (MQ=255)
                                                                                                           tCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcac                                        <  2:3525731/150‑1 (MQ=255)
                                                                                                           tCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcac                                        >  2:8286254/1‑150 (MQ=255)
                                                                                                            ccTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcaca                                       <  1:9352747/150‑1 (MQ=255)
                                                                                                             cTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcacac                                      <  1:6056106/150‑1 (MQ=255)
                                                                                                             cTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAGGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcacac                                      <  1:6598847/150‑1 (MQ=255)
                                                                                                               gcgcGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAg                                    >  2:8215221/1‑150 (MQ=255)
                                                                                                                    accaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCCACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                               <  1:6458718/150‑1 (MQ=255)
                                                                                                                    accaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                               >  2:4614718/1‑150 (MQ=255)
                                                                                                                    accaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                               >  1:451245/1‑150 (MQ=255)
                                                                                                                       accTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAcccc                            <  1:299118/150‑1 (MQ=255)
                                                                                                                        ccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACACGAAACCCCa                           <  2:7452359/150‑1 (MQ=255)
                                                                                                                        ccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAGGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCa                           <  2:606944/150‑1 (MQ=255)
                                                                                                                         ctctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCAt                          >  1:5755483/1‑150 (MQ=255)
                                                                                                                           ctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGa                        >  1:2544785/1‑150 (MQ=255)
                                                                                                                           ctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGa                        <  2:2796498/150‑1 (MQ=255)
                                                                                                                            agCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCAATGaa                       <  1:6761310/149‑1 (MQ=255)
                                                                                                                             gCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAAt                      <  1:4377643/150‑1 (MQ=255)
                                                                                                                               cgcgGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCCTGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCTACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATAc                    <  2:9299200/150‑1 (MQ=255)
                                                                                                                                  ggCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCa                 <  2:8656698/150‑1 (MQ=255)
                                                                                                                                   gCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAg                <  2:2109204/150‑1 (MQ=255)
                                                                                                                                   gCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAg                <  1:8431146/150‑1 (MQ=255)
                                                                                                                                   gCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACGCAg                >  2:6466295/1‑150 (MQ=255)
                                                                                                                                     cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAt              >  2:2570193/1‑150 (MQ=255)
                                                                                                                                       gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTg            <  1:9992204/150‑1 (MQ=255)
                                                                                                                                       gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTg            >  1:5959753/1‑150 (MQ=255)
                                                                                                                                       gTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTg            <  1:7227098/150‑1 (MQ=255)
                                                                                                                                        tCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGc           >  2:8951323/1‑150 (MQ=255)
                                                                                                                                         cGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGcc          >  2:4158930/1‑150 (MQ=255)
                                                                                                                                         cGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGcc          <  2:8438424/150‑1 (MQ=255)
                                                                                                                                         cGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAGGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCTGATTGcc          >  1:196210/1‑150 (MQ=255)
                                                                                                                                          ggACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGccc         <  1:9599050/150‑1 (MQ=255)
                                                                                                                                           gACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCa        >  1:9516961/1‑150 (MQ=255)
                                                                                                                                           gACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCa        <  1:207866/150‑1 (MQ=255)
                                                                                                                                             cccGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGa      <  1:412409/150‑1 (MQ=255)
                                                                                                                                              ccGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACATGAAACCCCATGAATACTCAGATTGCCCAGAt     <  2:7209471/150‑1 (MQ=255)
                                                                                                                                               cGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATc    >  2:8392416/1‑150 (MQ=255)
                                                                                                                                                 ttCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAc  <  1:1125148/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGG‑ATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAC  >  NC_008463/3566673‑3566966

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: ‑