Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,001,918	0	A	C	41.9%	123.9 / 12.5	74	intergenic (‑129/+192)	fliA/fliC	RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (37/6);  new base C (0/31);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 6.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001770‑2002066                                                                                                                                                           |                                                                                                                                                      gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑c                                                                                                                                                     >  2:158859/1‑150 (MQ=255)   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑cccc                                                                                                                                                  >  2:1144799/1‑151 (MQ=255)      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCccgcc                                                                                                                                               >  2:1155231/1‑151 (MQ=255)      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCccgc                                                                                                                                                <  1:1059090/150‑1 (MQ=255)              aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCgg                                                                                                                                       >  1:544615/1‑151 (MQ=255)               aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCggg                                                                                                                                      >  2:1083861/1‑151 (MQ=255)                ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                                     >  1:1002465/1‑151 (MQ=255)                  cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGGTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                                    >  1:831640/1‑150 (MQ=255)                     gaaaaaCGGCCCTGTTTTATAGGATTTATGGGAATGTTTATCCTTTTCTGCGGGGTTTTTAATCGGAGGTTAGGGGGGTGAAATGAGGGTTTTTTTGGGGT‑TT‑A‑AG‑GATAAATT‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGc                                                                                                                                 <  1:490696/147‑1 (MQ=255)                      aaaaaGAGCCCGGTGTTATATGATTAAGCGATATTTTTACCCTTTTCGGGGCGTTTTTTAACCGGACTTTTAGTGGGTGAAATGAGGGGTAGTGGGGGGG‑TG‑A‑CA‑GCTAAATT‑CA‑AGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                               <  2:42956/149‑1 (MQ=255)                        aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                             >  1:159396/1‑151 (MQ=255)                              cccTGCTTTATAGGATTTATGGGAATATTTATCCTTTTCTGCGGGTTTTTAAATCGGACGATTGGTGGGTAAAAGGGGGGGTTTTTTGGGGG‑TA‑A‑CA‑GAAAATTC‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                          <  1:79841/148‑1 (MQ=255)                                       atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                              >  2:675943/1‑151 (MQ=255)                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                            >  2:159639/1‑151 (MQ=255)                                          tgttttACCGGCATGATTGTCCGTTTCGGCAGGTTTTTTAAGCGGACATTTGGTGGGT‑AAAAGGGGGGTTATTTGGGGG‑GT‑T‑CAGGGAAAATT‑CC‑AGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                <  2:287932/142‑1 (MQ=255)                                              ttATGGGCTTGATTTCCCTTTTCTGCGGGGTTTTTAACCGGACGTTGAGGGGGTAAAATGGGGGTTTATTGGGGGT‑TT‑A‑CA‑GGTAAATT‑CC‑GGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                        <  2:695177/150‑1 (MQ=255)                                               tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                      >  2:754276/1‑151 (MQ=255)                                               tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCAAGTTGCTGACAAATTGCGCTTTAt                                                                                                        >  1:315964/1‑149 (MQ=255)                                                   ggaaTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTAAAATGAGGGTTTATTGGGGGT‑TT‑A‑AA‑GGTAAATT‑CC‑AGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGATTTATGTTg                                                                                                    <  2:453514/146‑1 (MQ=255)                                                   ggCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                  >  1:416816/1‑151 (MQ=255)                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACACATTGCGCTTTATGTTGcc                                                                                                  >  1:693116/1‑150 (MQ=255)                                                     cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                >  2:212574/1‑151 (MQ=255)                                                          ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                           >  2:677616/1‑151 (MQ=255)                                                          ttATCCGTTTCGGCAGGTTTTTTAATCGGAGGATTAGTGGGTGAAATGAGGGTTTATTGGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                              <  2:254413/148‑1 (MQ=255)                                                            aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                         >  1:251858/1‑151 (MQ=255)                                                            aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                         >  1:1185665/1‑151 (MQ=255)                                                                 tttCGGCTGGTTTTTTAATCGGACGATTAGGGGGTGAAAGGGGGGTTTTTTTGGGGT‑TT‑A‑CA‑GGAAAATC‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                         <  1:687240/146‑1 (MQ=255)                                                                      gcgggttttttAATAGGACGATTGGTGGGAAAAAGGAGGGGTTTTTGGGGGT‑TT‑A‑CA‑GATAATTC‑CC‑AGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                                <  2:415052/144‑1 (MQ=255)                                                                      gcggggTTTTTAATCGGACGATTAGGGGGTGAAAGGAGGGTTTATTGGGGGT‑TT‑A‑CA‑GGTAAATT‑CC‑AGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                               <  2:307954/148‑1 (MQ=255)                                                                         ggttttttAATCGGAGATTTGGGGGGTGAAAGGAGGGTTTATTAGTGGG‑TA‑A‑AA‑GTTAAATT‑CA‑AACAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                            <  2:11924/148‑1 (MQ=255)                                                                         ggttttttAATCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTTGGGGG‑TT‑A‑‑A‑GGGAAATTCCC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                   <  1:1143985/141‑1 (MQ=255)                                                                             ttttAATCGGACGATTAGTGGGTAAAAGGAGGGTTTATTTGGGGG‑TT‑A‑CA‑GGAAAATT‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         <  2:754946/150‑1 (MQ=39)                                                                              tttAATCGGACGTTTAGTGGGTAAAATGGGGGTTTATTTGGGGT‑TT‑A‑CA‑GGAAAATC‑CC‑AGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        <  2:320852/150‑1 (MQ=38)                                                                              ttaaaTGGGATGATTAGGGGGAGAAATGGGGGTTTTTTTGGGGG‑TT‑A‑CG‑GATAAATC‑CC‑GGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        <  1:211727/147‑1 (MQ=37)                                                                                 aaTCGGACGATTGGTGGGTAAAATGAGGGTTTTTTTGGGGG‑TT‑A‑GA‑GGTAAATC‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                    <  1:104204/151‑1 (MQ=38)                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                  >  1:585187/1‑150 (MQ=255)                                                                                      gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                                >  2:60106/1‑150 (MQ=255)                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  1:1087782/1‑151 (MQ=255)                                                                                           taggggggTAAAATGAGGGTTTTTTTGGGGG‑TT‑C‑CA‑‑GTAAATT‑CC‑AAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                                 <  1:461889/139‑1 (MQ=25)                                                                                            tAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         >  2:761135/1‑151 (MQ=255)                                                                                             aGTGGGGGAAATGAGGGGTTATTTGGGGG‑AG‑A‑CA‑GGTAAATT‑CC‑AGGAAGAAAAAAA‑CTCCGCCGGAGGCGGGGGAGAACGTTAATGACAAATTGAGCGTTATGTTGCCGGGTGCGGAGTAAGCGCATTATCCGGCCTATAAAAAtgtg                                                          >  2:85365/1‑149 (MQ=255)                                                                                              ggggggaaaaTGGGGGTTTATTGGGGTG‑TA‑A‑CG‑GGAAAATT‑CC‑AGGGAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       <  2:1115888/144‑1 (MQ=25)                                                                                              gggggTGAAATG‑GGGGTTTTTTGGGGGTTT‑A‑CA‑GGAAAATC‑CC‑GGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       <  2:1002611/149‑1 (MQ=37)                                                                                               tGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       >  1:582275/1‑150 (MQ=255)                                                                                               ggggTAAAAGGGGGGATTTTTGGGGGA‑TT‑A‑CA‑GGAAAATC‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                          <  2:192203/146‑1 (MQ=25)                                                                                                gggaaaaaaggggggTTATTTGGGGT‑TT‑A‑GA‑GGAAAATC‑CC‑GGGAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       <  2:1073020/138‑1 (MQ=37)                                                                                                gggTAAAAGGAGGGTTTTTTTGGGGG‑TC‑A‑AA‑GATAATTC‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                  <  1:294512/138‑1 (MQ=25)                                                                                                   tGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                   >  2:466099/1‑150 (MQ=255)                                                                                                     aaagggggggTTTTTTGGGGG‑TA‑A‑GA‑GGAAAATC‑CC‑AGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACTCCTTATCCGGCCTACaaaaa                                                              <  1:1011171/131‑1 (MQ=37)                                                                                                     aaaTGAGGGTTTATTGGGGGG‑TT‑A‑CA‑GGAAAATT‑CC‑AGGCAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                  <  2:98986/149‑1 (MQ=255)                                                                                                     aaaTGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                >  2:260017/1‑151 (MQ=255)                                                                                                     aaaTGAGGGGTT‑TTTGGGGG‑TTAA‑CA‑GAAAAATT‑CA‑AGGAAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                <  2:302542/151‑1 (MQ=25)                                                                                                      aaTGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                >  2:11501/1‑150 (MQ=255)                                                                                                      aaTGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                >  1:475441/1‑150 (MQ=255)                                                                                                       aTGAGGGGTTCTTGGGGGG‑TA‑A‑CA‑GAAAAATT‑CA‑ACAAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCGCGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                  <  2:257763/147‑1 (MQ=25)                                                                                                        tGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              >  1:135577/1‑150 (MQ=255)                                                                                                         gAGGGTTTATTGGGGGT‑TT‑A‑CA‑GGAAAAGT‑CC‑AGGCGAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              <  2:701890/149‑1 (MQ=37)                                                                                                           ggggTTTTTGGGGGG‑TT‑A‑CA‑GGTAAATC‑CC‑AGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                           <  2:959610/150‑1 (MQ=37)                                                                                                            ggttttttttGGGG‑TT‑T‑AA‑GGAAAAAT‑CCGAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          <  1:218380/145‑1 (MQ=21)                                                                                                            gggTTTTTTTGGGT‑TT‑A‑AA‑AGAAAAAT‑CCGCGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          <  1:204888/151‑1 (MQ=21)                                                                                                                  ttgggggg‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    <  2:1159033/147‑1 (MQ=255)                                                                                                                    tGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                 >  1:897002/1‑151 (MQ=255)                                                                                                                     tgggg‑TT‑ATCA‑GGTAATTT‑CC‑AGGGAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                  <  1:277866/149‑1 (MQ=255)                                                                                                                     ggggt‑TT‑A‑CA‑GGAAAATC‑CC‑GGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    <  2:975943/147‑1 (MQ=37)                                                                                                                     ggggg‑TT‑A‑CA‑GTAAAATC‑CC‑AGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                 <  1:523914/150‑1 (MQ=37)                                                                                                                           tt‑A‑CG‑GGTAAATC‑CC‑GGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                            <  2:824418/150‑1 (MQ=37)                                                                                                                              a‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAATGCAATTCAACTTGTAGGc                         >  1:1013416/1‑151 (MQ=255)                                                                                                                                cA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                        >  2:984023/1‑151 (MQ=255)                                                                                                                                 a‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                       >  1:1062405/1‑151 (MQ=255)                                                                                                                                    gTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                      >  2:48433/1‑150 (MQ=255)                                                                                                                                    gTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                      >  1:1172298/1‑150 (MQ=255)                                                                                                                                      aaaTT‑CC‑AGGAAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                      <  1:622231/148‑1 (MQ=255)                                                                                                                                             cagggaaaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag                 <  1:954424/140‑1 (MQ=255)                                                                                                                                               aGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  2:28436/1‑149 (MQ=255)                                                                                                                                               aGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTCACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  1:201275/1‑149 (MQ=255)                                                                                                                                                ggaaaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag            <  2:62638/145‑1 (MQ=255)                                                                                                                                                ggaaaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  2:849131/144‑1 (MQ=255)                                                                                                                                                         aaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  2:410778/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                      GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CA‑GGTAAATT‑CC‑AGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001770‑2002066

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: ‑