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

TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCG  >  NC_000913/1222086‑1222374
                                                                                                                                                        |                                                                                                                                              
tgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaaa                                                                                                                                             <  1:687397/151‑1 (MQ=255)
      tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAg                                                                                                                                      <  2:2411/151‑1 (MQ=255)
        gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                    >  1:1717245/1‑151 (MQ=255)
            aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                <  1:212329/151‑1 (MQ=255)
               aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATc                                                                                                                              <  1:1737305/150‑1 (MQ=255)
                 ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGCAACCGCGGTGTGGGGGGGCCTGAAATATAAGTTCTAATccc                                                                                                                            >  1:1068324/1‑150 (MQ=255)
                           aTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTa                                                                                                                 >  1:759294/1‑151 (MQ=255)
                                  gaaaaTTAACCTGGGGGAAAAGGGCACCCGAAACCCAGCGGCCAGCGTCGGGGGAATGGGGGGGTGAGGCTGGGTGATAATGGCTAAAATGCCCCCGCAGGGAGGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAAc                                                                                                           <  2:1419890/148‑1 (MQ=255)
                                  gaaaaTAAAATGTGGCGTAAAGGGCAACAGGAACCCAGGGCCAGGATTTTGGGGAATGGGGGCGTGCAGAGGGATGATATTGGCTACAAAGACCCCGCTGGGATGGGGGTGTTTAAAAAAAACTTCGCATCCCGCCGTTAGCTAAAAAAc                                                                                                           <  1:1098531/148‑1 (MQ=255)
                                  gaaaaTAAAACTTGGCGAAAGGGCCACCCTAAACCAGGCGGCCAGCGTCTGGGAAATGGGGGCGGGCGGCGGGGTAATAATGCCTAAAAGGACCCCGCAGTAATGGGGGGCTTAAAATAAAAGTTCTAACCCCGCCGTTAGCTAAAAAAc                                                                                                           <  1:1143428/148‑1 (MQ=255)
                                  gaaaaTAAAACTGGGCTAAAAGGCCACCCTGAGTCGGGCGGCCAGGGTGTGGGTAATGTGGGCGTGCAGGGGGGTAATATCGTATAAAATGCCCCCGTGTTGATGGGGGGCTGAAATTATAAGTTTTGATCCCGCCGTTAGCTAAAAAAc                                                                                                           <  1:571844/148‑1 (MQ=255)
                                    aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcg                                                                                                         >  2:1081089/1‑150 (MQ=255)
                                    aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGATACAATGACACCGCAATGATGGTGGGCCTGAAATATAAGTTCTAATCCCGCCGTTAGCTAAAAAACcg                                                                                                         >  2:597525/1‑150 (MQ=255)
                                    aaaTTAACTTGGGGAAAAATGAAACCAGAAACCAGGCGGCCGGCGTTGGGAGATTTAGGGGGTGGAGCTGGGTAATAATGAAAGCAATACCCCCGAAGGGAGGGAGGGCTTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACcg                                                                                                         <  2:1380055/150‑1 (MQ=255)
                                      tttAAGCTTGGGAAAAATGACACCCGGATCCCAGCGGCCAGCGTTTGGGAAGGGTGGGGGTGGAGGTGGGGAAAAATGTCTACAATGACACCGCAGTGGTGGGGGGCATAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgt                                                                                                      <  2:1406308/150‑1 (MQ=255)
                                             ttGGCGAAAATGGCACCCTAAACCCGGCGGCCGGCGTCGGGGGAATGGGGGCGTGCAGCTGGGTGAAAATGGCTAAAATGACCCCGAAGTGAGGGTGGGCCTAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                               <  2:98837/151‑1 (MQ=255)
                                              gggCGAAAAGGGAAACCGGAACCGAGCGGCCAGCGTGGGGGGAAGGGGGGGGGGCGGGGGGGTAATAGGGGCAAAAAGGACCCCGCAGGGAGGGGGGGCCTAAAAAAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                               <  2:1074924/149‑1 (MQ=255)
                                               ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                             >  2:1548560/1‑151 (MQ=255)
                                                gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGCGGGCCTGAAATATAAGTTCTGACCCCGCCGTTAGCTTAAAAACCGCGTCGTcatcatg                                                                                            >  1:1621786/1‑150 (MQ=255)
                                                 ctaaaaTGGCACCCTGAACCCGGCGCCCGGGGTCTGGGGATTGTGGGCGTGGAGCGGGGTGATATTGCCTACAATGCCCCCGCGTTGAGGGTGGGCCTAAAATAAAATTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                            <  2:453725/147‑1 (MQ=255)
                                                  gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGa                                                                                          >  2:699142/1‑151 (MQ=255)
                                                   aaaaTGGAACCCTAAATCCGGGGCCCGGGCTCGGGGGAGTGGGGGGGTGCGGCTGGTATAAAAGAGCTCAAATGCCCCGGCAGTGAGGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                         <  2:259775/150‑1 (MQ=255)
                                                     aaTGGAAACCAGAAACCGGGGCCCAGCGTTGGGGGAATGGGGGGGGGCAGTTGGGTAAAAATGGCTAAAAAGCCCCCGGAGTGTGGGGGGCCCTAAAAAATAATTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                       <  1:921435/151‑1 (MQ=255)
                                                     aaTGACAACCTGAACCCAGCGGCCAGCGTCGGGGAAATGGGGGCGGGCGGGTGGGGAATAAGGGCTACAAGACCCCCGAAGTGAGGGTGGCCTTAAAATAAAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                       <  1:1283082/151‑1 (MQ=255)
                                                       tGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                     >  1:956389/1‑151 (MQ=255)
                                                          cacccgaaacccAGGGGCCAGCGTCGGGGGAGTGGGGGCGTGCCGCTGGGTGATAAGGGCTACAATGACCCCGCGGTGAGGGGGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                  <  1:240193/141‑1 (MQ=255)
                                                             ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                               >  1:804087/1‑151 (MQ=255)
                                                               taaacccAGGGGCAAGGGTCGGGGGAAGGTGGGCGTGCTGCTGGGTAAAAAGGGCACAAAAGACCCCGCAGGGAGGGGGGCCCAAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  2:1692108/146‑1 (MQ=255)
                                                               tGAAGCCGGCGGCCAGCGGGGGGGGAATGGGGGGGGGCTGCGGGGTGATATTGCCTAAAATTACCCGGCAGTGAGGGGGGGCCTAAAAAATAAGTTGAGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  2:589313/151‑1 (MQ=255)
                                                               tGAACCCGGCGGCCAGGGTCGGGGGAAGGTGGGGTGGGGGGGGGGTGATATGGGCAACAAGGACGCCGCGGGGATGGGGGGCAGGAAAAATAGTTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  2:1320811/151‑1 (MQ=255)
                                                                tgagcccggcgcCCGGCGTCGGGGGAATGTGGGCGTGCAGGTGGGTGATAATGGCAACAAGGACCCCGCAGTGATGGTGGGCTTAAAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  1:818421/143‑1 (MQ=255)
                                                                taaacccggCGGCCGGCGTGGGGGGAGGGGGGGCGGGAGGCGGGTTAAAAAGGGTTACAATGACCCCGCGGTGAGGGGGGCCCGGAAAAAAAAGTCCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  2:810140/143‑1 (MQ=255)
                                                                taaaccagggggCCAGCGTCGGGGAATGGGGGGGTGGAGGTGGGGTGATAATGCCAAAAATGCCACCGAAGGGGGGGGGGGCATAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  1:923098/141‑1 (MQ=255)
                                                                gAACCCAGCGCCCGGCTTCGGGGGATTGGGGCGGGCCAGGGGGGTGATATGGGCTACAATCACCCCGCAGGGGTGGGGGGCTGAAAATAAAAGTCTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  2:1155054/150‑1 (MQ=255)
                                                                aaacccAGCGGCCGGCGTGTGGGGAAGGGGGGCGTGCGGGTGGGTGAAAATGGCAACAATGCCCCCGCAGTGTTGGGGGGCCTAAAATATAATTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  2:770600/146‑1 (MQ=255)
                                                                 gacgccggCGGCAGGCGTTGGGGGAAGGGGGGCGGGGAGGTGGGTGATAAGGGCAACAAGGACCCGGTGGTGAGGGTGGGCCAAAAATAAAAGTCCGGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  1:489457/143‑1 (MQ=255)
                                                                 gaccccggcgcCCAGCGTGGGGGGAATGGGGGGGTGGAGTGGGGAAAAAAGGGCTACAATGCCCCCGGAGTGATGGGGGCCATAAAATAAAGTTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  1:46154/143‑1 (MQ=255)
                                                                 aacccAGCGGCCAGCGTCGGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCGGTGATGGGGGGCTTAAAATATATGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                             <  1:1481849/146‑1 (MQ=255)
                                                                  aacccggCGGCAAGCTTCGGGGAAGTGTGGGCGTGGAGGTGGTAAATAATGCCTAAATTGCCCCCGCAGTGAGGGTGGGCCTAAAATATAAGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                           <  1:937993/145‑1 (MQ=255)
                                                                            taaaaagagaaaggATTGGGGGGGTGCAGGGGGTTAAAAAGGCCTAAAATGACCCCGAGGTGAGGGGGGGCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                  <  2:492472/139‑1 (MQ=255)
                                                                            ccAGCGTGGGGGGAATGGGGGGGGGCAGGGGGGTGATAATGGCAAAAATGACCCCGGAGGGAGGGGGGGCCTAAAAAATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                 <  2:659442/151‑1 (MQ=255)
                                                                                   tCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                           >  2:1552422/1‑151 (MQ=255)
                                                                                    cTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGa                                                          >  2:779620/1‑151 (MQ=255)
                                                                                         gAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                  >  1:492285/1‑138 (MQ=255)
                                                                                         gAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                     >  2:32879/1‑151 (MQ=255)
                                                                                         gAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                     >  2:280864/1‑151 (MQ=255)
                                                                                         gAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                     >  2:1650295/1‑151 (MQ=255)
                                                                                                   gggCAGCTGGGTAATAATGGCTACAATGCCCCCGAAGTGATGGTGGGCATAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                            <  1:1272118/148‑1 (MQ=255)
                                                                                                   gTGAAGTTGGGTGAAAATGGTAAAAATGCCCCCGCAGTGTTGGGGGGCCTAAAATAAAATTTTTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                            <  2:1472065/150‑1 (MQ=255)
                                                                                                     gtgcggtggggTGATAATGGCTAAAATGACACCGAAGGTAGGGGGGGCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                            <  2:582352/142‑1 (MQ=255)
                                                                                                          gtGGGTAAAAATGGCTACAAGGACCCGGAAGGGAGGGGGGCCAAAAAATATAAGTGTGAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                       <  1:440788/148‑1 (MQ=255)
                                                                                                          cTGGTTGATAAGGGCAACAAGGACCCCGAGGTGAGGGGGGGCTTAAAATAAATGTTCTGGTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                       <  1:1347788/149‑1 (MQ=255)
                                                                                                             gggTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt                                   >  1:1044814/1‑151 (MQ=255)
                                                                                                             gggTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt                                   >  1:1541061/1‑151 (MQ=255)
                                                                                                               gTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGGGTCATCAATTg                                  >  2:1593934/1‑150 (MQ=255)
                                                                                                               gTGATAATGGCTAAAATGCCCCCGCAGGGAGGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                  <  2:467176/150‑1 (MQ=255)
                                                                                                                ggaaaaggccTAAAATGACCCCGCAGTGAGGGGGGGCCTAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                 <  2:1330006/141‑1 (MQ=255)
                                                                                                                      tGGCAACAATGACACCGCAGTGAGGGGGGGCTTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                          <  2:957741/151‑1 (MQ=255)
                                                                                                                       ggCCACAATGACACTGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGGATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGCGTCATCAATTGCCGAGAGCa                         >  1:407122/1‑151 (MQ=255)
                                                                                                                        gCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCaa                        >  1:386925/1‑151 (MQ=255)
                                                                                                                         caaaaaTAACACCGAAGTGAGGGGGGGCATAAAATATATTTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                       <  1:1173174/147‑1 (MQ=255)
                                                                                                                         cTAAAATGACACCGCAGTGGGGGGGGCCCAAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                       <  2:1308692/151‑1 (MQ=255)
                                                                                                                          tACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                      >  2:1438975/1‑151 (MQ=255)
                                                                                                                            cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                     >  2:336696/1‑150 (MQ=255)
                                                                                                                                 acacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg               >  2:538278/1‑151 (MQ=255)
                                                                                                                                    ccGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg            >  1:1086142/1‑151 (MQ=255)
                                                                                                                                    ccGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg            >  2:935595/1‑151 (MQ=255)
                                                                                                                                    ccGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg            >  1:483190/1‑151 (MQ=255)
                                                                                                                                    ccGCAGTGATGGCGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg            >  1:1424000/1‑151 (MQ=255)
                                                                                                                                     cGCAGGGGTGGGGGGCGGAAAAGATAATTTCGGGCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg           <  2:1723511/151‑1 (MQ=255)
                                                                                                                                       cAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca                                      >  1:345332/1‑122 (MQ=255)
                                                                                                                                       cAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt          >  2:350628/1‑150 (MQ=255)
                                                                                                                                         gTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt       >  1:974783/1‑151 (MQ=255)
                                                                                                                                            aTGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg    >  2:981244/1‑151 (MQ=255)
                                                                                                                                                tggGCCTAGAAATAAATTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcg  <  1:568983/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                              
TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCG  >  NC_000913/1222086‑1222374

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