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

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222083‑1222382
                                                                                                                                                       |                                                                                                                                                      
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                                        >  2:532163/1‑151 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAata                                                                                                                                                   <  1:1209068/151‑1 (MQ=255)
                      ggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgcc                                                                                                                                 >  1:400525/1‑151 (MQ=255)
                        cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGt                                                                                                                               >  1:1213855/1‑151 (MQ=255)
                          aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  2:261303/1‑150 (MQ=255)
                              aTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCt                                                                                                                          >  2:764597/1‑150 (MQ=255)
                               tttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTa                                                                                                                         >  1:284478/1‑150 (MQ=255)
                               ttgggggaAAATTAAACTTGGCGAAAAGGGCACCCGAAACCCAGCGGCCAGCGTCGGGGGAATGTGGGGGTGCGGCTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTa                                                                                                                         <  2:400518/147‑1 (MQ=255)
                                        aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               >  2:425995/1‑151 (MQ=255)
                                        aaTTAAACTTGGCGTAAAGGGCACCCTAAACCCAGCGGCCAGCGTCTGGGGATTGGGGGGGGGCAGGTGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCATAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               <  1:975644/151‑1 (MQ=255)
                                           tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  1:314654/1‑151 (MQ=255)
                                             aaCTTGGCGTAAATGGCACCCTAAACCCGGCGGCCAGCGTCGGGGGAATGGGGGGGTGCAGGGGGGTAATAAGGGCTAAAATGACCCCGCAGTGATGGGGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                          <  1:580973/151‑1 (MQ=255)
                                               cTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        >  2:906021/1‑151 (MQ=255)
                                                ttGGCGAAAATGCAACCCTAAACCCAGCGGCCAGCGTCGGGGGAAGGGGGGCGGGCAGCTGGTTAAAAATGGCTAAAACGACCCGGAGGTGATGGTGGCCTTAAAATAAATTTTTTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  2:1220718/150‑1 (MQ=255)
                                                 gggggaaaaggccacccTGAACCCAGGGTCCGGCGTCGGGGGAATGTGGGCGGGCAGGTGGGTGATAATGCCTAAAATGACCCCGCAGTGTTGGTGGCCCTGAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                       <  1:147485/135‑1 (MQ=255)
                                                  ggctaaaaGGGCAACCGAAATCCGGCGCCCAGCGTCGGGGGAATGGGGGCGTGCAGTGGGGTAATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTGAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  1:1280288/143‑1 (MQ=255)
                                                  ggCGAAAAGGACACCCTGAACCCAGCGGCCAGCGTCGGGGAAATGGGGGCGGGCAGCGGGGTGATAATGTCTAAAAGGACCCCGCAGTGGGGGGGGCCCAAAAAAAAAGGTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  1:1143481/148‑1 (MQ=255)
                                                   gCGTAAATGGCACCCTGAACCCAGCGGCCAGCGTCGGGGAAATGGGGGGGTGCAGCTGGGTAATAATGGCTACAATGCCCCCGCAGTGAGGGGGGGCCTAAAAAATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  1:142822/151‑1 (MQ=255)
                                                   gCGTAAATGGCACCCTAAACCCAGCGGCCAGCGTCTGGGGAATGGGGGCGTGCAGGTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  1:754351/151‑1 (MQ=255)
                                                       aaagggaacccggAACCCAGCGGCCAGGGTCTGGGAATGGGGGGCGTGCAGTTGGGTAATATTGGTTACAAGGCCCCGGCGGTGATGGTGGGCTTAAAATATAGGTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   <  2:985173/136‑1 (MQ=255)
                                                       aaaTGGGCACCTGAAACCCAGCGCCCAGGGTGGGGGAAGGGTGGGCGGGCGGTTGGGTGATAATGCCTAAAATGACCCCGCGGTAGGGGTGGGCCTAAAATATAATTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  2:688321/148‑1 (MQ=255)
                                                       aaaTGGCACCCTGAACCCGGCGGCCAGCGTCTGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCTACAATGCCCCCGCGGTGATGGTGGCCCTGAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                <  1:760013/151‑1 (MQ=255)
                                                       aaaTGGCACCCTAAACCCGGCGCCCGGCGTTGGGGAAATGGGGGGGGGCGGTGGGGTAAAAATGGCTAAAATACCCCCGAGGTGATGGTGCCCAAAAAAATATTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                 <  1:44401/149‑1 (MQ=255)
                                                       aaaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                >  1:1018466/1‑151 (MQ=255)
                                                                ccTGAACCCAGCGGCCGGCGTGTGGGAAATGGGGGCGGGCAGCGGGGTTATAATGGCTAAAAGGCCCCCGCGGGGAGGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  1:484274/147‑1 (MQ=255)
                                                                 cgagattCCAGGGGCCAGCTTCTGGGGAATGGGGGCGTGCGGGTGGTTAAAAAGGGCTGAAAGCACACCGTAAGGATGGTGGCCTTAAAATAAATTTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                         <  1:1077984/143‑1 (MQ=255)
                                                                     aacccAGGGCCCAGCTTCGGGGGAAGGGGGGGGTCAAGCGGGTTAATAATGGCAACAATCACCCCGTTGTGATGGTGGGCTTAAAAAAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                                    <  2:321538/147‑1 (MQ=255)
                                                                     aacccAGCGGCCAGCGTGTGGGGAATGGGGGCGGGCAGCTGGGTGATAAGGGCTAAAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   <  1:1134688/148‑1 (MQ=255)
                                                                      acccggCGGCCAGCGTCGGGGAAATGGGGGCGTGCAGCTGGGTGATAATGCCTACAATGACACCGGAGTGTTGGGGGCCTTAAAAAAAATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        <  1:532171/140‑1 (MQ=255)
                                                                           gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              >  1:872636/1‑150 (MQ=255)
                                                                            cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            >  1:327074/1‑151 (MQ=255)
                                                                            cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            >  2:142747/1‑151 (MQ=255)
                                                                               ccAGCGTCGGGGGAATGTGGGCGTGCAGCTGGGTGATAAGGGCAACAATGACCCCGCGGTGAGGGGGGCCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  2:1037379/150‑1 (MQ=255)
                                                                                    gTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    >  1:1221079/1‑151 (MQ=255)
                                                                                     tggggggATGGTGGGGGTGCAGCTGGTTAATAATGTCTAAAATGCCCCCGCAGTGTGGGGGGGCTTGAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                     <  2:572877/146‑1 (MQ=255)
                                                                                     tcgggggATTGTGGGGGTGCGGCGGGGTAAAAATGCCTAAAATGACCCCGAAGTGAGGGGGGGAATGAAATAAAAGTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       <  1:350978/144‑1 (MQ=255)
                                                                                        ggggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                 >  2:795908/1‑150 (MQ=255)
                                                                                        ggggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                 >  2:458727/1‑150 (MQ=255)
                                                                                          ggAATGGGGGCGGGGAGCTGGGTAAAATTGACTAAATTGACCCCGAAGTAAGGGTGGGCATAAAAAATATTTTCTGATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       <  2:481502/142‑1 (MQ=255)
                                                                                                 ggggGTGCAGGTGGGGAAAAAGGGCTACAATGCCCCCGCAGTGATGGGGGGCCTGAAATATAAGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                                        <  1:31876/150‑1 (MQ=255)
                                                                                                     gTGCAGGGGGGTGATAATGGCTACAATGACCCGGCAGTGATGGTGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt                                                     <  2:282891/149‑1 (MQ=255)
                                                                                                              gggaaaaaggcataaaaTGCCCCCGGAGGGTTGGGGGGCTTAAAATAAATTTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                          <  2:969716/136‑1 (MQ=255)
                                                                                                               gtaaTAATGGCTACAATCACACCGGGGAGATGGTGGGCCAGAAAAATATTTGTTGATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa                                             <  2:86163/144‑1 (MQ=255)
                                                                                                                tGATAATGGCTACAATGACCCCGCAGTGTGGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         <  1:799575/150‑1 (MQ=255)
                                                                                                                  aaaaTGGCTAAAATGCCCCCGCAGTGATGGGGGGCTTAAAATATATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         <  2:947393/146‑1 (MQ=255)
                                                                                                                  aTAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg                                       <  1:10051/150‑1 (MQ=255)
                                                                                                                   tAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  2:479937/1‑151 (MQ=255)
                                                                                                                    aagggCTAAAATCACCCCGAGGTGATGGGGGGCTTAAAATATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     <  2:486911/147‑1 (MQ=255)
                                                                                                                           aCAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                               >  2:1274268/1‑149 (MQ=255)
                                                                                                                             aaTGACCCCGCAGTGAGGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                             <  1:1129662/149‑1 (MQ=255)
                                                                                                                              aTGACCCCGAAGTGAGGGGGGGCAAAAAATATAGGTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                          <  2:566264/151‑1 (MQ=255)
                                                                                                                               tGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt                         >  1:436688/1‑151 (MQ=255)
                                                                                                                                gACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        >  2:213490/1‑151 (MQ=255)
                                                                                                                                  cacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      >  2:1096703/1‑151 (MQ=255)
                                                                                                                                     cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                   >  2:205833/1‑151 (MQ=255)
                                                                                                                                      gcggtggtggtgGGTTTAAAATATAATTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        <  2:200866/142‑1 (MQ=255)
                                                                                                                                         ttGATGGGGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  <  2:314649/147‑1 (MQ=255)
                                                                                                                                             tggtggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc           >  1:1104963/1‑151 (MQ=255)
                                                                                                                                              ggtggCCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               <  1:1219494/146‑1 (MQ=255)
                                                                                                                                                    ccTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAcc    <  1:818065/151‑1 (MQ=255)
                                                                                                                                                       gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  1:1179423/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222083‑1222382

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