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

GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222085‑1222382
                                                                                                                                                    |                                                                                                                                                      
gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaa                                                                                                                                                     <  1:1156393/151‑1 (MQ=255)
    aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat                                                                                                                                                 <  2:1076070/151‑1 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                            >  1:702349/1‑151 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                            >  2:302317/1‑151 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                            <  1:934625/151‑1 (MQ=255)
             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                        <  1:1224208/151‑1 (MQ=255)
              aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAt                                                                                                                                       >  2:119175/1‑151 (MQ=255)
                    ggCGCGTAATTTGGGGGAAATTAAACCTGGCGTAAATGGCACCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTTAGCGCGCAGCTGGGTGATAATGGCCACAATGACACCGCAGCGATGGTGGGCCTGACACATAAGTTCCGATcc                                                                                                                                     >  2:498121/1‑147 (MQ=255)
                      cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGt                                                                                                                               >  2:74943/1‑151 (MQ=255)
                       gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  1:778334/1‑151 (MQ=255)
                       gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  2:1034022/1‑151 (MQ=255)
                        aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  1:488104/1‑150 (MQ=255)
                           aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCt                                                                                                                          >  1:1159521/1‑151 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  1:305277/1‑150 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  2:1313623/1‑150 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  2:998566/1‑150 (MQ=255)
                                     aaaTTAAACTGGGGGAAAATGGCAACCTGAACCCAGCGGCCGGCGTCGGGGGAATGTGGGGTTGCAGCTGGGTGATAATGGCTAAAATGACCCCGAAGTGATGGTGGGCTTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACcg                                                                                                                 <  2:697339/150‑1 (MQ=255)
                                      aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               >  1:496168/1‑151 (MQ=255)
                                         tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  1:518490/1‑151 (MQ=255)
                                            aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt                                                                                                         >  1:333066/1‑151 (MQ=255)
                                             cTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        >  1:544359/1‑151 (MQ=255)
                                                ggCGAAAAGGGCACCCTGAACCCAGCGGCCGGCGTGGGGGGAAGGGGGGGGTGCAGTGGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                          <  1:759172/146‑1 (MQ=255)
                                                ggCGAAAAGGCCACCCTAAACCCAGCGGCCAGGGTGGGGGGAAGGGGGGCGGGCAGGGGGGTAATAAGGGCTAAAATGACCCCGCGGTGAGGGGGGGCTTAAAATAAAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                                     <  1:556239/151‑1 (MQ=255)
                                                 gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    >  2:1232023/1‑151 (MQ=255)
                                                      aaggccacccTGAACCCAGCGGCCGGCGTCGGGGGAATGGGGGGGTGAGGCGGGGTGATAATGGCTAAAATGACCCCGCGGTGAGGGGGGCCCTAAAATATAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                               <  1:917309/143‑1 (MQ=255)
                                                      aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                >  2:1206302/1‑150 (MQ=255)
                                                           cAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  2:851977/1‑151 (MQ=255)
                                                            acccaaaaTACAGGGGCCAGCGTGTGGGGATTGGGGGGGTGCAGGTGGTTAAAAATGTCTACAACGACCCCCCAGTGATGGTGGGCCTAAAATAAATGTTCTCATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  2:211240/143‑1 (MQ=255)
                                                            aaCCTGAACCCAGGGCCGGGCGTGGGGGGAATGGGGGCGTGGGGCTGGGTGATAAGGGCTAAAATGACCCCGCGGTGAGGGTGGGCCTAAAATATATGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  2:39209/149‑1 (MQ=255)
                                                            aaCCTGAACCCAGCGGCCAGGGTCTGGGAATTGTGGGGGTGCAGCTGGTTGATAAGGGCTACAAGGACCCCGCGGTGGTGGGGGGCCTAAAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  1:613136/149‑1 (MQ=255)
                                                              ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        >  1:928169/1‑150 (MQ=255)
                                                              ccTAAACCCAGCGCCCGGCGTGTGGGGAATGGGGGCGTGCAGCTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCTTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  1:345268/151‑1 (MQ=255)
                                                               ctaaatcaggggcccAGCCTCGGGGGTTTGTGGGGGTGCAGTGGGGGTAAAATGGGCAAAATGCCCCCGCATGGTGGGGGGGCCTAAAATAAAAGTTCATACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  2:940798/137‑1 (MQ=255)
                                                                 gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                     >  2:934152/1‑150 (MQ=255)
                                                                  aaccaaGCGGCCGGCTTCGGGGGAATGGGGGCGGGCAGCGGGGTAATAATGCTAACAACAACCCGGCAGTGATGGGGGGCCGAAAATAAATTTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   <  2:378391/146‑1 (MQ=255)
                                                                    tCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  2:867751/1‑151 (MQ=255)
                                                                        gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGGTGGGTGATAAGGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                             <  1:99888/151‑1 (MQ=255)
                                                                          ggCCAGCGTCTGGGGATGGGGGGTGGGCGGTTGGGTGAAAATGGCTAAAAGGACCCGGCAGTGGTGGGGGGCTGAAAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATTGACGCGGTACATATGAAAtatta                                                                             <  1:342567/149‑1 (MQ=255)
                                                                          ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                           >  1:226249/1‑151 (MQ=255)
                                                                           gCCAGCGTCTGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            <  1:912598/149‑1 (MQ=255)
                                                                            ccAGCTTCTGGGGATGGTGGGCGTGCAGCGGGGTAATAATGGTTACAATCCCCCCCCAGTGAGGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  2:408732/150‑1 (MQ=255)
                                                                             cAGCGTGTGGGTAATGGGGGCGTGCAGCGGGGTGATAATGGCTACAAGGACCCCTCGGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                         <  2:226245/150‑1 (MQ=255)
                                                                               gCGTCTGGGGAATGTGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg                                                                      <  1:1322474/151‑1 (MQ=255)
                                                                                  tgtGGGGAGGGTGGGAGTGAGGGGGGGTAATATTGGCAAAATTAACACCCGAGGGATGGTGGGACTAAAAAATAGTTATTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAGTATTATTTTTTg                                                                      <  1:1213144/146‑1 (MQ=255)
                                                                                     ggggAATGGGGGGGTGCGGCGGGTTAATAATGTCTACAATGCCCCCGTAGGGATGGGGGCCTGAAAATATAAGTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                                <  2:837794/151‑1 (MQ=255)
                                                                                              ggggGTGCAGCGGGGTTATAATGCCTAAAATGACCCCGGAGTTGTGGTGGCTATGAAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc                                                            <  2:197919/146‑1 (MQ=255)
                                                                                              gggcgggcAGGGGGGTAAAAATGGCAACAAAGCCCCCGAGGTGAGGGGGGGCAAAAATTATAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg                                                       <  2:111538/151‑1 (MQ=255)
                                                                                                  gTGGAGCGGGGTGATAATGGTAAAAATGACCCCGGGGTGTGGGTGGGCGTAAAAAATAGTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt                                                     <  2:517840/149‑1 (MQ=255)
                                                                                                      ggCTGGGTGATAAGGGCTAAAAGGCCCCCGCGGTGAGGGGGGGCATAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca                                                 <  1:74943/148‑1 (MQ=255)
                                                                                                               aaaaagccTACATTGCCCCCGCAGTGTGGGGGGGCTTAAAAAATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                                      <  1:1018052/144‑1 (MQ=255)
                                                                                                                tAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  1:272596/1‑151 (MQ=255)
                                                                                                                 aaTGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  1:747197/1‑150 (MQ=255)
                                                                                                                    ggCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                                  >  1:1032391/1‑150 (MQ=255)
                                                                                                                        aCAATTACCCCCAAGTTATGGTGGGCCAAAAATATAATTTTTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                                 <  1:596674/147‑1 (MQ=255)
                                                                                                                         cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  1:1338223/1‑151 (MQ=255)
                                                                                                                                        gagggggggCTTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  <  1:451503/140‑1 (MQ=255)
                                                                                                                                           ggtggtCCTGAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               <  2:1111522/146‑1 (MQ=255)
                                                                                                                                           ggtggGCATAAAAAATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg            <  1:302323/149‑1 (MQ=255)
                                                                                                                                           ggggggCTTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc           <  1:119177/147‑1 (MQ=255)
                                                                                                                                                   tGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:1334169/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222085‑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: