Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA55.9% 28.5 / 87.5 68pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (24/6);  new base A (0/38);  total (24/44)
Fisher's exact test for biased strand distribution p-value = 3.95e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-09

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222083‑1222382
                                                                                                                                                         |                                                                                                                                                     
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑Tg                                                                                                                                                       >  2:1543/1‑151 (MQ=255)
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGAGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑Tg                                                                                                                                                       >  1:115913/1‑151 (MQ=255)
             aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCt                                                                                                                                          >  2:116159/1‑151 (MQ=255)
              cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTg                                                                                                                                         <  1:241644/151‑1 (MQ=255)
                 gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATc                                                                                                                                      >  2:689104/1‑151 (MQ=255)
                    ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATccc                                                                                                                                    >  2:618478/1‑150 (MQ=255)
                     gggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCccgc                                                                                                                                  <  2:125321/151‑1 (MQ=255)
                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                           >  2:457171/1‑151 (MQ=255)
                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                           >  1:734095/1‑151 (MQ=255)
                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                           >  2:906420/1‑151 (MQ=255)
                                ttGGGGGAAATTAACCTTGGCGTAAAGGGCACCCTGAACCCGGCGGCCGGCTTCGGGGAGAGGTGGGCGTGCAGGTGGGTGAAAAGGGCT‑ACAATGACACC‑GCAGTGGTGGTGGACT‑TAAAATATATTTTCTGATCCCGCCGTTAGCTaa                                                                                                                        <  1:445965/150‑1 (MQ=255)
                                 gggggaAAATAAAACTTGGCGAAAAGGGCACCCTAAACCCAGCGCCCGGCGTCGGGGGAATGTGGGCGTGCAGGGGGGTAAAAATGGCT‑ACAATGACCCC‑GAGGTGAGGGTGGGCT‑GAAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        <  1:429995/148‑1 (MQ=255)
                                            aacctgggCGAAAAGGGCACCCGGAACCCGGCGGCCAGGTTCGGGGGAATGGGGGCTTGCAGGTGGGTGATAATGGCT‑AAAATGCCCCC‑GGAGTGAGGGGGGGCT‑TAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                             <  2:762854/143‑1 (MQ=255)
                                            aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                           >  2:362658/1‑151 (MQ=255)
                                                   gcgaaaaggcaaaCCGGAACCCAGCGGCCGGGTTTGGGGGAAGGGGGGGGTGGAGGTGGGTGATAATGGCT‑AAAAGGCCCCG‑GAAGTGATGGGGGGCT‑TAAAATTTAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  1:847420/140‑1 (MQ=255)
                                                      aaaaTGGCACCCGAAATCCGGCGGCCAGCGTGGGGGAAATGGGGGGGTGCAGTTGGGTAATAATGCCT‑AAAATCCCCCC‑GAAGTGATGGGGGGCT‑AAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   <  1:119848/148‑1 (MQ=255)
                                                       aaaTGGCCACCTGAATCCAGCGGCCAGCGTCGGGGAATTGGGGGCGTGCGGCTGGGTGATAATGCCT‑AAAATGCCACC‑GCAGTGATGGTGGGCT‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   <  2:274311/148‑1 (MQ=255)
                                                       aaaGGGCAACCTGACTCCGGCGGCGAGCGTGTGGGGAGTGGGGGGGTGCAGCTGGGTGATATTGGCT‑ACAATGACCCC‑GCAGTGATGGTGGGCC‑TGAAATATATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGa                                                                                                  <  2:1073751/149‑1 (MQ=255)
                                                          tGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                             >  2:471027/1‑151 (MQ=255)
                                                           ggCACCCTGAATCCAGCGGCCGGCGTCTGGGGAAGGGGGGGGTGCAGCGGGGTGAAAATGCCT‑AAAATGCCACC‑GCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                            <  2:850835/151‑1 (MQ=255)
                                                              acccgaaacccAGCGGCCGGCGTGGGGGGAATGTGGGCGTGCAGCTGGGTAATAATGGCT‑AAAACACCCCC‑GAATTGATGGTGGGCC‑TGAAATATAATTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          <  2:991041/141‑1 (MQ=255)
                                                              acccTGAATCCAGCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCT‑ACAAGGCCCCC‑GCGGTGATGGGGGCCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                         <  2:990287/149‑1 (MQ=255)
                                                               aCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        >  1:830147/1‑151 (MQ=255)
                                                                ccgaaacccAGGGCCCGGCGTCGGGGGAATGGGGGGGGGCAGTTGGGTGATAATGGCT‑ACAATGCCCCC‑GAAGTGAGGGGGGGCT‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        <  1:230697/143‑1 (MQ=255)
                                                                ccTGAACCCGGCGGCCGGCGTCTGGGAAGTGGGGGCGTGCAGCTGGGTAATAAGGGCT‑AAAATGACCCC‑GCAGTGAGGGGGGCCC‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        <  2:549016/150‑1 (MQ=255)
                                                                ccTGAACCCAGCGGCCGGGGTCGGGGGATTGTGGGCGTGCAGCTGGGTAATAATGGCT‑ACAATGACCCC‑GAAGTGGTGGTGGGCC‑TAAAATATGAGTTCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                         <  1:1041837/149‑1 (MQ=255)
                                                                ccTAAACCCAGCGGCCAGCGTCTGGGAAATGTGGGCGGGCAGCTGGGTGATAAGGGCA‑AAAATGACCCC‑GAGGTGTTGGTGGCCC‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  2:549141/151‑1 (MQ=255)
                                                                 cTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       >  2:630031/1‑150 (MQ=255)
                                                                   gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                     >  2:908842/1‑150 (MQ=255)
                                                                   aacccccggggCCGGCGTCGGGGGAATGGGGGCGCGCAGTGGGTTAATATGGGCT‑AAAATAACACC‑GAAGTGATGGGGGCCA‑TGAAATATATTTCCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                                    <  2:38763/142‑1 (MQ=255)
                                                                    aacccAGCGGCCAGCGTCTGGGGATTGGGGGCGTGCAGCTGGGTGATAATGGCT‑AAAATGACCCC‑GCGGTGATGGGGGGCC‑TAAAATAAAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   <  2:1040827/148‑1 (MQ=255)
                                                                    aacccAGCGGCCAGCGTCTGGAGAGTGGGGGCGTGCAGCTGGGTGATAATGCCT‑AAAACGCCCCC‑GCATTGGGGGTGGGCT‑TAAAATATAGTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                      <  1:987047/145‑1 (MQ=255)
                                                                    aaTCCGGCGCCCAGCGTCGGGGGAATGTGGGCGTGCAGCGGGGTGATAAGGGCT‑ACAATGACCCC‑GCAGTGAGGGGGGGCC‑GAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                      <  2:1034986/148‑1 (MQ=255)
                                                                          gCGGCAGGCGTCGGGGGAATGGGGGGGTGCAGCTGGGTGATAATGGCT‑AAAATGACACC‑GCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              <  1:498794/150‑1 (MQ=255)
                                                                            ggCAAGCGTCGGGGGAATGTGGGGGTGAGGTTGGGTGATAAGGGCT‑ACAAGGACCCC‑GCAGGGATGGTGGGCT‑TAAAATATAAGTTCGGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            <  2:430090/150‑1 (MQ=255)
                                                                              caaGCGTCGGGGGAATGGGGGCGTGCAGTGGGGTGAAAAGGGCT‑AAAATGACCCC‑GAGGTGAGGGGGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  1:551607/148‑1 (MQ=255)
                                                                               ctgggtgtGGGGAGGGGGGGCGTGCAGCGGGGTTA‑AAAGGCTAAAAAAGACACC‑GAAGTAATGGGGGGCG‑TAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  1:67581/142‑1 (MQ=255)
                                                                               cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                        >  1:1029803/1‑151 (MQ=255)
                                                                               cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                        >  1:56728/1‑151 (MQ=255)
                                                                                aGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       >  2:162761/1‑151 (MQ=255)
                                                                                 gCGTGTGGGGAGTGTGGGGGTCGAGCGGGGTGAAAATGGCT‑ACAATAACCCC‑GCGGTGGGGGGGGGCA‑TAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  2:737382/147‑1 (MQ=255)
                                                                                 gCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       >  1:654113/1‑150 (MQ=255)
                                                                                       ggggAATGTGGGCGTGCAGTGGGGTGATAATGGCT‑AAAATGCCCCC‑GCGGTG‑TGGGGGGCCTTAAAATATAATTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                 <  2:1032002/150‑1 (MQ=255)
                                                                                         ggAAGGGGGGGGTGCAGCGGGGTGATAATGGCT‑AAAATGACCCC‑GCAGGGAGGGTGGGCC‑AAAAATAAAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc                                                              <  2:630712/151‑1 (MQ=255)
                                                                                              gtgGGGGTGCAGCGGGGTGATAAGGGCT‑AAAATGACACC‑GCAGTGAGGGTGGGCC‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc                                                            <  2:768816/148‑1 (MQ=255)
                                                                                              gtgGGGGGGCGGTTGGGTGATAATGCCT‑ACAAGGCCCCG‑GGGGGGGGGGGGGGCC‑TAAAATAAAGTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc                                                            <  1:642537/148‑1 (MQ=255)
                                                                                              ggggggGTGCAGCTGGGTAATAATGCCT‑AAAATGACCCC‑GCAGTGAGGGGGGGCC‑TAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCa                                                         <  2:624835/149‑1 (MQ=255)
                                                                                              gggggCGTGCAGCTGGGTAATAATGGCT‑ACAATGACCCC‑GCGGTGATGGGGGGCT‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                          <  2:371033/148‑1 (MQ=255)
                                                                                               tgGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                                        >  2:703330/1‑151 (MQ=255)
                                                                                                 ggttggCAGCTGGGTAATAA‑GGGT‑AAAATGGCCCCGGCAGTGAGGGTGGACC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                      <  1:314427/146‑1 (MQ=255)
                                                                                                 ggCTTGCGGCGGGGTAATAAGGGCT‑ACAATCACCCC‑GGTGTGGTGGTGGGCC‑TAAAATAAATGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                      <  2:885067/151‑1 (MQ=255)
                                                                                                      gcagtgggttGAAAATGGCA‑AAAATGACCCC‑GGAGTGATGGTGGGCC‑TAAAATATATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca                                                 <  2:159221/142‑1 (MQ=255)
                                                                                                       cAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcat                                                >  1:957004/1‑151 (MQ=255)
                                                                                                         ggtGGGTGATAATGCCT‑AAAATGACCCC‑GAAGTGATGGGGGGCC‑TAAAATATAATTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca                                              <  1:952133/149‑1 (MQ=255)
                                                                                                          cTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca                                              >  2:79459/1‑150 (MQ=255)
                                                                                                            ggttGATAATGGCA‑AAAATGACCCG‑GAGGTGATGGGGGGCA‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                            <  2:653489/147‑1 (MQ=255)
                                                                                                               tGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGcc                                        >  2:475257/1‑151 (MQ=255)
                                                                                                                   aaTGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  2:247720/1‑150 (MQ=255)
                                                                                                                   aaTGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  2:331558/1‑150 (MQ=255)
                                                                                                                      ggCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                                  >  1:358418/1‑150 (MQ=255)
                                                                                                                         t‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                               >  1:295754/1‑150 (MQ=255)
                                                                                                                           aaaaggccccc‑GAAGGGAGGGTGGGCT‑GAAAATAAAATTCCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                               <  1:408727/140‑1 (MQ=255)
                                                                                                                              aataccccggCAGTGAGGGTGGGCT‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGCGATCAACAGc                            <  1:116148/141‑1 (MQ=255)
                                                                                                                                   acC‑GCGGTGGTGGTGGGCT‑TAAAATATAATTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      <  1:689004/150‑1 (MQ=255)
                                                                                                                                    cc‑GCGGTGATGGGGGGCC‑TAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                       <  2:1094139/148‑1 (MQ=255)
                                                                                                                                      cgggttgtgggtggGCA‑AAAAATAAAAGTCCGAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg                    <  1:625802/141‑1 (MQ=255)
                                                                                                                                            gATGGGGGGCT‑TAAAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  <  1:618386/146‑1 (MQ=255)
                                                                                                                                               ggtggGCC‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc                 <  2:245582/144‑1 (MQ=255)
                                                                                                                                               ggtggGCC‑AGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcg          <  1:817451/151‑1 (MQ=255)
                                                                                                                                                        tGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:569284/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCT‑ACAATGACACC‑GCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  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: ‑