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

GAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222087‑1222382
                                                                                                                                                  |                                                                                                                                                      
gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    <  1:202505/151‑1 (MQ=255)
      cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                             <  1:247794/151‑1 (MQ=255)
       gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                            <  1:916586/151‑1 (MQ=255)
                ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATccc                                                                                                                                    <  1:386809/150‑1 (MQ=255)
                ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGGGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTCGGCCTGAAATATAAGTTCTGATCccg                                                                                                                                   <  2:914390/151‑1 (MQ=255)
                           tttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTa                                                                                                                         >  2:1176432/1‑150 (MQ=255)
                               ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaaa                                                                                                                     >  2:1293285/1‑150 (MQ=255)
                                   aaattaaattGGGTGAAAGGGGCACCCGGAACCCGGCGGCGGGCGTGGGGGGAATGTGGGGGTGGGGTTGGTTGATAATGGCTACAATGCCCCCGCGGTGATGGTGGGCATGAAAAATAAGTTCTGACCCCGCCGTTAGCTAAAAAAcc                                                                                                                  <  2:117990/149‑1 (MQ=255)
                                    aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               >  2:990861/1‑151 (MQ=255)
                                     aTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgt                                                                                                              >  1:1097266/1‑151 (MQ=255)
                                     aTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgt                                                                                                              >  1:263721/1‑151 (MQ=255)
                                        aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                           >  1:730503/1‑151 (MQ=255)
                                          aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt                                                                                                         >  2:169001/1‑151 (MQ=255)
                                             tGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                       >  1:933852/1‑150 (MQ=255)
                                             gggCGAAAAGGGCAACCGTAAGCCAGCGGCGAGCGTCTGGGCAGTGGGGGGGTGAGTCTGGGTGAAAATGGCTACAATGACCCGGGAGTAAGGGTGGCCCAGAAATTTCTTTTCATACCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                          <  1:1210546/146‑1 (MQ=255)
                                              gggttAAAGGGCACCCTGAACCCGGCGGCCGGCGTGGGGGGAATGTGGGCGTGCGGCGGGGTGATAAGGCTTAAAATGACCCCGCAGTGAGGGTGGGCCTAAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                                     <  2:1081689/147‑1 (MQ=255)
                                                  tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                 >  1:1081224/1‑151 (MQ=255)
                                                       ggaacccTGAACCCAGCGGCCAGCGTCGGGGGAATGGGGGGGTGCAGCTGGGTGAAAATGGCTACAATGCCACCGCAGTGATGGGGGGCCTGAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg                                                                                              <  1:1442639/144‑1 (MQ=255)
                                                       ggCACCCTGACCCCGGCGGCCGGCGTCGGGGGAGTGGGGGCGTGCAGCGGGGTGATAATGCCTACAATGACCCCGCGGTGATGGGGGGCTTAAAATATAAGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                             <  2:418856/150‑1 (MQ=255)
                                                          accccgAATCAAGCGACCAGGGTCGGGGGAATGTGGGCGTGCGGCTGGGTGAAAACGGCCACAAGGACCCCGCAGTGATGGGGGGCCAAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          <  1:261382/145‑1 (MQ=255)
                                                          aaCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  2:1440149/1‑150 (MQ=255)
                                                            ccTGAACCCGGCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGCTGGGTAATAAGGGCTAAAATGCCCCCGCAGTGATGGTGGCCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        <  1:1257157/150‑1 (MQ=255)
                                                                atcccggCGGCCGTCTTCGGGGGAATGTGGGCGTGCATCTGGGTGAAAATGGCTAAAATCACCCCGGAGGTTTGGTGGGCCTAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   <  1:1146730/145‑1 (MQ=255)
                                                                aacccAGCGGCCGGCGTCGGGGGAATGTGGGCGTACAGGTGGGTAATAATGGTTAAAATGCCACCGAGTTGATGGGGGGCCTAAAATATAAGTTCTACTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          <  2:1264715/141‑1 (MQ=255)
                                                                 acccAGCGGCCAGGGTTTGGGGAATGGGGGCGTGAGGGTGGGTTATAATAAATAAAATGACCCCGCAGTGTGGGGGGGCTTAAAATATAATTTCGAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa                                                                                  <  2:544514/149‑1 (MQ=255)
                                                                 acccAGCGGCCAGCGTCGGGGGAATGGGGGCGGGCAGGTGGGTGATAATGGCTACAATGACCCCGCAGTGAGGGTGGGCCTGAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa                                                                                  <  2:74984/149‑1 (MQ=255)
                                                                 aTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   >  1:444905/1‑150 (MQ=255)
                                                                   ccgcgggCCGGCGTCGGGGGATTGTGGGTGTCCCGTGGGTTAATAATGGCCAAAATCACCCCGCAGTGTTGGTGGTTTTAAAATATAAGTTCGGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  2:570938/146‑1 (MQ=255)
                                                                   ccAGCGGCCAGCTTTTGGGGAATGTGGGGGTGCAGGTGGTTAAAAATGACTAAAAGGCCACCGCGGGGAGGGGGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  2:1078002/151‑1 (MQ=255)
                                                                   ccAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  2:405386/1‑150 (MQ=255)
                                                                    cggggcccggggttgggggAATGGGGGGGTGCAGTGGGGTAAAAAGGGCTAAAATGACCCCGAGGTGAGGGTGGGCATAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  1:589064/135‑1 (MQ=255)
                                                                      ggggCCAGTGTTTGGAGATGGGGGGGGTGCAGCTGGGTGATAATGTCTATAATCACTCCGGAGTGGTGGGGGGCTTCAAATAAAGCTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                             <  2:177978/149‑1 (MQ=255)
                                                                      gcgcCCAGCTTCGGGGAAAGGTGGGCGTGCGGCGGGGAAAAAGTGCCTACAGGACCCCCGCAGTGAGGGGGGGCTTAAAATATATGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                               <  1:218169/149‑1 (MQ=255)
                                                                        ggCCGGCGTCGGGGGAATGGGGGCGTGCAGGTGGGTGATAATGGCTAAAATGCCCCCGCAGTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            <  1:416015/150‑1 (MQ=255)
                                                                        ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                           >  2:648056/1‑151 (MQ=255)
                                                                        ggCAAGGGTGTGGGGAATGGGGGGTTGCAGGTGGTTAATATTGGTTACAAGGACACCGAATGGATGGTGGGCTTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  2:572317/146‑1 (MQ=255)
                                                                          ccggCGTCGGGGGATTGTGGGGGTGCAGTTGGGTAATAATGGCTAAAATGACACCGCAGTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                         <  1:1318492/148‑1 (MQ=255)
                                                                           cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                        >  1:113541/1‑151 (MQ=255)
                                                                             gCGTGTGGGAAGTGTGGGCGTGCAGGTGGGTGAAATTGGCTAAAATACCCCCGAAGTGATGGGGGGCATAAAATATAAGTTCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg                                                                      <  2:209149/151‑1 (MQ=255)
                                                                               gTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    >  2:804522/1‑151 (MQ=255)
                                                                                   ggggAGTGTGGGCGTGCAGTTGGGTGATAATGGCTACAATGACCCCGTAGTGATGGTGGGCTTAAAATATAATTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                     <  2:469137/146‑1 (MQ=255)
                                                                                   ggggAGTGGGGGGTTGAAGCGGGTGGATATTGGCTACAACGACGCCCGGGTGAGGGGGGCCAAAAATAAAAAGTTCTTCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    <  1:719476/147‑1 (MQ=255)
                                                                                   ggggAATGGGGGGGGGCAGCGGGGTAATAATGGTTAAAATGACCCCGCGGTGAGGGTGGGCCGAAAATAAAAGTTCTTTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                                <  2:1158625/151‑1 (MQ=255)
                                                                                   ggggAATGGGGGCGTGCGGTTGGGTAATAATGGCAACAATGACCCCGCAGTGATGGTGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                 <  2:263692/150‑1 (MQ=255)
                                                                                     ggAATGTGGGCGTGCAGCTGGGTGATAATGGCTAAAATGACACCGCAGTGATGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc                                                              <  2:742118/151‑1 (MQ=255)
                                                                                       aaTGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGCGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc                                                              >  1:158375/1‑149 (MQ=255)
                                                                                          gtgGGCGGCCAGTTGGGTGATAATGGCTAAAATACCCCGGCAGTGGTGGGGGGCATAAAAAATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                           <  1:830465/149‑1 (MQ=255)
                                                                                          ggggggGTGCGGCGGGGTAATAATGGCTAAAAGGACCCCGAAGTGATGGGGGGCTTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                           <  1:1332929/147‑1 (MQ=255)
                                                                                          gggggCGTGCAGTTGGGTAAAAATGGCTACAATGCCCCCGCAGTGATGGTGGCCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCa                                                         <  1:624565/149‑1 (MQ=255)
                                                                                            gggCGTGCGGCGGGGTAATAATGGTAACAATGCCCCCGAAGTGATGGGGGGCTAAAATTAAAATTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                           <  2:698813/147‑1 (MQ=255)
                                                                                             ggCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                      >  1:976209/1‑151 (MQ=255)
                                                                                               agTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                    >  2:245426/2‑151 (MQ=255)
                                                                                                   catcggggTAATAAGGGCTAAAATGCCACCGCGTTGATGTTGGCCTTAAAATATAATTTCTTCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt                                                   <  2:1019613/143‑1 (MQ=255)
                                                                                                    aGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGccatc                                               >  2:602720/1‑151 (MQ=255)
                                                                                                             aaaagggCTACAATGACACCGAAGTGAGGGGGGGCCAAAAATAAATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                                      <  1:888327/146‑1 (MQ=255)
                                                                                                               aaTGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  1:574914/1‑150 (MQ=255)
                                                                                                                        aaTGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  2:992837/1‑150 (MQ=255)
                                                                                                                          tGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt                         >  1:855809/1‑151 (MQ=255)
                                                                                                                             ccccGCGGTGAGGGTGGGCATAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt                         <  2:556622/146‑1 (MQ=255)
                                                                                                                             cacCGCAGTGAGGGGGGGCCAAAAAAATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      <  1:609916/151‑1 (MQ=255)
                                                                                                                                cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                   >  2:298557/1‑151 (MQ=255)
                                                                                                                                  cAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc                 >  2:1247711/1‑151 (MQ=255)
                                                                                                                                    gTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               >  2:500696/1‑151 (MQ=255)
                                                                                                                                         ggtggCCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg            <  1:1270057/149‑1 (MQ=255)
                                                                                                                                           tggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:696060/1‑150 (MQ=255)
                                                                                                                                            gggCGTAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa        <  1:648130/150‑1 (MQ=255)
                                                                                                                                             ggCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa       >  2:235998/1‑150 (MQ=255)
                                                                                                                                             ggCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg            <  1:731121/145‑1 (MQ=255)
                                                                                                                                                cTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCt   >  1:380868/1‑151 (MQ=255)
                                                                                                                                                  gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:518810/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222087‑1222382

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