Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA40.3% 93.1 / 37.7 77intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (38/8);  new base A (0/31);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 2.31e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC  >  NC_000913/2001763‑2002054
                                                                                                                                                          |                                                                                                                                                
acCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGa                                                                                                                                                 >  2:218859/1‑150 (MQ=255)
acCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAGCGGCATGATTATCCGTTGCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGCGTTATTTGGGGGTTACAGGTAAATTCCATGCAGaa                                                                                                                                                <  1:772474/151‑1 (MQ=255)
     cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaa                                                                                                                                            <  1:491416/150‑1 (MQ=255)
       gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                         >  2:158859/1‑150 (MQ=255)
         ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                      >  2:1144799/1‑151 (MQ=255)
            ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                   >  2:1155231/1‑151 (MQ=255)
            ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                    <  1:1059090/150‑1 (MQ=255)
                    aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgg                                                                                                                           >  1:544615/1‑151 (MQ=255)
                     aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                          >  2:1083861/1‑151 (MQ=255)
                      ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                         >  1:1002465/1‑151 (MQ=255)
                        cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGGTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                        >  1:831640/1‑150 (MQ=255)
                           gaaaaaCGGCCCTGTTTTATAGGATTTATGGGAATGTTTATCCTTTTCTGCGGGGTTTTTAATCGGAGGTTAGGGGGGTGAAATGAGGGTTTTTTTGGGGTTTAAGGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                     <  1:490696/147‑1 (MQ=255)
                            aaaaaGAGCCCGGTGTTATATGATTAAGCGATATTTTTACCCTTTTCGGGGCGTTTTTTAACCGGACTTTTAGTGGGTGAAATGAGGGGTAGTGGGGGGGTGACAGCTAAATTCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                   <  2:42956/149‑1 (MQ=255)
                              aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                 >  1:159396/1‑151 (MQ=255)
                                    cccTGCTTTATAGGATTTATGGGAATATTTATCCTTTTCTGCGGGTTTTTAAATCGGACGATTGGTGGGTAAAAGGGGGGGTTTTTTGGGGGTAACAGAAAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                              <  1:79841/148‑1 (MQ=255)
                                             atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                  >  2:675943/1‑151 (MQ=255)
                                               atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                >  2:159639/1‑151 (MQ=255)
                                                tgttttACCGGCATGATTGTCCGTTTCGGCAGGTTTTTTAAGCGGACATTTGGTGGGTAAAAGGGGGGTTATTTGGGGGGTTCAGGGAAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                    <  2:287932/142‑1 (MQ=255)
                                                    ttATGGGCTTGATTTCCCTTTTCTGCGGGGTTTTTAACCGGACGTTGAGGGGGTAAAATGGGGGTTTATTGGGGGTTTACAGGTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                            <  2:695177/150‑1 (MQ=255)
                                                     tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                          >  2:754276/1‑151 (MQ=255)
                                                     tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAAGTTGCTGACAAATTGCGCTTTAt                                                                                            >  1:315964/1‑149 (MQ=255)
                                                         ggaaTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTAAAATGAGGGTTTATTGGGGGTTTAAAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGATTTATGTTg                                                                                        <  2:453514/146‑1 (MQ=255)
                                                         ggCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                      >  1:416816/1‑151 (MQ=255)
                                                          gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACACATTGCGCTTTATGTTGcc                                                                                      >  1:693116/1‑150 (MQ=255)
                                                           cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                    >  2:212574/1‑151 (MQ=255)
                                                                ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                               >  2:677616/1‑151 (MQ=255)
                                                                ttATCCGTTTCGGCAGGTTTTTTAATCGGAGGATTAGTGGGTGAAATGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                  <  2:254413/148‑1 (MQ=255)
                                                                  aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                             >  1:1185665/1‑151 (MQ=255)
                                                                  aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                             >  1:251858/1‑151 (MQ=255)
                                                                       tttCGGCTGGTTTTTTAATCGGACGATTAGGGGGTGAAAGGGGGGTTTTTTTGGGGTTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                             <  1:687240/146‑1 (MQ=255)
                                                                            gcgggttttttAATAGGACGATTGGTGGGAAAAAGGAGGGGTTTTTGGGGGTTTACAGATAATTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                    <  2:415052/144‑1 (MQ=255)
                                                                            gcggggTTTTTAATCGGACGATTAGGGGGTGAAAGGAGGGTTTATTGGGGGTTTACAGGTAAATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                   <  2:307954/148‑1 (MQ=255)
                                                                               ggttttttAATCGGAGATTTGGGGGGTGAAAGGAGGGTTTATTAGTGGGTAAAAGTTAAATTCAAACAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                <  2:11924/148‑1 (MQ=255)
                                                                               ggttttttAATCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTTGGGGGTTAAGGGAAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                       <  1:1143985/141‑1 (MQ=255)
                                                                                   ttttAATCGGACGATTAGTGGGTAAAAGGAGGGTTTATTTGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                             <  2:754946/150‑1 (MQ=39)
                                                                                    tttAATCGGACGTTTAGTGGGTAAAATGGGGGTTTATTTGGGGTTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                            <  2:320852/150‑1 (MQ=38)
                                                                                    ttaaaTGGGATGATTAGGGGGAGAAATGGGGGTTTTTTTGGGGGTTACGGATAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                            <  1:211727/147‑1 (MQ=37)
                                                                                       aaTCGGACGATTGGTGGGTAAAATGAGGGTTTTTTTGGGGGTTAGAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                        <  1:104204/151‑1 (MQ=38)
                                                                                          cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                      >  1:585187/1‑150 (MQ=255)
                                                                                            gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                    >  2:60106/1‑150 (MQ=255)
                                                                                               gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                >  1:1087782/1‑151 (MQ=255)
                                                                                                 taggggggTAAAATGAGGGTTTTTTTGGGGGTTCCAGTAAATTCCAAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                     <  1:461889/139‑1 (MQ=25)
                                                                                                  tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                             >  2:761135/1‑151 (MQ=255)
                                                                                                   aGTGGGGGAAATGAGGGGTTATTTGGGGGAGACAGGTAAATTCCAGGAAGAAAAAAACTCCGCCGGAGGCGGGGGAGAACGTTAATGACAAATTGAGCGTTATGTTGCCGGGTGCGGAGTAAGCGCATTATCCGGCCTATAAAAAtgtg                                              >  2:85365/1‑149 (MQ=255)
                                                                                                    ggggggaaaaTGGGGGTTTATTGGGGTGTAACGGGAAAATTCCAGGGAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                           <  2:1115888/144‑1 (MQ=25)
                                                                                                    gggggTGAAATGGGGGTTTTTTGGGGGTTTACAGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                           <  2:1002611/149‑1 (MQ=37)
                                                                                                     tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                           >  1:582275/1‑150 (MQ=255)
                                                                                                     ggggTAAAAGGGGGGATTTTTGGGGGATTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                              <  2:192203/146‑1 (MQ=25)
                                                                                                      gggaaaaaaggggggTTATTTGGGGTTTAGAGGAAAATCCCGGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                           <  2:1073020/138‑1 (MQ=37)
                                                                                                      gggTAAAAGGAGGGTTTTTTTGGGGGTCAAAGATAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                      <  1:294512/138‑1 (MQ=25)
                                                                                                         tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                       >  2:466099/1‑150 (MQ=255)
                                                                                                           aaagggggggTTTTTTGGGGGTAAGAGGAAAATCCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACTCCTTATCCGGCCTACaaaaa                                                  <  1:1011171/131‑1 (MQ=37)
                                                                                                           aaaTGAGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                      <  2:98986/149‑1 (MQ=255)
                                                                                                           aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                    >  2:260017/1‑151 (MQ=255)
                                                                                                           aaaTGAGGGGTTTTTGGGGGTTAACAGAAAAATTCAAGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                    <  2:302542/151‑1 (MQ=25)
                                                                                                            aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                    >  2:11501/1‑150 (MQ=255)
                                                                                                            aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                    >  1:475441/1‑150 (MQ=255)
                                                                                                             aTGAGGGGTTCTTGGGGGGTAACAGAAAAATTCAACAAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCGCGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                      <  2:257763/147‑1 (MQ=25)
                                                                                                              tGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                  >  1:135577/1‑150 (MQ=255)
                                                                                                               gAGGGTTTATTGGGGGTTTACAGGAAAAGTCCAGGCGAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                  <  2:701890/149‑1 (MQ=37)
                                                                                                                 ggggTTTTTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                               <  2:959610/150‑1 (MQ=37)
                                                                                                                  ggttttttttGGGGTTTAAGGAAAAATCCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                              <  1:218380/145‑1 (MQ=21)
                                                                                                                  gggTTTTTTTGGGTTTAAAAGAAAAATCCGCGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                              <  1:204888/151‑1 (MQ=21)
                                                                                                                        ttggggggTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                        <  2:1159033/147‑1 (MQ=255)
                                                                                                                          tGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                     >  1:897002/1‑151 (MQ=255)
                                                                                                                           tggggTTATCAGGTAATTTCCAGGGAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                      <  1:277866/149‑1 (MQ=255)
                                                                                                                           ggggtTTACAGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                        <  2:975943/147‑1 (MQ=37)
                                                                                                                           gggggTTACAGTAAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                     <  1:523914/150‑1 (MQ=37)
                                                                                                                                 ttACGGGTAAATCCCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                <  2:824418/150‑1 (MQ=37)
                                                                                                                                    aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAATGCAATTCAACTTGTAGGc             >  1:1013416/1‑151 (MQ=255)
                                                                                                                                      cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc            >  2:984023/1‑151 (MQ=255)
                                                                                                                                       aGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt           >  1:1062405/1‑151 (MQ=255)
                                                                                                                                          gTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg          >  2:48433/1‑150 (MQ=255)
                                                                                                                                          gTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg          >  1:1172298/1‑150 (MQ=255)
                                                                                                                                            aaaTTCCAGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg          <  1:622231/148‑1 (MQ=255)
                                                                                                                                                     aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  >  2:28436/1‑149 (MQ=255)
                                                                                                                                                     aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTCACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  >  1:201275/1‑149 (MQ=255)
                                                                                                                                                          |                                                                                                                                                
ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC  >  NC_000913/2001763‑2002054

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