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

TGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222084‑1222382
                                                                                                                                                       |                                                                                                                                                     
tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑g                                                                                                                                                       >  1:364729/1‑150 (MQ=255)
tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑g                                                                                                                                                       <  1:119522/150‑1 (MQ=255)
 gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑Ga                                                                                                                                                      >  1:599503/1‑150 (MQ=255)
     aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAatat                                                                                                                                                 <  2:1162032/151‑1 (MQ=255)
       ggcGAACAAGATGGAGCGAGTAAGTTTGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCACCGTATGGGGAATGTGGGCATGCAGCTAGGAGATAAAGGCAACATTGTCA‑CCGCAGTGATGGTGGGCCT‑GAAATATaa                                                                                                                                               <  2:1122117/149‑1 (MQ=255)
          gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGtt                                                                                                                                            >  1:17584/1‑151 (MQ=255)
           aaCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTc                                                                                                                                           <  1:683764/151‑1 (MQ=255)
            aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCt                                                                                                                                          <  1:287968/151‑1 (MQ=255)
             cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTg                                                                                                                                         >  1:9696/1‑151 (MQ=255)
              aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTg                                                                                                                                         >  2:203035/1‑150 (MQ=255)
               aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGa                                                                                                                                        >  2:1464644/1‑150 (MQ=255)
                 aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGGGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATcc                                                                                                                                     >  1:1458688/1‑151 (MQ=255)
                   ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATccc                                                                                                                                    <  1:402440/150‑1 (MQ=255)
                    gggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCccgc                                                                                                                                  >  1:138716/1‑151 (MQ=255)
                     ggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCccgcc                                                                                                                                 <  1:952759/151‑1 (MQ=255)
                        gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  1:1174915/1‑151 (MQ=255)
                         aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  2:356287/1‑150 (MQ=255)
                               ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  2:537908/1‑150 (MQ=255)
                                 ggggaAAATTAAACTTGGCGAAAATGGCAACCTGAATCCGGCGGCCAGCGTCGGGGGAATGTGGGCGTGCAGCTGGTTAATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑AAAATATAAGTTCGGATCCCGCCGTTAGCTaaa                                                                                                                       <  1:520328/149‑1 (MQ=255)
                                     gAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTTTGGGGAATGTGGGCGCGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcg                                                                                                                 >  2:175368/1‑151 (MQ=255)
                                      aaaTAAAACTGGGCGAAAGGGGCACCCTGAACCCAGGGGCCAGCGTCGGGGGAAGGTGGGGGTGCAGCGGGGGGATAAGGGCTACAAGGACC‑CCGAAGTGATGGGGGGCCT‑AAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACcg                                                                                                                 <  1:910994/150‑1 (MQ=255)
                                       aaTAAAACTGGGCGAAAAGGCAAACCTAAACCCAGCGGCCAGCGTGGGGGGATTGGGGGGGGGCGGCTGGTTAATAAGGGCTAAAATGCCC‑CCGCAGTGAGGGTGGGCGT‑AAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               <  2:88079/151‑1 (MQ=255)
                                           aaaCTTGGGGAAAATGGCCCCCGGAACCCAGGGGCCGGCGTCGGGGGAATGGGGGGGTGAGGCGGGGTAAAAATGGCTACAAGGCCC‑CCGTGGTGAGGGGGGCCAT‑GAAATATAATTTCTTTCCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            <  2:271212/150‑1 (MQ=255)
                                           aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  2:91963/1‑150 (MQ=255)
                                           aaaCTGCGCGAAAAGGGCACCCTGATCCCGGCGGCCGGCGTGTGGGGAAGGGGGGGGGGCAGGTGGGTGAAAAGGCCTAAAATGACC‑CGGCAGGGTTGGTGGCCTT‑AAAATATAATTTCTGTTCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                           <  1:1074878/151‑1 (MQ=255)
                                            aaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                           >  1:613913/1‑150 (MQ=255)
                                                tGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                       >  2:678739/1‑150 (MQ=255)
                                                 ggCGAAAAGGGCACCCGGAACCCGGCGGCCAGCGTCGGGGGAAGGTGGGCGTGAAGTTGGGTGATAATGGCTACAAGGACA‑CCGCAGTGAGGGTGGGCCT‑GAAATATAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTAt                                                                                                         <  1:48361/147‑1 (MQ=255)
                                                  gCGTAAAGGGCACCCTGAACCCAGCGGCCAGCGTCGGGGGATTGGGGGCGTGCAGCGGGGTAATAATGGCTAAAATGACA‑CCGCAGTGATGGTGGGCTT‑AAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  1:407378/151‑1 (MQ=255)
                                                   cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  1:721248/1‑151 (MQ=255)
                                                    gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  1:1407108/1‑150 (MQ=255)
                                                      aaaggccacccTAAACCCGGCGCCCAGCGTCTGGGGAATGTGGGGGTGAAGGTGGGTAATAATGTCTACAATGACA‑CCGCAGTGATGGGGGGCCG‑AAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                 <  1:139550/141‑1 (MQ=255)
                                                      aaaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                >  1:645344/1‑151 (MQ=255)
                                                      aaaGGGCAACCTAAATCCAGCGGCCAGCGTGGGGGGAGGGGGGGGGGGCAGGGGGGTGATAAGGGCTAAAAGACCC‑CCGGAGTGGGGGGGGGCAG‑AAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   <  2:910651/148‑1 (MQ=255)
                                                       aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                               >  2:303165/1‑151 (MQ=255)
                                                         gggCAACCTGAACCCGGCGGCCAGCGTCGGGGGAATGTGGGGGGGCAGGGGGGTGATAGGGGCTACAAGGCCC‑CCGCAGTGATGGTGGGCTT‑AAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                             <  2:608596/150‑1 (MQ=255)
                                                           gCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                            >  1:693574/1‑150 (MQ=255)
                                                             acccTGACTCCGGCGGCCGGCGTCGGGGAAATGGGGGCGGGCAGCTGGGTAAAAATGGCTAAAATGACC‑CCGCAGTGAGGGTGGGCCT‑AAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                         <  1:777786/149‑1 (MQ=255)
                                                             acccTAAACCCAGCGGCCAGCGTCGGGGAAAGGTGGGCGGGCAGGTGGGTAATAAGGGCTACAATGCCC‑CCGCAGTGAGGGGGGCCCT‑AAAATAAAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          <  1:533128/148‑1 (MQ=255)
                                                               cctaacgccAGCGGCCGGGTTGGGGGGAAGGTGGGCGTGCATATGGGTAATAATGGCTACAGTGCCC‑CCGATGTGAGGGGGGCCAT‑AAAAAATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                         <  2:168125/142‑1 (MQ=255)
                                                               ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        >  1:439503/1‑150 (MQ=255)
                                                                  gAACCCGGCGGCCGGCGTCGGGGGAAGGGGGGCGTGCAGCTGGGTAATAATGGCTACAATGACC‑CCGCAGTGGTGGTGGGCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                      <  1:417248/149‑1 (MQ=255)
                                                                     tCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  1:866871/1‑151 (MQ=255)
                                                                      ccggCGGCCGGCGTCGGGGGAAGGGGGGGGTGAGGGTGGGTGATAATGGCTACAATGACC‑CCGCGGTGGGGGTGGGCAT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 <  1:358067/147‑1 (MQ=255)
                                                                      ccAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  1:778238/1‑150 (MQ=255)
                                                                        aGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              >  1:323744/1‑151 (MQ=255)
                                                                         gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              >  2:916637/1‑150 (MQ=255)
                                                                          cgcacgtcGTCGGGGGAGTGTGGGGGTGCGGCGGGGGAAAAGTGGATAAAAGGACC‑CCGCGGGAAGGGTGGGACT‑AAAAAAAATGTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  2:1431901/140‑1 (MQ=255)
                                                                           ggCCAGCGTCGGGGGAAGGTGGGCGTGCGGCGGGGTGATAATGGCTAAAATGACC‑CCGCAGTGATGGTGGCCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              <  1:438032/148‑1 (MQ=255)
                                                                             ccggCGTCGGGGGAATGTGGGGGTGCGGTGGGGTAATAATGCCAACAATACCC‑CCGATGGGATGG‑GGGCCTAAAAATATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattattc                                                                          <  2:1101090/147‑2 (MQ=255)
                                                                                  gtttGGGTAGTGTGGGGGTGCGGCGGGGTAATAGTGGCTAAAATGACC‑CCGCAGTGTTGGTGGGCCT‑AAAATAAAGTTTCTGTACCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  1:4706/142‑1 (MQ=255)
                                                                                  gTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    >  1:643840/1‑151 (MQ=255)
                                                                                  gTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    >  2:537554/1‑151 (MQ=255)
                                                                                   ttggggaaagggggggGTGCAGTGGGGTAATAATGGCTACAAGAACC‑CCGCAGGGAGGGGGGCCAT‑AAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    <  2:1222790/138‑1 (MQ=255)
                                                                                   tcgggggATTGGGGGCGGGCAGCGGGGTGATAATGGCTACAAGAACC‑CCGCAGTGATGGTGGGCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       <  1:160076/144‑1 (MQ=255)
                                                                                   tcgggggAATGTGGGCGTGGAGCTGGTTGAAAATTCAAAAAAGGCCC‑CCGCAGTGAGGGTGGGCTT‑AAAAAATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                     <  1:115300/146‑1 (MQ=255)
                                                                                      ggggAATGGGGGCGTGCAGCGGGGTGATAATGGCTAAAATGACC‑CCGCAGTGATGGTGGGCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                                <  1:642203/151‑1 (MQ=255)
                                                                                         gAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc                                                              >  1:318740/1‑150 (MQ=255)
                                                                                                     gCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtc                                                  >  1:567708/1‑150 (MQ=255)
                                                                                                       aGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCAATAGCACGCATGGTGtcatc                                               >  2:1408010/1‑151 (MQ=255)
                                                                                                        gttGGGTGATAATGGCTAAAATGACC‑CCGCGGTGATGGTGGGCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                    <  1:535625/143‑1 (MQ=255)
                                                                                                           ggttGATAAGGGCTAC‑ATTACACCCGCAGGGGTGGTGTGGCT‑AAAATATAATTTCGGTCACCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                            <  2:537219/147‑1 (MQ=255)
                                                                                                           gggTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt                                           >  2:1127166/1‑151 (MQ=255)
                                                                                                             gTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         >  1:390591/1‑151 (MQ=255)
                                                                                                                aTAATGGTTAAAATGACC‑CCGCGGTGAGGGTGGGCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         <  1:962444/148‑1 (MQ=255)
                                                                                                                  aaTGGCTACAATGCCA‑CCGAAGTGAGGGTGGGCCT‑AAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                                    <  1:1399972/151‑1 (MQ=255)
                                                                                                                     ggCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                                  >  2:1350132/1‑150 (MQ=255)
                                                                                                                          cAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  2:553790/1‑151 (MQ=255)
                                                                                                                           aaTGACA‑CCGGAGTGGGGGTGGGCTT‑AAAATATAAGTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            <  1:642631/150‑1 (MQ=255)
                                                                                                                            aTGACA‑CCGCAGTGATGGGGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATATCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                          >  2:93172/1‑151 (MQ=255)
                                                                                                                                ca‑cCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      >  2:1068934/1‑151 (MQ=255)
                                                                                                                                   ccGGAGTGATGGTGGCCCT‑AAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                       <  1:158821/148‑1 (MQ=255)
                                                                                                                                    cGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                   >  1:371448/1‑151 (MQ=255)
                                                                                                                                     gaaGTGAGGGAGGGCCA‑AAAATATAATTTGTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  <  1:1138365/149‑1 (MQ=255)
                                                                                                                                     gCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  >  1:1075794/1‑151 (MQ=255)
                                                                                                                                         tGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt              >  1:441754/1‑151 (MQ=255)
                                                                                                                                          gATGGTGGGCAT‑AAAAAAAAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      <  1:185085/142‑1 (MQ=255)
                                                                                                                                             ggtggGCCT‑AAAATATAATTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               <  2:230703/146‑1 (MQ=255)
                                                                                                                                             ggggggCCT‑AAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg            <  1:1434743/146‑1 (MQ=255)
                                                                                                                                                gggCTT‑AAAATATAAATTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa             <  1:1362287/145‑1 (MQ=255)
                                                                                                                                                gggCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa        >  2:952018/1‑150 (MQ=255)
                                                                                                                                                  gCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc     >  2:193269/1‑151 (MQ=255)
                                                                                                                                                    cT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCt   >  2:233283/1‑151 (MQ=255)
                                                                                                                                                     t‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:1376751/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222084‑1222382

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