Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,105,4990AC41.8% 151.6 / 12.3 92Y356S (TAC→TCC) mltCmembrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base A (37/16);  new base C (0/38);  total (37/55)
Fisher's exact test for biased strand distribution p-value = 5.39e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATC  >  NC_000913/3105353‑3105645
                                                                                                                                                  |                                                                                                                                                  
gcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccg                                                                                                                                                 <  1:964307/150‑1 (MQ=255)
gcGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAAGCCTAccg                                                                                                                                                 <  1:1215946/150‑1 (MQ=255)
    gTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgcc                                                                                                                                            <  2:575342/151‑1 (MQ=255)
     tCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTAccgccgccg                                                                                                                                           >  1:910952/1‑151 (MQ=255)
            gAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTc                                                                                                                                    <  1:440023/151‑1 (MQ=255)
             aaggATAAGATTCAGACTGCGAATATAATTAACCCCATGACGCCGGCCGATGTTTACCAAACGCTGACGACCCGCCATCCCCCTGCGGAATCTCGCCGTTATCTTTAAAAAGTTAATACCGCGCAAAAATCCTACCGCCGCCGATAATTcc                                                                                                                                   <  2:973124/148‑1 (MQ=255)
                  tAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTa                                                                                                                               >  1:503098/1‑150 (MQ=255)
                  tAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTa                                                                                                                               >  2:1369918/1‑150 (MQ=255)
                   aaGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAc                                                                                                                             >  2:1281740/1‑151 (MQ=255)
                    aGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTaa                                                                                                                              >  2:39021/1‑149 (MQ=255)
                        tCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGccc                                                                                                                        >  2:465555/1‑151 (MQ=255)
                           ggCTGCCAATATTATTAACACCATGCCGCCGGGCGTTGTTTATCAGCCGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCTTTATCTTTATAAAGTGAATCCCGCCCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                                     <  1:967266/151‑1 (MQ=255)
                              tGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACg                                                                                                                  >  2:990372/1‑151 (MQ=255)
                                caaaTATTATTAACCCCATGCCGCCGGGCGATGTTTATCAGCCGCTGACGACCCGCCATCCCTCGGCGGAATCTCGCCGTTATCTTTATAAAGGGAATCCCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGAc                                                                                                                   <  1:744230/146‑1 (MQ=255)
                                caaaTATTATTAACCCAATGAGCCCGCGCTTTTTTTATCAGACGCGAACGCCCCGCCCTCCCTCTGCGGAATCCCCCCTTTATATTTAAAAAGTAAACCCCGCGCAAAAACCTTCCCGCCCCCGATATTTCCATTAACCGCCCCTGACGAt                                                                                                                <  2:1350972/149‑1 (MQ=255)
                                   atatTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGc                                                                                                              >  2:1219112/1‑150 (MQ=255)
                                    tattatTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                            >  1:762012/1‑151 (MQ=255)
                                         ttaACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgggg                                                                                                       >  1:384149/1‑151 (MQ=255)
                                         taaccacCATGCCCCGGGGGAATTTTTATAAGCCCCTGACGCCCCCCCATCCCTCGGCGAAACCTCGCCTTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAggg                                                                                                        <  1:591898/146‑1 (MQ=255)
                                           aaCCCCATGACGCCGGGCGATTTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCTTTATCTTTATAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCa                                                                                                     <  1:844759/151‑1 (MQ=255)
                                           aaCACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGc                                                                                                      >  2:129963/1‑150 (MQ=255)
                                               ccATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGTCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGTTGCTCAGGGGCaaa                                                                                                   >  2:679434/1‑149 (MQ=255)
                                               caatgcccccGGGGGTTTTTTACAAGACGCTGACGACCCCCCACCCCTCTGCGAAATCCCCCCTTTATCTTTATAAATTGAATCCCGCGAAAAAATCCGCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaa                                                                                                  <  2:1089060/142‑1 (MQ=255)
                                               caatgcccccGGGCGTTGTTTATCAGACTCTGACGCCCCGCCACCCCTCTGCGAAATCCCGCCGTTATCTTTATAAATTGAATCCCGCGAAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaa                                                                                                    <  2:995753/140‑1 (MQ=255)
                                                    aCGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGt                                                                                              >  1:1356964/1‑149 (MQ=255)
                                                     cGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                                            <  1:317044/150‑1 (MQ=255)
                                                         gggCGATGTTTATCGGCCGCTGACGACCCGCCACCCCTCTGCGGAATCTCGCCTTTATTTTTATAAAGTGAATACCGCGCAAAAACCCTACCGCCGACGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCca                                                                                        <  2:281852/150‑1 (MQ=255)
                                                         gggCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                       >  1:253809/1‑151 (MQ=255)
                                                          ggCGATTTTTACCAGCCGCTGACGCCCCGCCATCCCTCGGCGGAATCTCGCCTTTATCTTTATAAAGTGAATACCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcaca                                                                                      <  1:386574/151‑1 (MQ=255)
                                                             gATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATc                                                                                    >  1:834082/1‑150 (MQ=255)
                                                              aTGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcac                                                                                  >  2:1140635/1‑151 (MQ=255)
                                                                 ttaaccggACCCTGACGCCCCGCCACCCCCCTGCAGATTCTAGCCGTTATCTTAATAAAGGAAATACCACACAAAAATCCACCCCCCGCTGATAATGCCATTACCCGACCCTGACGATGCTCAGGGGCAAAAATGTTATCCTCATCACaa                                                                                <  2:274741/143‑1 (MQ=255)
                                                                  ttATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                              >  1:322682/1‑151 (MQ=255)
                                                                   tATCAGCCGCTGCCGCCCCCCCACCCCTCTGGGAAATCTCGCCGTTATTTTTAAAAAGGAAAACCCGCGCAAAAATCCTCCCGCCGCCGTTAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                              <  2:836427/150‑1 (MQ=255)
                                                                       aGACGCTGACGCCCCGCCATCCCTATGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                          <  2:679570/150‑1 (MQ=255)
                                                                       aGACGCTGACGACCCGCCACCCCTCTGCGGAATCTCGCCTTTATCTTTATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                                           <  1:13179/149‑1 (MQ=255)
                                                                        gACGCTGCCGCCCCGCCACCCCCCTGCGGAATCCCGCGTTTATTTTTATAAAGTAAATCCCGCGAAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                                                                         <  2:1005154/150‑1 (MQ=255)
                                                                          cGCTGACGCCCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                                            <  1:114553/145‑1 (MQ=255)
                                                                          cGCTGACGCCCCGCCACCCCTCTGCGGAATCTCGCCGTTATTTTTATAAAGTGAATACCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                                                           <  2:904610/146‑1 (MQ=255)
                                                                            cTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                     >  2:1275402/1‑150 (MQ=255)
                                                                             tGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGc                                                                     >  2:740819/1‑149 (MQ=255)
                                                                              gacgacCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                   >  2:1022927/1‑150 (MQ=255)
                                                                               acaacCCGCCTTCCGATTGCGGAATCTCCCCGTTATCTTAATAAAGTAAATACCGCGAAAAAATCCTCCGGCCGCGGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTGTGCAAAt                                                                 <  2:980040/148‑1 (MQ=255)
                                                                                   ccccccATCCCTCTGGGGATTCTCGCCGTTATTTTTATAAAGGAAATCCCGCGCAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                             <  1:489898/151‑1 (MQ=255)
                                                                                     cGCCATCCCTCTGCGGAAGCTCGCCGTTATCTTTATAAATTGAATACCGCGCAAAAATCCTACCCCCGCCGATAATTCCATTACCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGaa                                                           <  2:289271/151‑1 (MQ=255)
                                                                                      cccACCCCTCTGGGAATTCCCGCCTTTATTTTTAAAAAGTAACTCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                                          <  2:767073/150‑1 (MQ=255)
                                                                                       ccATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATg                                                         >  1:1194604/1‑151 (MQ=255)
                                                                                       ccATCACTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                                          >  1:714461/1‑150 (MQ=255)
                                                                                         aTCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt                                                       >  1:644830/1‑151 (MQ=255)
                                                                                           cccTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                     <  1:60863/151‑1 (MQ=255)
                                                                                               ctGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCa                                                   <  1:258468/149‑1 (MQ=255)
                                                                                                tGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                                                >  1:181660/1‑151 (MQ=255)
                                                                                                tGCGGAATCTCCCCGTTATTTTTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCaa                                                  <  1:756707/149‑1 (MQ=255)
                                                                                                 gCGTAATCTCGCCTTTATTTTTAAAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatta                                               <  2:1375351/151‑1 (MQ=255)
                                                                                                 gCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                                                >  1:993328/1‑150 (MQ=255)
                                                                                                  cgaaaTCTCGCCGTTATTTTTATAAAGTGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatta                                               <  2:881064/147‑1 (MQ=255)
                                                                                                   ggatccccGCCTTTATCTTTAAAAAGTGAATACCCCGCAAAAACCCAACCGCCGCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattaa                                              <  2:761455/143‑2 (MQ=255)
                                                                                                   ggAATCTCGCCGTTTACTTTATAAAGTGAATCCCTGGCAAAAGCCCTCCCGCCGCCGATAATTCCACTAACCGCCCCTGACGATGCTCAGGGGCAAAAAGGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt                                                        <  2:404348/140‑1 (MQ=255)
                                                                                                   ggAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                                             <  1:56879/151‑1 (MQ=255)
                                                                                                   ggAATCTCGCCGTTATCTTTAAAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAatt                                                <  1:573602/148‑1 (MQ=255)
                                                                                                   gaaaTCTCGCCGTTATCTCTATAAAGTGAATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                              <  1:333133/148‑1 (MQ=255)
                                                                                                     aaTCTCGCCGTTATCTTTATAAAGTGAATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTAttt                                            <  1:662530/150‑1 (MQ=255)
                                                                                                      atttcGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGc                                          >  2:581081/4‑151 (MQ=255)
                                                                                                      atctgcccGTTATCTTTATAAAGAGAATCCCTGGAAAAAATCCTCCCGCAGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt                                             <  2:993417/142‑1 (MQ=255)
                                                                                                          cGCCTTTATTTTTATAAAGTAAATCCCGCCCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCca                                        <  1:1104827/149‑1 (MQ=255)
                                                                                                          cGCCGTTATCTTTATAAAGTGAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                       <  1:1282469/150‑1 (MQ=255)
                                                                                                           gCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCcac                                       >  1:832908/1‑149 (MQ=255)
                                                                                                                  tCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCAATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACa                               >  2:650578/1‑150 (MQ=255)
                                                                                                                    tttATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaa                             >  2:466305/1‑150 (MQ=255)
                                                                                                                    tttATAAAGTAAATACCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaa                             <  1:302693/150‑1 (MQ=255)
                                                                                                                    tttATAAAGAGAATACCGCGAAAAAATGCTACAGACGCCGATAATTCCATTAATAGCGCCTGAAGATGCACAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGAAATTATTTGCCACAGGTAACaaa                             >  2:291943/1‑150 (MQ=255)
                                                                                                                     tgatatAGTGAACCCCGCGCAAAAACCCTTCCGCCGCCTATAAGCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAAGTATTTGCCACAGGTAACaaaa                            <  2:490806/148‑1 (MQ=255)
                                                                                                                       ataAAGTGAATCCCGCGAAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaaa                           <  1:911114/149‑1 (MQ=255)
                                                                                                                        aaaaGTAAATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                        <  1:572262/150‑1 (MQ=255)
                                                                                                                         aaaGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                        >  1:827814/1‑150 (MQ=255)
                                                                                                                             tGAATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTc                    <  1:545196/150‑1 (MQ=255)
                                                                                                                             tGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc                   >  1:1269447/1‑151 (MQ=255)
                                                                                                                             cgAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGTCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACTTCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGACACAGGTAACAAAAAACCAGTc                    >  2:986497/2‑150 (MQ=255)
                                                                                                                               aaTCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGTGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcg                  <  1:110056/150‑1 (MQ=255)
                                                                                                                               aaTACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                 >  1:160947/1‑151 (MQ=255)
                                                                                                                                atcccGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                <  1:129989/148‑1 (MQ=255)
                                                                                                                                   ccGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                <  1:1310534/148‑1 (MQ=255)
                                                                                                                                   ccGCGCAAAAACCCCCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGaa              <  2:205194/150‑1 (MQ=255)
                                                                                                                                   ccGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa               <  1:575732/149‑1 (MQ=255)
                                                                                                                                     gcgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGtt           >  2:264211/1‑151 (MQ=255)
                                                                                                                                     ccgaaaaaaCCCTCCCCCCTCAAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGt            <  2:689569/146‑1 (MQ=255)
                                                                                                                                         aaaaaTCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGtt           <  1:4931/147‑1 (MQ=255)
                                                                                                                                         aaaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt        >  1:828456/1‑150 (MQ=255)
                                                                                                                                         aaaaaCCGTCCCGCGGCGGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa       <  2:1298816/151‑1 (MQ=255)
                                                                                                                                            aaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGCATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGa     >  2:316238/1‑150 (MQ=255)
                                                                                                                                               ccTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAAt   <  2:557794/149‑1 (MQ=255)
                                                                                                                                               ccTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATc  <  1:1116457/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCGCCGTTATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATC  >  NC_000913/3105353‑3105645

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