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

GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222085‑1222382
                                                                                                                                                    |                                                                                                                                                      
gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑aa                                                                                                                                                     <  1:1156393/151‑1 (MQ=255)
    aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAatat                                                                                                                                                 <  2:1076070/151‑1 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt                                                                                                                                            >  1:702349/1‑151 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt                                                                                                                                            >  2:302317/1‑151 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt                                                                                                                                            <  1:934625/151‑1 (MQ=255)
             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGa                                                                                                                                        <  1:1224208/151‑1 (MQ=255)
              aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGAt                                                                                                                                       >  2:119175/1‑151 (MQ=255)
                    ggCGCGTAATTTGGGGGAAATTAAACCTGGCGTAAATGGCACCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTTAGCGCGCAGCTGGGTGATAATGGCCACAATGACACCGCAGCGATGGTGGGCCTG‑ACACATAAGTTCCGATcc                                                                                                                                     >  2:498121/1‑147 (MQ=255)
                      cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGt                                                                                                                               >  2:74943/1‑151 (MQ=255)
                       gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  1:778334/1‑151 (MQ=255)
                       gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  2:1034022/1‑151 (MQ=255)
                        aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  1:488104/1‑150 (MQ=255)
                           aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCt                                                                                                                          >  1:1159521/1‑151 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  1:305277/1‑150 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  2:1313623/1‑150 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  2:998566/1‑150 (MQ=255)
                                     aaaTTAAACTGGGGGAAAATGGCAACCTGAACCCAGCGGCCGGCGTCGGGGGAATGTGGGGTTGCAGCTGGGTGATAATGGCTAAAATGACCCCGAAGTGATGGTGGGCTTA‑AAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACcg                                                                                                                 <  2:697339/150‑1 (MQ=255)
                                      aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               >  1:496168/1‑151 (MQ=255)
                                         tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  1:518490/1‑151 (MQ=255)
                                            aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt                                                                                                         >  1:333066/1‑151 (MQ=255)
                                             cTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        >  1:544359/1‑151 (MQ=255)
                                                ggCGAAAAGGGCACCCTGAACCCAGCGGCCGGCGTGGGGGGAAGGGGGGGGTGCAGTGGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCCTA‑AAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                          <  1:759172/146‑1 (MQ=255)
                                                ggCGAAAAGGCCACCCTAAACCCAGCGGCCAGGGTGGGGGGAAGGGGGGCGGGCAGGGGGGTAATAAGGGCTAAAATGACCCCGCGGTGAGGGGGGGCTTA‑AAATAAAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                                     <  1:556239/151‑1 (MQ=255)
                                                 gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    >  2:1232023/1‑151 (MQ=255)
                                                      aaggccacccTGAACCCAGCGGCCGGCGTCGGGGGAATGGGGGGGTGAGGCGGGGTGATAATGGCTAAAATGACCCCGCGGTGAGGGGGGCCCTA‑AAATATAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                               <  1:917309/143‑1 (MQ=255)
                                                      aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                >  2:1206302/1‑150 (MQ=255)
                                                           cAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  2:851977/1‑151 (MQ=255)
                                                            acccaaaaTACAGGGGCCAGCGTGTGGGGATTGGGGGGGTGCAGGTGGTTAAAAATGTCTACAACGACCCCCCAGTGATGGTGGGCCTA‑AAATAAATGTTCTCATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  2:211240/143‑1 (MQ=255)
                                                            aaCCTGAACCCAGGGCCGGGCGTGGGGGGAATGGGGGCGTGGGGCTGGGTGATAAGGGCTAAAATGACCCCGCGGTGAGGGTGGGCCTA‑AAATATATGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  2:39209/149‑1 (MQ=255)
                                                            aaCCTGAACCCAGCGGCCAGGGTCTGGGAATTGTGGGGGTGCAGCTGGTTGATAAGGGCTACAAGGACCCCGCGGTGGTGGGGGGCCTA‑AAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  1:613136/149‑1 (MQ=255)
                                                              ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        >  1:928169/1‑150 (MQ=255)
                                                              ccTAAACCCAGCGCCCGGCGTGTGGGGAATGGGGGCGTGCAGCTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCTTA‑AAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  1:345268/151‑1 (MQ=255)
                                                               ctaaatcaggggcccAGCCTCGGGGGTTTGTGGGGGTGCAGTGGGGGTAAAATGGGCAAAATGCCCCCGCATGGTGGGGGGGCCTA‑AAATAAAAGTTCATACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  2:940798/137‑1 (MQ=255)
                                                                 gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                     >  2:934152/1‑150 (MQ=255)
                                                                  aaccaaGCGGCCGGCTTCGGGGGAATGGGGGCGGGCAGCGGGGTAATAATGCTAACAACAACCCGGCAGTGATGGGGGGCCGA‑AAATAAATTTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   <  2:378391/146‑1 (MQ=255)
                                                                    tCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  2:867751/1‑151 (MQ=255)
                                                                        gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGGTGGGTGATAAGGGCTACAATGACCCCGCAGTGATGGTGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                             <  1:99888/151‑1 (MQ=255)
                                                                          ggCCAGCGTCTGGGGATGGGGGGTGGGCGGTTGGGTGAAAATGGCTAAAAGGACCCGGCAGTGGTGGGGGGCTGA‑AAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATTGACGCGGTACATATGAAAtatta                                                                             <  1:342567/149‑1 (MQ=255)
                                                                          ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                           >  1:226249/1‑151 (MQ=255)
                                                                           gCCAGCGTCTGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCTTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            <  1:912598/149‑1 (MQ=255)
                                                                            ccAGCTTCTGGGGATGGTGGGCGTGCAGCGGGGTAATAATGGTTACAATCCCCCCCCAGTGAGGGGGGGCCTA‑AAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  2:408732/150‑1 (MQ=255)
                                                                             cAGCGTGTGGGTAATGGGGGCGTGCAGCGGGGTGATAATGGCTACAAGGACCCCTCGGTGATGGTGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                         <  2:226245/150‑1 (MQ=255)
                                                                               gCGTCTGGGGAATGTGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGATGGGGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg                                                                      <  1:1322474/151‑1 (MQ=255)
                                                                                  tgtGGGGAGGGTGGGAGTGAGGGGGGGTAATATTGGCAAAATTAACACCCGAGGGATGGTGGGACTA‑AAAAATAGTTATTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAGTATTATTTTTTg                                                                      <  1:1213144/146‑1 (MQ=255)
                                                                                     ggggAATGGGGGGGTGCGGCGGGTTAATAATGTCTACAATGCCCCCGTAGGGATGG‑GGGCCTGAAAATATAAGTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                                <  2:837794/151‑1 (MQ=255)
                                                                                              ggggGTGCAGCGGGGTTATAATGCCTAAAATGACCCCGGAGTTGTGGTGGCTATG‑AAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc                                                            <  2:197919/146‑1 (MQ=255)
                                                                                              gggcgggcAGGGGGGTAAAAATGGCAACAAAGCCCCCGAGGTGAGGGGGGGCAAA‑AATTATAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg                                                       <  2:111538/151‑1 (MQ=255)
                                                                                                  gTGGAGCGGGGTGATAATGGTAAAAATGACCCCGGGGTGTGGGTGGGCGTA‑AAAAATAGTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt                                                     <  2:517840/149‑1 (MQ=255)
                                                                                                      ggCTGGGTGATAAGGGCTAAAAGGCCCCCGCGGTGAGGGGGGGCATA‑AAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca                                                 <  1:74943/148‑1 (MQ=255)
                                                                                                               aaaaagccTACATTGCCCCCGCAGTGTGGGGGGGCTTA‑AAAAATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                                      <  1:1018052/144‑1 (MQ=255)
                                                                                                                tAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  1:272596/1‑151 (MQ=255)
                                                                                                                 aaTGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  1:747197/1‑150 (MQ=255)
                                                                                                                    ggCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                                  >  1:1032391/1‑150 (MQ=255)
                                                                                                                        aCAATTACCCCCAAGTTATGGTGGGCCAA‑AAATATAATTTTTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                                 <  1:596674/147‑1 (MQ=255)
                                                                                                                         cAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  1:1338223/1‑151 (MQ=255)
                                                                                                                                        gagggggggCTTA‑AAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  <  1:451503/140‑1 (MQ=255)
                                                                                                                                           ggtggtCCTG‑AAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               <  2:1111522/146‑1 (MQ=255)
                                                                                                                                           ggtggGCATA‑AAAAATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg            <  1:302323/149‑1 (MQ=255)
                                                                                                                                           ggggggCTTA‑AAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc           <  1:119177/147‑1 (MQ=255)
                                                                                                                                                   tG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:1334169/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  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: ‑