Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,884 T→C 5.2% T146A (ACG→GCG)  glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,8840TC 5.2% 337.6 / 6.9 115T146A (ACG→GCG) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base T (54/54);  new base C (4/2);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 117 total aligned reads displayed.

CGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGC  >  NC_000913/3570737‑3571033
                                                                                                                                                   |                                                                                                                                                     
cgGTAATGCCCACTTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTgc                                                                                                                                                     <  2:178610/150‑1 (MQ=255)
  gTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTgcgc                                                                                                                                                   <  1:403537/150‑1 (MQ=255)
  gTAATGCCCAAATTTTTCAAATAGTTGATCATCACCGGATTCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTgcgc                                                                                                                                                   <  2:442245/150‑1 (MQ=255)
   tAATGCCCAATTGTTTCAAATAGTTGATCATCACCGCATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCAGGGTGCAAGTACGGTAATCCTTAGACATGTGCTTCAAAAATAATCGTGCTGCCCCACGGCGTGcgcg                                                                                                                                                  >  1:366670/1‑150 (MQ=255)
       gCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGcggcgg                                                                                                                                              >  1:1452422/1‑150 (MQ=255)
        cccAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCggg                                                                                                                                             >  2:263424/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)
              tGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatca                                                                                                                                       <  1:922901/150‑1 (MQ=255)
              tGTTTCAAATAGTGGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCGCCGCGCACGCCGTGGGgcagca                                                                                                                                       <  2:1534151/150‑1 (MQ=255)
                 ttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCt                                                                                                                                    >  1:1779117/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                                                                                                                                 <  2:1393435/150‑1 (MQ=255)
                    aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                                 >  1:462431/1‑150 (MQ=255)
                    aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                                 >  1:1359445/1‑150 (MQ=255)
                    aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTAATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                                 <  1:492375/150‑1 (MQ=255)
                     aaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTccc                                                                                                                                >  2:840064/1‑150 (MQ=255)
                       tAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAg                                                                                                                              >  1:118769/1‑150 (MQ=255)
                       tAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAg                                                                                                                              <  2:1309002/150‑1 (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)
                                           cccGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACcaccac                                                                                                          >  2:156244/1‑150 (MQ=255)
                                                   cTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGGAACTACGTTAATCCTTTGACATGGTATTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCTTCTTGCGAGTCATAGTGGTCGGCCCCCAGTACGCAtt                                                                                                  >  1:298507/1‑150 (MQ=255)
                                                      tataAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCg                                                                                               >  1:1350377/1‑150 (MQ=255)
                                                       ataAGTGCCACGGATCTCGACCGGGATCTCCGGGGGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCgg                                                                                              <  2:1393266/150‑1 (MQ=255)
                                                          aGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGgcgc                                                                                           <  1:49167/150‑1 (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)
                                                                  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)
                                                                                  ctcCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTc                                                                   >  2:470838/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)
                                                                                             aGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTAt                                                        >  1:1269053/1‑150 (MQ=255)
                                                                                             aGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTAt                                                        >  1:1551348/1‑150 (MQ=255)
                                                                                              gTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATg                                                       <  2:1390007/150‑1 (MQ=255)
                                                                                                     aaTCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCg                                                <  1:460295/150‑1 (MQ=255)
                                                                                                         cTTTGACATGGGCTGCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTgca                                            <  2:868156/150‑1 (MQ=255)
                                                                                                               cATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCg                                      >  2:1707897/1‑150 (MQ=255)
                                                                                                                 tGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCggg                                    <  1:1664758/150‑1 (MQ=255)
                                                                                                                 tGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCggg                                    <  1:743665/150‑1 (MQ=255)
                                                                                                                  gggCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGt                                   <  2:1452422/150‑1 (MQ=255)
                                                                                                                  gggCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGt                                   <  1:1108035/150‑1 (MQ=255)
                                                                                                                    gCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCAGTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAGTATGACCGGCGTGCAGCAGCGGGttt                                 >  1:984638/1‑149 (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)
                                                                                                                          tAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATGAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATGATGACCGGCGTGCAGCAGCGGGTTATCTTTa                           >  1:277192/1‑150 (MQ=255)
                                                                                                                            aaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTATaaa                         >  2:170728/1‑150 (MQ=255)
                                                                                                                                   gTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCg                  >  1:288803/1‑150 (MQ=255)
                                                                                                                                   atgatgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCg                  >  1:214888/5‑150 (MQ=255)
                                                                                                                                    tgctgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGGCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTTTCGCGATAGGCAGGTGCATTATGCCCGGCGTGCAGCAGCGGGGTATCTTTATACTCcgcgg                 >  2:259712/1‑149 (MQ=255)
                                                                                                                                    tgatgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGt                 >  2:1255488/4‑150 (MQ=255)
                                                                                                                                     gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc                >  1:1163823/1‑150 (MQ=255)
                                                                                                                                     gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc                >  1:1236255/1‑150 (MQ=255)
                                                                                                                                     gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc                >  2:1217852/1‑150 (MQ=255)
                                                                                                                                     gctgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTc                <  2:544164/150‑1 (MQ=255)
                                                                                                                                       tgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCaa              <  2:1806730/150‑1 (MQ=255)
                                                                                                                                       tgcCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCaa              >  2:1578495/1‑150 (MQ=255)
                                                                                                                                          cccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATATTTAAACTCCCCGTAAAttt           >  2:1695090/1‑150 (MQ=255)
                                                                                                                                          cccACGGAGTGCGCGGCGCGGCATCATCTTCACATTCATATTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTAGCGATAGCCAGGTTCATTATGACCGGCGGGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAttt           <  2:298507/150‑1 (MQ=255)
                                                                                                                                           ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTg          >  1:1184666/1‑150 (MQ=255)
                                                                                                                                             aCGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGcc        <  2:1350377/150‑1 (MQ=255)
                                                                                                                                                   tGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCgcgcgc  >  1:1066245/1‑150 (MQ=255)
                                                                                                                                                   tGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAATATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCgcgcgc  >  1:581445/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGC  >  NC_000913/3570737‑3571033

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: