Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,878 A→G 5.1% W148R (TGG→CGG)  glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,8780AG 5.1% 348.1 / 6.9 117W148R (TGG→CGG) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (57/54);  new base G (4/2);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-02

Only 100 of 118 total aligned reads displayed.

CCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCC  >  NC_000913/3570731‑3571027
                                                                                                                                                   |                                                                                                                                                     
ccAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGAGATGGGCTTCGTAAATGATGGTGCTGCCCCACg                                                                                                                                                     >  1:283813/1‑150 (MQ=255)
ccAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACg                                                                                                                                                     >  1:541921/1‑150 (MQ=255)
ccAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACg                                                                                                                                                     >  2:612286/1‑150 (MQ=255)
      cgGTAATGCCCACTTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTgc                                                                                                                                               <  2:178610/150‑1 (MQ=255)
        gTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTgcgc                                                                                                                                             <  1:403537/150‑1 (MQ=255)
        gTAATGCCCAAATTTTTCAAATAGTTGATCATCACCGGATTCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTgcgc                                                                                                                                             <  2:442245/150‑1 (MQ=255)
             gCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGcggcgg                                                                                                                                        >  1:1452422/1‑150 (MQ=255)
              cccAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCggg                                                                                                                                       <  2:144764/150‑1 (MQ=255)
               ccAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCgggg                                                                                                                                      >  2:797821/1‑150 (MQ=255)
               ccAATTGTTTCAAATAGTTGAGCATCACCGGATGCCCGACGGCTTTATAAGTGCCACGGAGCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCGTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCgggg                                                                                                                                      <  2:1581117/150‑1 (MQ=255)
                cAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGc                                                                                                                                     >  1:560893/1‑150 (MQ=255)
                  aTTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGATAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCGCCGCGCACGCCGTGGGGcag                                                                                                                                   <  1:1419101/150‑2 (MQ=255)
                  aTTGTTGCAAATACTTGAGCCTCCCCGGCTGCCCGAGGGCCTTCTAAGTGCCCCGGAGCTCGCCCGGGATCTCCGGGGGCACGTCCGGTCATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcat                                                                                                                                   <  1:1681917/150‑1 (MQ=255)
                   ttGTTTCAAATAGTTGATCCTCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAGCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                  <  1:1669647/150‑1 (MQ=255)
                   ttGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                  <  2:1807791/150‑1 (MQ=255)
                   ttGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                  <  1:840304/150‑1 (MQ=255)
                    tGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatca                                                                                                                                 <  1:922901/150‑1 (MQ=255)
                       ttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCt                                                                                                                              >  1:1779117/1‑150 (MQ=255)
                       ttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCt                                                                                                                              >  2:1634560/1‑150 (MQ=255)
                        tCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGAGTGCGAGGCGGGGCATCATCtt                                                                                                                             >  2:751086/1‑150 (MQ=255)
                          aaaTAGTTGATCATCTCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                           <  1:499246/150‑1 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCTCGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                           <  2:627373/150‑1 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCGTCATCTTcc                                                                                                                           >  1:1203428/1‑150 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                           >  1:1359445/1‑150 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                           <  2:1393435/150‑1 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                           >  1:462431/1‑150 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTAATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                           <  1:492375/150‑1 (MQ=255)
                           aaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTccc                                                                                                                          >  2:840064/1‑150 (MQ=255)
                             tAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAg                                                                                                                        <  2:1309002/150‑1 (MQ=255)
                             tAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAg                                                                                                                        >  1:118769/1‑150 (MQ=255)
                             tAGTTGATCATCACCGGATGCACGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAg                                                                                                                        <  2:895539/150‑1 (MQ=255)
                                  gATCATCCCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATa                                                                                                                   >  2:1252682/1‑150 (MQ=255)
                                  gATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATa                                                                                                                   >  2:45486/1‑150 (MQ=255)
                                   atcatcACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAg                                                                                                                  <  2:819033/150‑1 (MQ=255)
                                     catcaCCGGATGCACGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTg                                                                                                                <  1:1634386/150‑1 (MQ=255)
                                       tcaCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGAt                                                                                                              >  2:305073/1‑150 (MQ=255)
                                        caCCGGATGCCCGAGGGCTTTATAAGTGCCCCGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATc                                                                                                             <  2:645974/150‑1 (MQ=255)
                                         aCCGGATGCCCGAGGGCTTTATAAGTTCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCa                                                                                                            <  2:1305272/150‑1 (MQ=255)
                                            ggATGCCCGAGGGGTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAacc                                                                                                         <  2:331543/150‑1 (MQ=255)
                                               tGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAaccacc                                                                                                      <  1:115893/150‑1 (MQ=255)
                                               tGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAaccacc                                                                                                      <  1:1060700/150‑1 (MQ=255)
                                               tGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAaccacc                                                                                                      <  1:1323089/150‑1 (MQ=255)
                                                 cccGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACcaccac                                                                                                    >  2:156244/1‑150 (MQ=255)
                                                  ccGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCACCCACCCCt                                                                                                   >  2:1433889/1‑150 (MQ=255)
                                                         cTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGGAACTACGTTAATCCTTTGACATGGTATTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCTTCTTGCGAGTCATAGTGGTCGGCCCCCAGTACGCAtt                                                                                            >  1:298507/1‑150 (MQ=255)
                                                         cTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCAtt                                                                                            >  2:1318456/1‑150 (MQ=255)
                                                             ataAGTGCCACGGATCTCGACCGGGATCTCCGGGGGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCgg                                                                                        <  2:1393266/150‑1 (MQ=255)
                                                                 gTGCCACGGATCTCGCCCGGGATCTCCGGGTGAAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTTCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCa                                                                                    >  1:1172794/1‑150 (MQ=255)
                                                                 gTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCa                                                                                    >  1:1614146/1‑150 (MQ=255)
                                                                   gCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAAt                                                                                  >  1:818938/1‑150 (MQ=255)
                                                                    ccACGGATCTCGACCGGGATCTCCGGGTGCAAGTACATTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATg                                                                                 >  1:1461316/1‑150 (MQ=255)
                                                                       cGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATggcg                                                                              <  2:118769/150‑1 (MQ=255)
                                                                       cGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTGCCCAGTCATTGTGATCAACCACCACTACGCATTTCGGCGCAATggcg                                                                              >  2:19352/1‑150 (MQ=255)
                                                                        ggATCTCGACCGGGATCTCCGGGTGCAAGTAGGTTAATCCTTTGTCATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTACCAGTCATAGTGATCAACCACCACTAGGCATTTCGGCGCAATggcgg                                                                             <  2:1041473/150‑1 (MQ=255)
                                                                          aTCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGgcggcg                                                                           <  1:22455/150‑1 (MQ=255)
                                                                             tcGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTg                                                                        >  1:258301/1‑150 (MQ=255)
                                                                                aCCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGTCGTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTcg                                                                     <  1:78694/150‑1 (MQ=255)
                                                                                 ccGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTcgc                                                                    <  1:515799/150‑1 (MQ=255)
                                                                                   gggATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGa                                                                  >  1:88862/1‑150 (MQ=255)
                                                                                    ggATCTCCGGTTGCAAGTCAGTTAATCCATTGACAAGGGCTTCATAAATTATGATGCTCCCCCACGGCGTGCGCGGCGTGGCATCATCTTCACAGTCATAGTGATCAACCACCACAACGCATTTCGGCGCAATGGCGGCGTTGTCGCGAt                                                                 <  1:307164/150‑1 (MQ=255)
                                                                                      aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGTTCAACCACCACTACGAATTTCGGCGCATTGGCGGCGTTGTCGAGCTAg                                                               >  1:287465/1‑150 (MQ=255)
                                                                                       tctcCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGt                                                              >  2:1342936/1‑150 (MQ=255)
                                                                                         tcCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGCGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCa                                                            <  1:136287/150‑1 (MQ=255)
                                                                                          ccGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAg                                                           <  2:541921/150‑1 (MQ=255)
                                                                                              gTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTc                                                       >  1:1301398/1‑150 (MQ=255)
                                                                                                gCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAt                                                     <  2:48527/150‑1 (MQ=255)
                                                                                                 cAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTCCGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAtt                                                    >  1:1168670/1‑150 (MQ=255)
                                                                                                 cAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAtt                                                    <  2:649291/150‑1 (MQ=255)
                                                                                                 cAAGTACGTTAATCCTTTGACATGGGCATCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAtt                                                    <  2:556271/150‑1 (MQ=255)
                                                                                                  aaGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTa                                                   >  2:733458/1‑150 (MQ=255)
                                                                                                  aaGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGGGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTa                                                   <  1:515619/150‑1 (MQ=255)
                                                                                                   aGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTAt                                                  >  1:1269053/1‑150 (MQ=255)
                                                                                                   aGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTAt                                                  >  1:1551348/1‑150 (MQ=255)
                                                                                                      aCGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGAc                                               <  1:407036/150‑1 (MQ=255)
                                                                                                           aaTCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCg                                          <  1:460295/150‑1 (MQ=255)
                                                                                                               cTTTGACATGGGCTGCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTgca                                      <  2:868156/150‑1 (MQ=255)
                                                                                                                       tGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCggg                              <  1:1664758/150‑1 (MQ=255)
                                                                                                                        gggCTTCATAAATGTTGGTGCTGCCGCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGAACAACCACCACTACGCATTTCGGCGCCATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGt                             <  1:285105/150‑1 (MQ=255)
                                                                                                                        gggCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGt                             <  1:1108035/150‑1 (MQ=255)
                                                                                                                        gggCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGt                             <  2:1452422/150‑1 (MQ=255)
                                                                                                                             tCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCt                        <  1:446115/150‑1 (MQ=255)
                                                                                                                              cgtAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCtt                       >  2:1005089/3‑150 (MQ=255)
                                                                                                                              cATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATAATCTTACCAGTAATACTGATCAACCCCCACTACGCATTTCGGCGCAATGGCGGCGTTGTAGCGATAATCAGGTTCATTATGACCCGCGTGAAGCAGCGGGGTATAtt                       >  2:1318504/1‑150 (MQ=255)
                                                                                                                              aaTAAATGATGGTGCTGCCCCAAGGCGTGCGGGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCtt                       >  2:669670/2‑150 (MQ=255)
                                                                                                                                tctatgatgGAGCTGCCCCACGGCGTGCGCGTCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACGACTCCGCCTTTCGGCGCAAGGGCGGCGTTGTCGCGATAGTAAGGGTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTa                     <  1:263424/147‑1 (MQ=255)
                                                                                                                                  aaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTATaaa                   >  2:170728/1‑150 (MQ=255)
                                                                                                                                        ggTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGAAAGTCAGGTTCATTAAGACCGGCGTGCAGCAGCGGGTTAGCTTTAAACTccgc             >  2:1116271/1‑148 (MQ=255)
                                                                                                                                         gTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCg            >  1:288803/1‑150 (MQ=255)
                                                                                                                                         atgatgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCg            >  1:214888/5‑150 (MQ=255)
                                                                                                                                          tgctgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGGCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTTTCGCGATAGGCAGGTGCATTATGCCCGGCGTGCAGCAGCGGGGTATCTTTATACTCcgcgg           >  2:259712/1‑149 (MQ=255)
                                                                                                                                           gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc          >  2:1217852/1‑150 (MQ=255)
                                                                                                                                           gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc          <  2:544164/150‑1 (MQ=255)
                                                                                                                                           gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc          >  1:1236255/1‑150 (MQ=255)
                                                                                                                                           gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc          >  1:1163823/1‑150 (MQ=255)
                                                                                                                                             tgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCaa        >  2:1578495/1‑150 (MQ=255)
                                                                                                                                             tgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCaa        <  2:1806730/150‑1 (MQ=255)
                                                                                                                                                cccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATATTTAAACTCCCCGTAAAttt     >  2:1695090/1‑150 (MQ=255)
                                                                                                                                                cccACGGAGTGCGCGGCGCGGCATCATCTTCACATTCATATTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTAGCGATAGCCAGGTTCATTATGACCGGCGGGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAttt     <  2:298507/150‑1 (MQ=255)
                                                                                                                                                 ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTg    >  2:447862/1‑150 (MQ=255)
                                                                                                                                                   aCGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGcc  <  2:1350377/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCC  >  NC_000913/3570731‑3571027

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