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

TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG  >  NC_000913/1222086‑1222372
                                                                                                                                                     |                                                                                                                                           
tgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaa                                                                                                                                           >  2:1222753/1‑150 (MQ=255)
    aGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                      <  3:640397/151‑1 (MQ=255)
    aGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                      <  2:959300/151‑1 (MQ=255)
      tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                      >  2:339318/1‑149 (MQ=255)
       cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                   <  2:906695/151‑1 (MQ=255)
           cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                               <  3:367101/151‑1 (MQ=255)
             aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAt                                                                                                                             >  3:1634301/1‑151 (MQ=255)
                tGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATccc                                                                                                                          >  3:1178365/1‑151 (MQ=255)
                 ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccg                                                                                                                         <  3:578015/151‑1 (MQ=255)
                  gggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgc                                                                                                                        <  2:516907/151‑1 (MQ=255)
                   ggCGAGTAATTTGGGGGAAATTAAACTGGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgcc                                                                                                                       <  3:1493266/151‑1 (MQ=255)
                      gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                    >  2:1557937/1‑151 (MQ=255)
                       aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                    >  4:10544/1‑150 (MQ=255)
                              tGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaa                                                                                                             >  2:1018998/1‑150 (MQ=255)
                                ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaaa                                                                                                           >  2:691468/1‑150 (MQ=255)
                                 gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAAc                                                                                                         >  3:307116/1‑151 (MQ=255)
                                    aaaTTAATCTTGGGGTAAGTGGCACCCTGAACCCAGCGGCCAGCGTGGGGGGAAGGGGGGCGTGCAGTGGGGTGATAATGTCTACAATGCCCCCGCGGGGATGGGGGCCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAAcc                                                                                                        <  3:304113/149‑1 (MQ=255)
                                    aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgc                                                                                                      >  2:1455368/1‑151 (MQ=255)
                                       ttAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                   >  3:696290/1‑151 (MQ=255)
                                       ttAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                   >  2:275108/1‑151 (MQ=255)
                                              tGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                             >  3:333600/1‑150 (MQ=255)
                                               ggggaaaaTGGCACCCTGAATCCAGCGGCCAGCGTCTGGGGAATGGGGGCGGGCAGCGGGTTAATAATGGCTACAATGACCCCGAGGTGAGGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                <  2:1217547/141‑1 (MQ=255)
                                               ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                           >  3:1093777/1‑151 (MQ=255)
                                                gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                          >  3:1590041/1‑151 (MQ=255)
                                                    aaaGGGCAACCTAAATCCAGCGGCCGGCGTCGGGGGAATGGGGGCGTGCGGCGGGGTGATAAGGGCTACAATGACCCCGCAGTGATGGTGGGCCGAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                      <  3:1157531/151‑1 (MQ=255)
                                                        ggCAGCCTGAACCCGGCGGCCAGCGTCTGGTGAGGGGGGGGGTGCGGGGGGGTGATAAGGGCTACAATGACCCCGCGTTGATGTTGGCCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                     <  2:1557694/148‑1 (MQ=255)
                                                             ccggAACCCGGCGGCCGGCTTCTGGGGAAGGGGGGGGTGCGGGGGGGTGATAATGGCTAAAATGACCCCGCGGTGTGGGGGGGCCTAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                              <  3:850427/147‑1 (MQ=255)
                                                             ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATAGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                              >  2:767702/1‑150 (MQ=255)
                                                             ccTGAACCCAGCGGCCGCCCTGGGGGGAGGGGGGGGGGGCAGTGGGGTGATAATGCCTAAAAAAACGCCAAAGTGGTGGGGGGTCTGAATTATAAGTCCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                               <  2:564589/149‑1 (MQ=255)
                                                              cTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGCGATAATGGCTACAATGACACCGCAGTGATGGCGGGCCTGAAATATAAGCTCTGATCCCGCCGCTAGCTAAAAAACCGCGTCGGATTCATCGACGCGGTACctat                                                                            >  3:1608280/1‑151 (MQ=255)
                                                              cTGAACCCGGCGGCCAGCGTCGGGGGAATGGGGGGGTGCAGTGGGGTGATAATGGCAAAAATGCCCCCGCAGGGTTGGTGGTCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                              <  2:855883/149‑1 (MQ=255)
                                                                gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                           >  2:1267860/1‑150 (MQ=255)
                                                                 aacccggCGGCCGGCGTCGGGGAATGGGGGGGGTGCGGCTGGGTAAAAATGGCTACAATGCCACCGAGGTGGTGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                            <  2:1187708/142‑1 (MQ=255)
                                                                 aaTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                          >  2:792444/1‑150 (MQ=255)
                                                                  aTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                         >  2:515534/1‑150 (MQ=255)
                                                                  aTCCAGCGGCCAGCGTCGGGGGAATGTGGGCGGGCGGCTGGGTGAAAATGCCTACAATGACCCCGCGGTGAGGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa                                                                        <  2:520829/151‑1 (MQ=255)
                                                                       gCGGCCGGCGTCTGGGGAAGGTGGGCGTGCAGCTGGGTAATAATGGCTAAAATGCCCCCGAGGTGATGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                   <  2:425802/151‑1 (MQ=255)
                                                                       gCGGCCGGCGTCGGGGGAAGGTGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                       <  2:1328728/147‑1 (MQ=255)
                                                                       gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                    >  3:802383/1‑150 (MQ=255)
                                                                        cggccggcGTCGGGGGAATGTGGGCGTGCAGTTGGGTGATATTGGCTACAATGACCCCGCAGTGTTGGTGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                    <  2:290092/149‑1 (MQ=255)
                                                                         ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTTATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                     <  3:1163038/147‑1 (MQ=255)
                                                                         ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                 >  2:1478196/1‑151 (MQ=255)
                                                                         gcccggCGTCGGGGGAAGGTGGGCGTGAAGGTGGGTGATAATGGCTACAATGACCCCGAGGTGATGGTGGGCTTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                  <  3:692377/145‑1 (MQ=255)
                                                                          ggccggcgtgggggAATGGGGGGCGGGCAGCGGGTTGAAAAGGGCTAAAATGCCCCCGCAGGGGTGGTGGGCTAAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                      <  3:89444/136‑1 (MQ=255)
                                                                              gCGTCGGGGGAATGGGGGCGTGCAGCTGGGTAATAATGGCTAAAATGACACCGCAGTGATGGGGGGCATAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                              <  3:686064/149‑1 (MQ=255)
                                                                                 tCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                         >  2:1505457/1‑151 (MQ=255)
                                                                                    ggggAGTGTGGGGGTGCAGGGGGGTAATAATGGCAAAAAGACCACCGAGGGGGGGGTGGGCTTAAAAAAAAAGTTCAAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                       <  3:401390/150‑1 (MQ=255)
                                                                                    ggggAAGTGTGGGCGGGCAGCTGGGTGAAAAGGCAAACAACCCCCCCATAGTGATGGGGGGCCTAAATTATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                       <  3:1441937/151‑1 (MQ=255)
                                                                                    ggggAATGGGGGCGTGCGGCTGGGTGAAAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                      <  2:723590/151‑1 (MQ=255)
                                                                                    ggggAATGGGGGCGGGCAGCTGGGTGAAAAGGTCTAAAAGGCCCCCGAGGTGAGGGGGGCCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                      <  2:116528/151‑1 (MQ=255)
                                                                                    ggggAATGGGGGCGGGCAGCGGGGTGATAATGGCTACAATGACCCCGCGGTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                         <  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)
                                                                                                                                                     |                                                                                                                                           
TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG  >  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: