Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC40.7% 96.8 / 12.7 59intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/24);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001772‑2002065
                                                                                                                                                     |                                                                                                                                                   
ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAAccc                                                                                                                                                  <  1:139545/150‑1 (MQ=255)
 gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAACAACCccg                                                                                                                                                <  1:190690/151‑1 (MQ=255)
   ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCccgcc                                                                                                                                              >  2:134151/1‑151 (MQ=255)
   ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCccgc                                                                                                                                               <  2:839798/150‑1 (MQ=255)
     agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCgg                                                                                                                                            >  2:63139/1‑151 (MQ=255)
      gggTGAATTCCCGATAAACGGCCCTGGGTTATAGGGGTTATGGGGAGGATTATCGGTTGCTGCGGGGTTTTTAAGGGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAACCCccgccg                                                                                                                                             <  1:740440/147‑1 (MQ=255)
          gAATTCACGATAAACAGCCCTGCGTTATATGAGTTCTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCGAGG‑CAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                       >  1:12640/1‑151 (MQ=255)
            aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                     >  2:540358/1‑151 (MQ=255)
               cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  2:318224/1‑150 (MQ=255)
                aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  2:197907/1‑149 (MQ=255)
                 cGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAg                                                                                                                                 >  1:115010/1‑150 (MQ=255)
                    tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGCTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTATAGGTAAATTCCCGG‑CAGAAAAAAATCACGCCGGTGGCGGGGAAGCAc                                                                                                                              >  1:902690/1‑150 (MQ=255)
                            ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                     >  1:37944/1‑151 (MQ=255)
                                  ttttatgatttcggggCTTGTTTTTCCTTTTTTGCGGGGTTTTTAAGCGGAGGTTTAGGGGGTGAAA‑GGGGGGTTTTTTGGGGTT‑TTACAGTTAATTCCCAGG‑AAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                                  <  2:341018/134‑1 (MQ=255)
                                    atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTATTGGGGG‑G‑TTACAGGAAAATTCCAGG‑CAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                              <  2:594828/150‑1 (MQ=255)
                                      atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  1:388957/1‑151 (MQ=255)
                                       tGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                          >  2:163763/1‑151 (MQ=255)
                                        gAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGAGCttt                                                                                                         >  1:505627/1‑151 (MQ=255)
                                         aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                          >  2:834920/1‑149 (MQ=255)
                                           ttATCGGCTTGATTATCCGTTTCTGCGGGGTTTTTAATCGGACAATTAGTGGGTGAAAGGAGGGGTTATTTGGGG‑G‑TTACAGGTAATTTCCAGA‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                        <  1:681153/149‑1 (MQ=255)
                                              tCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGCAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCCCGTTGCTGACAAATTGCGCTTTATGtgg                                                                                                   >  1:954595/1‑149 (MQ=255)
                                                       ttATCTTTTTCTGCGGGTTTTTTAACAGGGCGTTAGGGGGGTGAAAAGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATCCCAGG‑CAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                             <  1:25115/148‑1 (MQ=255)
                                                        tttCCGTTTCGGGTGGTTTTTAAATCGGACGATTGGGGGTTGAAATGGGGGGTTATTTGGGG‑T‑TAAGAGGAAAATTCCAGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                          <  1:933727/148‑1 (MQ=255)
                                                         aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                         >  1:95368/1‑150 (MQ=255)
                                                              tttCTGCAGGTTTTTTAATCGAACGATTAGTGGGTAAAATGAGGGGTTTTTGGGGG‑T‑TTACGGGAAAATTCCAGG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                      <  2:215615/148‑1 (MQ=255)
                                                                   gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                              >  2:223229/1‑151 (MQ=255)
                                                                    cAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                             >  1:370516/1‑151 (MQ=255)
                                                                          ttttAATCGGACGATTAGGGGGTGAAATGAGGGTTTTTTTGGGG‑G‑TTACAGGTAAATCCCAGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                        <  1:846745/150‑1 (MQ=39)
                                                                              aaTCGGACGATTAGTGGGTGCAATGAGGGGCTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGCTGCGGCGTAAACGCCTTATCCGGCCt                                                                    >  1:655325/1‑150 (MQ=255)
                                                                              aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    >  1:81311/1‑150 (MQ=255)
                                                                                  ggACGATTGGTGGGTGAAATGAGGGGTTTTTTGGGG‑G‑TTACAGGTAAATCCCAGC‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                  <  1:318035/148‑1 (MQ=38)
                                                                                  ggACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGG‑G‑TTACAGTTAAATCCCAGG‑CAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                               <  2:530724/151‑1 (MQ=37)
                                                                                         tAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                        >  2:359565/1‑151 (MQ=255)
                                                                                          cttGGGTGAAATGAGTGGTTATTTCGGG‑G‑TTTCAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCTGGGAAGGACGTTGCTGACATAATGCGCTTTATGTTGCCTGATGCGGCGTAAACGCCTTATACGGCCTACaaaaa                                                             >  1:470133/3‑145 (MQ=38)
                                                                                            ggggTGAAAGGGGGGGTTATTTGGGG‑G‑TTACAGGTAATTTCCGGC‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                      <  1:720582/149‑1 (MQ=38)
                                                                                                tGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                  >  1:784222/1‑150 (MQ=255)
                                                                                                  aatgggggttttttttGGGG‑G‑TTACAGGTAAATCCCAGA‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               <  1:914781/138‑1 (MQ=37)
                                                                                                  aaaTGAGGGGTTTTTTGGGG‑G‑TTACAGG‑AAATTCCCGGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               <  2:313048/151‑1 (MQ=37)
                                                                                                  aaaTGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               >  1:506077/1‑151 (MQ=255)
                                                                                                  aaaTGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               >  1:856994/1‑151 (MQ=255)
                                                                                                   aaTGAGGGTTTTTTTGGGG‑G‑TTAAAGTTAAATTCCAGG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              <  1:970423/151‑1 (MQ=37)
                                                                                                      gAGGGTTTATTGGGGG‑T‑TTACAGGAAAATCCCAGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                           <  2:490453/151‑1 (MQ=37)
                                                                                                      gAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                           >  1:252850/1‑151 (MQ=255)
                                                                                                        ggggTTATTGGGGG‑G‑TTACAGGTAAATCCCAGG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                         <  2:324371/151‑1 (MQ=255)
                                                                                                        ggggTTATTGGGGG‑G‑TTACAGGTAAATCCCAGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                         <  2:727769/151‑1 (MQ=37)
                                                                                                         gggTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                        >  2:233526/1‑151 (MQ=255)
                                                                                                            ttATGTGGGG‑T‑ATACAGTTAACTTCCAGC‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                        <  1:473194/148‑1 (MQ=25)
                                                                                                             tttttGGGG‑GTTTACAGGTAAATTCCAGG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                      <  2:37952/148‑1 (MQ=255)
                                                                                                              aTTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    >  1:774352/1‑150 (MQ=255)
                                                                                                               ttggggg‑g‑TTACAGGTAAATCCCAGG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                                     <  2:243331/145‑1 (MQ=255)
                                                                                                                  gggg‑g‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                               >  2:137176/1‑151 (MQ=255)
                                                                                                                          tACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                         >  1:363345/1‑151 (MQ=255)
                                                                                                                                 aattcccggg‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                      <  2:972741/139‑1 (MQ=255)
                                                                                                                                 aaaTTCCAGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  <  1:420417/151‑1 (MQ=255)
                                                                                                                                  aaTTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  2:248969/1‑150 (MQ=255)
                                                                                                                                  aaTTCCAGG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                       <  1:771889/145‑1 (MQ=255)
                                                                                                                                         gg‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca            <  1:197795/149‑1 (MQ=255)
                                                                                                                                          g‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg         >  1:110942/1‑151 (MQ=255)
                                                                                                                                             aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc        >  1:647260/1‑150 (MQ=255)
                                                                                                                                               aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc     >  2:465150/1‑151 (MQ=255)
                                                                                                                                                aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa    <  2:126721/151‑1 (MQ=255)
                                                                                                                                                 aaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg   <  1:721833/151‑1 (MQ=255)
                                                                                                                                                  aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  <  1:63118/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGG‑G‑TTACAGGTAAATTCCAGG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001772‑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: ‑