Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC45.6% 111.4 / 17.5 69intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (35/2);  new base C (0/31);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 2.41e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTATAGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065
                                                                                                                                                             |                                                                                                                                                   
aGCGGTATAGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGaaaaaaa                                                                                                                                                     >  2:430714/1‑150 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAAcccc                                                                                                                                                 <  2:608152/151‑1 (MQ=255)
    gTATAGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCccg                                                                                                                                                >  2:383912/1‑151 (MQ=255)
    gTATAGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCccg                                                                                                                                                >  1:183699/1‑151 (MQ=255)
      ataGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCccgcc                                                                                                                                              >  2:147601/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCgg                                                                                                                                            >  2:339703/1‑151 (MQ=255)
            tGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGc                                                                                                                                        >  1:564245/1‑151 (MQ=255)
             gAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑GTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                       >  1:112751/1‑151 (MQ=255)
                   aCGATAAACTG‑CCTTGTGTTATATGAGTTATCGGCATGATT‑TCCGTTTTCGGA‑GGGTTTTTT‑AATCGGACGTTTAGT‑GGGTAAAAGGGGGGT‑TTTTTT‑GGGGTTTACAGAAAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                  <  2:1219316/150‑1 (MQ=255)
                   aCGATAAACAG‑CCCTGCTTTATATGATTTAGCGGCATGATTACCCGTTTCTGCA‑GGT‑TTTTT‑AATCGGACGATTAGT‑GGGTAAAATGAGGGT‑TTTTTT‑GGGGGTTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCggg                                                                                                                                     <  2:777517/147‑1 (MQ=255)
                   aCGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGCTTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                    >  1:516840/1‑148 (MQ=255)
                   aCGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  2:1102063/1‑149 (MQ=255)
                        aaaCAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                            >  1:577205/1‑151 (MQ=255)
                         aaCAGCCCCTGTTTTATATGAGTTATCGGCATGATTAT‑CGTTTCCGGAGGGG‑TTTTT‑AACGGGACAATTGGG‑GGGGAAAAGGGGGGT‑TTTTTT‑GGGGGTAACGGGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                  <  2:890320/145‑1 (MQ=255)
                                  tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                    >  2:1073779/1‑150 (MQ=255)
                                   gctttATATGAGTTATCGGCATGATTATCCGTTTCTGCG‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTAAAATGAGGGG‑TTATTT‑GGGGTTTACAGGTAAATTCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                    <  1:38640/146‑1 (MQ=255)
                                      ttATATGATTTATCGGCATGATTATCCTTTTCTGCG‑GGG‑TTTTT‑AATCGGACGATTAGG‑GGGTGAAATGGGGGG‑TTTTTT‑GGGGGTTACAGTTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                                  <  2:1126459/148‑1 (MQ=255)
                                       tataTGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                              >  2:770319/1‑151 (MQ=255)
                                        agaagattaaggggCAGGTTTACCCGTTTC‑GAA‑GGT‑TTTTTAAATGGGAAAATTGGT‑GGGTAAAGGGAGGGT‑TTTTTG‑GGGGTTTACAGGGAATTCCC‑AGAGAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                               <  1:735033/137‑1 (MQ=255)
                                          atGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  2:1065679/1‑151 (MQ=255)
                                          atGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  1:1182562/1‑151 (MQ=255)
                                               ttacgggCATGATTATCCTTTTCGGAG‑GGT‑TTTTT‑AATCGGAATATTGGG‑GGGTGAAAGGAGGTT‑TTTTTT‑GGGGGTTAAGGTTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                            <  2:798178/140‑1 (MQ=255)
                                               ttATCGGCATGATTATCCGTTTCTG‑A‑GGG‑TTTTTAAATGGGACGATTAGG‑GGGTAAAATGAGGGG‑TTATTG‑GGGGGTTACAGGAAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           <  1:1104996/146‑1 (MQ=255)
                                                tATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                      >  2:166917/1‑150 (MQ=255)
                                                 aTCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAAAAATTGCGCTTTATg                                                                                                      >  2:299739/1‑149 (MQ=255)
                                                 aTCGGCATGATTATCCGTTTCTGCA‑GGG‑GTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                    >  2:344780/1‑151 (MQ=255)
                                                    ggaatatttACCCTTTTCTGAA‑GGG‑TTTTT‑ATGCGGACGATTAGG‑GGGGAAAAGGAGGGT‑TTTTTT‑GGGGGTTACAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                 <  2:174795/144‑1 (MQ=255)
                                                     gttgattttcccttttCTGGT‑GGT‑TTTTT‑AATCGGACAATAAGG‑GGGGAAAAGGGGGGT‑TTATTG‑GGGGGTTAAGGGTAAATACA‑GGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                     <  1:173030/133‑1 (MQ=255)
                                                     gaaTGATTGGCGGTTTCTATA‑GGG‑TTTTT‑AAGCGGAGGTTTAGG‑AGGTGAAAGGGGGCG‑TTTTTT‑GGGGGTTACGGATATTGTCA‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  <  2:105181/147‑1 (MQ=255)
                                                      cATGATTATCCTTTTCTGCG‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAAGGAGGGG‑TTTTTG‑GGGGGTTACAGGAAAATTCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCAgg                                                                                               <  1:383685/151‑1 (MQ=255)
                                                      cATGATTACCCGTTTCTGCA‑GGT‑TTTTT‑AACCGGACGATTAGG‑GGGAAAAATGAGGGT‑TTTTTT‑GGGGGTTACAGGTAAATCCC‑AGGCAAAAAACCCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                <  2:658204/150‑1 (MQ=255)
                                                       aTGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                               >  2:368715/1‑150 (MQ=255)
                                                       aTGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                >  1:190565/1‑149 (MQ=255)
                                                        tgtttATCCGTTTCTGCA‑GGT‑TTTTA‑AATCGGACGATTAGT‑GGGTAAAAGGGGGGT‑TTTTTT‑GGGGGTTACAGGTAATTTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                               <  2:367737/146‑1 (MQ=255)
                                                        tGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                             >  1:61804/1‑151 (MQ=255)
                                                           ttATCCGTTTCTGCA‑GGT‑TTTTT‑AACCGGACGATTAGG‑GGGTAAAATGAGGGG‑TTTTTG‑GGGGTTTACAGGAAAATTCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                              <  1:958523/147‑1 (MQ=255)
                                                               ccGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCACGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                       >  2:572584/1‑150 (MQ=255)
                                                                       gga‑GGT‑TTTTA‑AATCGGACGATTGGG‑GGGTAAAAGGGGGGT‑TTATTT‑GGGGGTAACAGGAAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                              <  2:722752/149‑1 (MQ=255)
                                                                            ggtttttt‑AATCGGACGATTAGT‑GGGTAAAATGAGGGT‑TTATTT‑GGGGTTTACAGGAAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                           <  1:870761/148‑1 (MQ=255)
                                                                            gg‑TTTTT‑AATCGGACGATTAGG‑GGGTAAAATGGGGGG‑TTATTG‑GGGGTTTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                          <  1:383185/151‑1 (MQ=255)
                                                                               ttttt‑AATCGGACGATTGGG‑GGGTAAAAGGGGGGT‑TTTTTT‑GGGGGTTACAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                         <  2:846144/150‑1 (MQ=255)
                                                                                tttt‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                       >  1:244677/1‑151 (MQ=255)
                                                                                  tt‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                     >  2:573764/1‑151 (MQ=255)
                                                                                   t‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    >  2:424598/1‑151 (MQ=255)
                                                                                   t‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    >  2:251543/1‑151 (MQ=255)
                                                                                     aaTCGGAGGATTGGG‑GGGTGAAATGAGGGT‑TTATTG‑GGGGTTTACAGTAAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    <  1:560324/150‑1 (MQ=37)
                                                                                     aaTCGGACGATTGGT‑GGGTAAAATGGGGGT‑TTATTG‑GGGGTTTACAGGTAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   <  1:941184/151‑1 (MQ=38)
                                                                                      aacgggagATTAGTGGGGGGAAAGGGGGGT‑TTTTTT‑GGGGGTGCCAGG‑AAATTCCGGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    <  2:747268/143‑1 (MQ=21)
                                                                                      aTCGGAGGATTGGG‑GGGTGAAAGGGGGGGTTTTTTT‑GGGGGTTACAGGTAAATCCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   <  2:1058763/151‑1 (MQ=37)
                                                                                          gACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                               >  2:1084181/1‑150 (MQ=255)
                                                                                            cGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                            >  2:121988/1‑151 (MQ=255)
                                                                                              aTTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                           >  1:956965/1‑150 (MQ=255)
                                                                                               ttAGT‑GGGTGAAATGGGGGT‑TTTTTT‑GGGGTTAACAGTTAATTCCC‑AGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                               <  2:254514/145‑1 (MQ=37)
                                                                                                tAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                        >  2:325276/1‑151 (MQ=255)
                                                                                                            aggggggg‑TTTTTTGGGGGGTTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                 <  1:430724/144‑1 (MQ=37)
                                                                                                            aTGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              >  1:353861/1‑150 (MQ=255)
                                                                                                              gggggt‑TTTTTT‑GGGGTTTACAGGTAATTTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                             <  2:334807/147‑1 (MQ=37)
                                                                                                              gAGGGT‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                            <  1:1184234/150‑1 (MQ=255)
                                                                                                                  ggttttttt‑GGGGGTTACAGGAAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  2:702622/143‑1 (MQ=37)
                                                                                                                             gggtTTAGAGGAAAATTCC‑AGGAAAAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                              <  1:7199/151‑1 (MQ=37)
                                                                                                                              gggTTACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                             <  2:619998/151‑1 (MQ=255)
                                                                                                                                   aCAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                        >  1:603422/1‑151 (MQ=255)
                                                                                                                                    cAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                       >  1:729812/1‑151 (MQ=255)
                                                                                                                                         aaaTTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   >  1:953356/1‑150 (MQ=255)
                                                                                                                                           aTTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  1:985567/1‑149 (MQ=255)
                                                                                                                                            ttCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa                 >  2:858142/1‑149 (MQ=255)
                                                                                                                                              cc‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             <  2:216203/151‑1 (MQ=255)
                                                                                                                                                  ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc               <  2:289898/141‑1 (MQ=255)
                                                                                                                                                  ggCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc          <  1:1102790/151‑1 (MQ=255)
                                                                                                                                                     aGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAgcgcagcgca       <  1:216646/151‑1 (MQ=255)
                                                                                                                                                      gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt      >  1:252778/1‑151 (MQ=255)
                                                                                                                                                      aaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt      <  2:1062826/150‑1 (MQ=255)
                                                                                                                                                       aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc     >  1:610188/1‑151 (MQ=255)
                                                                                                                                                        aaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa    <  1:526641/151‑1 (MQ=255)
                                                                                                                                                          aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  1:870148/1‑151 (MQ=255)
                                                                                                                                                             |                                                                                                                                                   
AGCGGTATAGAGTGAATTCACGATAAACAG‑CCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCA‑GGG‑TTTTT‑AATCGGACGATTAGT‑GGGTGAAATGAGGGG‑TTATTT‑GGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑