Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA44.4% 105.1 / 34.7 81pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (37/8);  new base A (0/36);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.42e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG  >  NC_000913/1222086‑1222372
                                                                                                                                                     |                                                                                                                                           
tgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaa                                                                                                                                           >  2:1222753/1‑150 (MQ=255)
    aGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                      <  3:640397/151‑1 (MQ=255)
    aGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                      <  2:959300/151‑1 (MQ=255)
      tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                      >  2:339318/1‑149 (MQ=255)
       cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                   <  2:906695/151‑1 (MQ=255)
           cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                               <  3:367101/151‑1 (MQ=255)
             aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAt                                                                                                                             >  3:1634301/1‑151 (MQ=255)
                tGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATccc                                                                                                                          >  3:1178365/1‑151 (MQ=255)
                 ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccg                                                                                                                         <  3:578015/151‑1 (MQ=255)
                  gggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgc                                                                                                                        <  2:516907/151‑1 (MQ=255)
                   ggCGAGTAATTTGGGGGAAATTAAACTGGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgcc                                                                                                                       <  3:1493266/151‑1 (MQ=255)
                      gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                    >  2:1557937/1‑151 (MQ=255)
                       aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                    >  4:10544/1‑150 (MQ=255)
                              tGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaa                                                                                                             >  2:1018998/1‑150 (MQ=255)
                                ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaaa                                                                                                           >  2:691468/1‑150 (MQ=255)
                                 gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAAc                                                                                                         >  3:307116/1‑151 (MQ=255)
                                    aaaTTAATCTTGGGGTAAGTGGCACCCTGAACCCAGCGGCCAGCGTGGGGGGAA‑G‑GGGGGCGTGCAGTGGGGTGATAATGTCTACAATGCCCCCGCGGGGATGGGGGCCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAAcc                                                                                                        <  3:304113/149‑1 (MQ=255)
                                    aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgc                                                                                                      >  2:1455368/1‑151 (MQ=255)
                                       ttAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                   >  3:696290/1‑151 (MQ=255)
                                       ttAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                   >  2:275108/1‑151 (MQ=255)
                                              tGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                             >  3:333600/1‑150 (MQ=255)
                                               ggggaaaaTGGCACCCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GGGGGCGGGCAGCGGGTTAATAATGGCTACAATGACCCCGAGGTGAGGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                <  2:1217547/141‑1 (MQ=255)
                                               ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                           >  3:1093777/1‑151 (MQ=255)
                                                gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                          >  3:1590041/1‑151 (MQ=255)
                                                    aaaGGGCAACCTAAATCCAGCGGCCGGCGTCGGGGGAA‑T‑GGGGGCGTGCGGCGGGGTGATAAGGGCTACAATGACCCCGCAGTGATGGTGGGCCGAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                      <  3:1157531/151‑1 (MQ=255)
                                                        ggCAGCCTGAACCCGGCGGCCAGCGTCTGGTGAG‑G‑GGGGGGGTGCGGGGGGGTGATAAGGGCTACAATGACCCCGCGTTGATGTTGGCCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                     <  2:1557694/148‑1 (MQ=255)
                                                             ccggAACCCGGCGGCCGGCTTCTGGGGAA‑G‑GGGGGGGTGCGGGGGGGTGATAATGGCTAAAATGACCCCGCGGTGTGGGGGGGCCTAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                              <  3:850427/147‑1 (MQ=255)
                                                             ccTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATAGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                              >  2:767702/1‑150 (MQ=255)
                                                             ccTGAACCCAGCGGCCGCCCTGGGGGGAG‑G‑GGGGGGGGGCAGTGGGGTGATAATGCCTAAAAAAACGCCAAAGTGGTGGGGGGTCTGAATTATAAGTCCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                               <  2:564589/149‑1 (MQ=255)
                                                              cTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGCGATAATGGCTACAATGACACCGCAGTGATGGCGGGCCTGAAATATAAGCTCTGATCCCGCCGCTAGCTAAAAAACCGCGTCGGATTCATCGACGCGGTACctat                                                                            >  3:1608280/1‑151 (MQ=255)
                                                              cTGAACCCGGCGGCCAGCGTCGGGGGAA‑T‑GGGGGGGTGCAGTGGGGTGATAATGGCAAAAATGCCCCCGCAGGGTTGGTGGTCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                              <  2:855883/149‑1 (MQ=255)
                                                                gAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                           >  2:1267860/1‑150 (MQ=255)
                                                                 aacccggCGGCCGGCGTC‑GGGGAA‑TGGGGGGGGTGCGGCTGGGTAAAAATGGCTACAATGCCACCGAGGTGGTGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                            <  2:1187708/142‑1 (MQ=255)
                                                                 aaTCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                          >  2:792444/1‑150 (MQ=255)
                                                                  aTCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                         >  2:515534/1‑150 (MQ=255)
                                                                  aTCCAGCGGCCAGCGTCGGGGGAA‑T‑GTGGGCGGGCGGCTGGGTGAAAATGCCTACAATGACCCCGCGGTGAGGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa                                                                        <  2:520829/151‑1 (MQ=255)
                                                                       gCGGCCGGCGTCTGGGGAA‑G‑GTGGGCGTGCAGCTGGGTAATAATGGCTAAAATGCCCCCGAGGTGATGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                   <  2:425802/151‑1 (MQ=255)
                                                                       gCGGCCGGCGTCGGGGGAA‑G‑GTGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                       <  2:1328728/147‑1 (MQ=255)
                                                                       gCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                    >  3:802383/1‑150 (MQ=255)
                                                                        cggccggcGTCGGGGGAA‑T‑GTGGGCGTGCAGTTGGGTGATATTGGCTACAATGACCCCGCAGTGTTGGTGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                    <  2:290092/149‑1 (MQ=255)
                                                                         ggCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTTATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                     <  3:1163038/147‑1 (MQ=255)
                                                                         ggCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                 >  2:1478196/1‑151 (MQ=255)
                                                                         gcccggCGTCGGGGGAA‑G‑GTGGGCGTGAAGGTGGGTGATAATGGCTACAATGACCCCGAGGTGATGGTGGGCTTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                  <  3:692377/145‑1 (MQ=255)
                                                                          ggccggcgtgggggAA‑TGGGGGGCGGGCAGCGGGTTGAAAAGGGCTAAAATGCCCCCGCAGGGGTGGTGGGCTAAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                      <  3:89444/136‑1 (MQ=255)
                                                                              gCGTCGGGGGAA‑T‑GGGGGCGTGCAGCTGGGTAATAATGGCTAAAATGACACCGCAGTGATGGGGGGCATAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                              <  3:686064/149‑1 (MQ=255)
                                                                                 tCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                         >  2:1505457/1‑151 (MQ=255)
                                                                                    ggggAG‑T‑GTGGGGGTGCAGGGGGGTAATAATGGCAAAAAGACCACCGAGGGGGGGGTGGGCTTAAAAAAAAAGTTCAAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                       <  3:401390/150‑1 (MQ=255)
                                                                                    ggggAAGT‑GTGGGCGGGCAGCTGGGTGAAAAGGCAAACAACCCCCCCATAGTGATGGGGGGCCTAAATTATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                       <  3:1441937/151‑1 (MQ=255)
                                                                                    ggggAA‑T‑GGGGGCGTGCGGCTGGGTGAAAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                      <  2:723590/151‑1 (MQ=255)
                                                                                    ggggAA‑T‑GGGGGCGGGCAGCTGGGTGAAAAGGTCTAAAAGGCCCCCGAGGTGAGGGGGGCCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                      <  2:116528/151‑1 (MQ=255)
                                                                                    ggggAA‑T‑GGGGGCGGGCAGCGGGGTGATAATGGCTACAATGACCCCGCGGTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                         <  2:1133455/148‑1 (MQ=255)
                                                                                             gtgGGCGTGCAGCTGGGTGATAATGGCAAAAATGACCCCGCAGTGATGGGGGCCATAAAATATAATTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc                                                  <  2:988674/148‑1 (MQ=255)
                                                                                             gggggCGTGCAGTGGGTTGATAATGGCTACAATGACCCCGTGGTGATGGGGGCCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc                                                    <  3:241423/144‑1 (MQ=255)
                                                                                               gggCGTGCAGGGGGGTAATAATGGCTAAAATGCCCCCGAGGGGGTGGGGGGCCGAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                              <  2:485253/150‑1 (MQ=255)
                                                                                                 gCGTGCGGCTGGGTAATAAGGTCTAAAATGCCCCCGCAGTGAGGGTGGCCTAAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt                                           <  3:1037554/151‑1 (MQ=255)
                                                                                                      cAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcat                                      >  3:610159/1‑151 (MQ=255)
                                                                                                      cAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcat                                      >  2:422161/1‑151 (MQ=255)
                                                                                                       aGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                     >  2:373380/1‑151 (MQ=255)
                                                                                                           gggTGATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt                                 <  2:1186791/151‑1 (MQ=255)
                                                                                                            ggTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                >  2:636872/1‑151 (MQ=255)
                                                                                                              tGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGcc                              >  3:346642/1‑151 (MQ=255)
                                                                                                                aTAATGCCTAAAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                            <  2:309310/151‑1 (MQ=255)
                                                                                                                  aagggCTAAAATGACCCCGCAGTGGGGGGGGCCCTAAAATATAATTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                            <  2:1175761/146‑1 (MQ=255)
                                                                                                                     ggCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                        <  2:277786/150‑1 (MQ=255)
                                                                                                                        tACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                     >  2:253326/1‑150 (MQ=255)
                                                                                                                             tGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt               >  3:391093/1‑151 (MQ=255)
                                                                                                                               acacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg             >  2:1397225/1‑151 (MQ=255)
                                                                                                                                cacCGCGGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                <  2:784992/147‑1 (MQ=255)
                                                                                                                                cacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc            >  2:1363716/1‑151 (MQ=255)
                                                                                                                                 acCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg           >  2:1445036/1‑151 (MQ=255)
                                                                                                                                  ccgatgtgggggtggGCTTAAAATAAAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg          <  2:301100/141‑1 (MQ=255)
                                                                                                                                  ccGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg          >  3:1526757/1‑151 (MQ=255)
                                                                                                                                  ccGCAGTGAGGGGGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg          <  2:97355/151‑1 (MQ=255)
                                                                                                                                   cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg         >  3:1417929/1‑151 (MQ=255)
                                                                                                                                   cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg         >  2:1587068/1‑151 (MQ=255)
                                                                                                                                       gtgggggtggGCCTAAAATATAATTTCGAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt        <  2:1630785/143‑1 (MQ=255)
                                                                                                                                        tGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt    >  2:759849/1‑151 (MQ=255)
                                                                                                                                         gATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt    <  2:1444031/150‑1 (MQ=255)
                                                                                                                                            ggtggGCCTAAAATATATTTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt     <  2:1120766/146‑1 (MQ=255)
                                                                                                                                            ggggggCCTAAAATATATTTCCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa   <  2:1335312/145‑1 (MQ=255)
                                                                                                                                            ggggggCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg  <  2:129555/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑T‑GTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG  >  NC_000913/1222086‑1222372

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