Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,827,198 G→C 5.4% A269A (GCG→GCC srlE → glucitol/sorbitol‑specific enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,827,1980GC 5.4% 286.0 / 6.3 92A269A (GCG→GCCsrlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (43/44);  new base C (2/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

ATCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATC  >  NC_000913/2827049‑2827346
                                                                                                                                                     |                                                                                                                                                    
aTCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAggcg                                                                                                                                                      >  1:1269991/1‑150 (MQ=255)
 tCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAggcgg                                                                                                                                                     <  2:1551782/150‑1 (MQ=255)
     ttccttccCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTg                                                                                                                                                 <  2:577397/150‑1 (MQ=255)
       ccttccCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCCGCCTGCg                                                                                                                                               >  2:203549/1‑150 (MQ=255)
         ttccCACAGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATAATCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGAc                                                                                                                                             <  2:1490658/150‑1 (MQ=255)
            ccACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGCCTTc                                                                                                                                          >  2:1613147/1‑150 (MQ=255)
            ccACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                          <  1:6043/150‑1 (MQ=255)
            ccACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                          >  1:1515458/1‑150 (MQ=255)
            ccACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCCGCCTGCgcgttg                                                                                                                                          <  1:1303133/150‑6 (MQ=255)
            ccACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCCGCCTGCgcgttg                                                                                                                                          <  2:1422729/150‑6 (MQ=255)
            ccACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATAATCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCACCGCGCAGGCGGCCTGCGCCTTc                                                                                                                                          >  1:1128439/1‑150 (MQ=255)
             cACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGCCTTCa                                                                                                                                         >  1:1456581/1‑150 (MQ=255)
             cACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCa                                                                                                                                         >  1:1232381/1‑150 (MQ=255)
             cACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCCGCCTGCgcgttga                                                                                                                                         <  1:203549/150‑7 (MQ=255)
                  cTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcg                                                                                                                                    >  1:56491/1‑150 (MQ=255)
                   attCACCTTACCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcgg                                                                                                                                   <  2:1045236/149‑1 (MQ=255)
                        ccTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                              >  2:585414/1‑150 (MQ=255)
                        ccTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcgcacggt                                                                                                                              >  2:782025/1‑150 (MQ=255)
                         cTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGtt                                                                                                                             >  2:662646/1‑150 (MQ=255)
                         cTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTAGCCATCAAAGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGtt                                                                                                                             >  1:1334019/1‑150 (MQ=255)
                             ccTCGGCCCAGGCGCAGTTATCGCACAGGTTATGGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTc                                                                                                                         >  1:1613743/1‑150 (MQ=255)
                                  gCCCAGGCGCAGTTATCGCACAGGTTATCGACGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGTCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCAGGTCGGTTTGTCGCTgg                                                                                                                    >  1:614/1‑150 (MQ=255)
                                       ggCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGaa                                                                                                               >  1:1481841/1‑150 (MQ=255)
                                        gcgcCGTTATCGCACAGGGTATCGGCGTATAGATTGGCGTGCCGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCCTCCCGGTCGGTTTGTCGCTGGCGGAAg                                                                                                              <  2:678311/150‑1 (MQ=255)
                                        gcgcAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAg                                                                                                              <  1:1616958/150‑1 (MQ=255)
                                        gcgcAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAg                                                                                                              <  1:941725/150‑1 (MQ=255)
                                          gcAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGcc                                                                                                            >  2:1674121/1‑150 (MQ=255)
                                           cAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGGCCTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCCTCCCGGGCGGTTTGTCGCTGGCGGACGccc                                                                                                           >  1:80924/1‑150 (MQ=255)
                                            aGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCg                                                                                                          <  2:1379966/150‑1 (MQ=255)
                                            aGTTATCGCACAGGTTATCGGCGTAATGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCg                                                                                                          >  2:664629/1‑150 (MQ=255)
                                             gTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGTCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGt                                                                                                         >  1:270320/1‑150 (MQ=255)
                                                aTCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAg                                                                                                      <  1:410488/150‑1 (MQ=255)
                                                aTCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAg                                                                                                      >  2:1668866/1‑150 (MQ=255)
                                                 tCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGTAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAgg                                                                                                     <  2:1249100/150‑1 (MQ=255)
                                                   gCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGac                                                                                                   >  1:930384/1‑150 (MQ=255)
                                                   gCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGac                                                                                                   >  2:991114/1‑150 (MQ=255)
                                                    cacaGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGACAGGaca                                                                                                  >  2:350710/1‑150 (MQ=255)
                                                     acaGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGacac                                                                                                 <  1:822999/150‑1 (MQ=255)
                                                     acaGGTTATCGGCGGATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGacac                                                                                                 <  1:534425/150‑1 (MQ=255)
                                                       aGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACgg                                                                                               <  1:46795/150‑1 (MQ=255)
                                                       aGGTTATCGGCGGATTGATTGGCGTGCAGATTGGTCTCGTCAATAATCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACgg                                                                                               <  1:167767/150‑1 (MQ=255)
                                                           tATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGATTTACCGGCACTGTTTGCCATCACCGCGCAGGCGGCCTGCGACTTCATACCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACAGGGTTcg                                                                                           >  1:1562214/1‑150 (MQ=255)
                                                               ggCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCAGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTc                                                                                       >  2:1207557/1‑150 (MQ=255)
                                                                 cGTATTGATTGGCGTGCAGATTGGGCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCgg                                                                                     <  1:1014231/150‑1 (MQ=255)
                                                                    tttgattgGCGTGCAGATTGTTCTCGGAAATATTCCGCCTCATCTGGATTTAACTTAACTGTTTGCCATCTACGCGCAGGCGGCCTGCGATTTCATCCCGGTCTGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGgtgt                                                                                  <  2:614/149‑1 (MQ=255)
                                                                        tttgGCGTGCTGATTTGTCTGGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCt                                                                              <  2:279877/149‑1 (MQ=255)
                                                                        attgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCt                                                                              <  2:56491/150‑1 (MQ=255)
                                                                         ttgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACTCGGTTCGCGTCGGTGTCCCtt                                                                             <  2:890635/150‑1 (MQ=255)
                                                                          tgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCCTCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCAGTCAGGACACGGTTCGCGTCGGTGTCCCTTc                                                                            <  1:713161/150‑1 (MQ=255)
                                                                          tgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTc                                                                            <  2:765643/150‑1 (MQ=255)
                                                                           ggCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTcg                                                                           >  1:401284/1‑149 (MQ=255)
                                                                           ggCGTGCAGATTGGGCTCGGCAATATTCCGCCGCATCTGGCTGTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCt                                                                           <  2:1480336/150‑1 (MQ=255)
                                                                            gcctGCAGATTGGGCGCGGCAATATGCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGGCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGGCGGTGTCCCTTCTg                                                                          <  2:915948/147‑1 (MQ=255)
                                                                            gCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTg                                                                          <  1:218194/150‑1 (MQ=255)
                                                                            gCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTg                                                                          >  1:1324655/1‑150 (MQ=255)
                                                                            gCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTGTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTg                                                                          <  2:915947/150‑1 (MQ=255)
                                                                            gCGTGCAGAATGGGCGCGTCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCCCGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGGCGGTGTCCCTTCTg                                                                          <  2:178618/150‑1 (MQ=255)
                                                                                gCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTTCt                                                                      >  1:597938/1‑150 (MQ=255)
                                                                                  aGATTGGTCTCGGCAATATTCCGCCGCAGCTGGCTTTACCGGCACTCTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGATCGCGTCGGTGGCACTTCAGGACTgt                                                                    >  1:258478/1‑149 (MQ=255)
                                                                                    aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGGAGGCGGCCTGCGACTTCTTCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTg                                                                  <  1:39529/150‑1 (MQ=255)
                                                                                        gTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGcc                                                              >  2:9904/1‑150 (MQ=255)
                                                                                         gctcGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCg                                                             <  2:512643/149‑1 (MQ=255)
                                                                                          ctcGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGc                                                            >  1:1650800/1‑150 (MQ=255)
                                                                                               cAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACCGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCCCttttt                                                       >  2:811471/1‑150 (MQ=255)
                                                                                                    ttCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCg                                                  >  2:640387/1‑150 (MQ=255)
                                                                                                     tCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCgg                                                 <  2:1232381/150‑1 (MQ=255)
                                                                                                           gcATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGcc                                           <  2:783034/150‑1 (MQ=255)
                                                                                                            cATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCg                                          >  1:901695/1‑150 (MQ=255)
                                                                                                                tGGCTTTACCGGCACTGTTGGCCATCATCTGGCTGGCGGCCTGGGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCGCTTCTGTACTGGTGAGCCGCTTTTTTACCGGCGCGCCGACTg                                      <  2:1334019/150‑1 (MQ=255)
                                                                                                                       aCCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTGTTTAACCGGCGCGCCGACTGTACTGAt                               >  1:20564/1‑150 (MQ=255)
                                                                                                                        gcGGAACTGTTTGCCATCAACGAGCAGGCGGCCGGCGACTTCATCCCGGTCGGTTTGGCGCTGGCGGAAGCCCGTCAGGACACGGGTCGCGGCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATc                              <  2:141869/149‑1 (MQ=255)
                                                                                                                        ccGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCACTTCTGTACTGGTGAGCCGCGTTTTAACCGGCGCGCCTACTGTAATGATc                              >  1:1267365/1‑150 (MQ=255)
                                                                                                                         cGGCACTGTTTGCCCTCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGGTCGCGTCGGTGTCCCTGCTGTACGGGTGCGCCGCTTTTTACCCGGCGCGCCGACTGTACTGCTCg                             >  1:1033142/1‑150 (MQ=255)
                                                                                                                          ggCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGc                            >  2:347560/1‑150 (MQ=255)
                                                                                                                          ggCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGc                            >  1:42653/1‑150 (MQ=255)
                                                                                                                             aCTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTg                         <  2:42653/150‑1 (MQ=255)
                                                                                                                              cTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTgg                        <  2:1481841/150‑1 (MQ=255)
                                                                                                                                 tttGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGttt                     <  2:1204956/150‑1 (MQ=255)
                                                                                                                                  tggCCATCAACGCGCAGGCGGCCTGCGACGTCATCCCGGTCGGTTTGTCGCGGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTg                    <  2:80924/148‑1 (MQ=255)
                                                                                                                                     ccATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGCCTGTACTGATCGCCTGGTTTGtct                 >  1:320726/1‑150 (MQ=255)
                                                                                                                                      cATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGtctc                <  1:1668866/150‑1 (MQ=255)
                                                                                                                                      cATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCGGGTTTGtctc                >  1:876125/1‑150 (MQ=255)
                                                                                                                                      cATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCACGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGGCTGGTTTGtctc                >  2:1651925/1‑150 (MQ=255)
                                                                                                                                       aTCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTAGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTcc               >  2:427560/1‑150 (MQ=255)
                                                                                                                                           agtcgcAGACCGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGGCCCTTCTGTACTGGTGAGCCGCTGATTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTACGGtt           >  1:595149/4‑150 (MQ=255)
                                                                                                                                            cgcgCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGttt          >  1:1664105/1‑150 (MQ=255)
                                                                                                                                            cgcgCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGttt          <  1:640387/150‑1 (MQ=255)
                                                                                                                                            cgcgCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGGACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGttt          >  2:207821/1‑150 (MQ=255)
                                                                                                                                              cgcCGGCGGCCTGCGACTTCATCCCGGTCGGGTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGGGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTta        <  2:164857/150‑1 (MQ=255)
                                                                                                                                               gcAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGACGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtat       <  1:1376915/150‑1 (MQ=255)
                                                                                                                                                cAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatc      <  1:664629/150‑1 (MQ=255)
                                                                                                                                                cAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGGCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatc      <  2:1684467/150‑1 (MQ=255)
                                                                                                                                                cAGGCGGCCTGCGACTTCATCCCGGTCGGGTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGGCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatc      <  2:491738/150‑1 (MQ=255)
                                                                                                                                                 aGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatct     <  2:1324655/150‑1 (MQ=255)
                                                                                                                                                  ggcggcCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatc      <  1:584421/148‑1 (MQ=255)
                                                                                                                                                  ggcggcCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatc      >  2:584421/1‑148 (MQ=255)
                                                                                                                                                  ggcggcCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGGCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatcta    <  2:1688336/150‑1 (MQ=255)
                                                                                                                                                    cggcCTGCGACTTCATCCCGGGCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatctatc  <  1:1396891/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATC  >  NC_000913/2827049‑2827346

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: