Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,486	0	C	A	39.7%	100.5 / 43.0	73	intergenic (+291/‑165)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (41/3);  new base A (0/29);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 2.61e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29338‑29635                                                                                                                                                         |                                                                                                                                                      tttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCc‑aca                                                                                                                                                      >  1:305166/1‑150 (MQ=255)  ttgttgTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACaa                                                                                                                                                     <  2:523017/150‑1 (MQ=255)      tgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAaata                                                                                                                                                 >  1:733225/1‑150 (MQ=255)                tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATaaaaaa                                                                                                                                       <  2:4639/150‑1 (MQ=255)                  aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATcc                                                                                                                                    >  2:254535/1‑151 (MQ=255)                   aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATccc                                                                                                                                   >  2:953965/1‑151 (MQ=255)                        gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCa                                                                                                                               <  1:387725/150‑1 (MQ=255)                        gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGAGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGGC‑ACAAAATAATATAAAAAATCCCGCCAt                                                                                                                              >  2:1244008/1‑151 (MQ=255)                           cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                            >  2:1131437/1‑150 (MQ=255)                                ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                       >  1:530305/1‑150 (MQ=255)                                ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                       >  1:207770/1‑150 (MQ=255)                                  tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                    >  2:90438/1‑151 (MQ=255)                                        tttttCTCCTTGTTTTTTTTTCAGACATTCTTCTTCCTGCGCAA‑CGTTTTCCAAGAACAGGTTAGATGACTTTTTTGCCACTTAGTGCTTGT‑AAAACTGGAAT‑GGCCA‑AAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAg                                                                                                               <  1:1131514/150‑1 (MQ=255)                                         ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                             >  2:816132/1‑151 (MQ=255)                                         ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                             >  1:603570/1‑151 (MQ=255)                                          tttCTGCTTTTTTTTTTTTCAGGCATTTTTTTTGCTGCGAAACCGTTTTCA‑AGAAAAGTTTAGATGTTTTTTTTGTCGTTTAAGGCTTTA‑AAAACATGCAG‑GAGCC‑AAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAgcg                                                                                                             <  1:1172443/150‑1 (MQ=255)                                             cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAACTCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                          >  1:706101/1‑150 (MQ=255)                                             cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                          >  2:177697/1‑150 (MQ=255)                                              tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCTCAAGCGTTTTAC‑AGAACAGGTTAGATGATCTTTTTTTTGCTTAAGGCGTGT‑AAAGCCTGCAT‑GAGCC‑ACAATATTATATATAAAATCCAGCTATAAAGTTGATTTTTAgtgtc                                                                                                           >  2:980272/1‑144 (MQ=255)                                               gCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                        >  2:171154/1‑150 (MQ=255)                                                   gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATTTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   >  1:247606/1‑151 (MQ=255)                                                    ttttttttCCATTCATTTTTCTTGGCGCGCAA‑CTTTTTCCAAAAACAGGTAAAATAATTTTTTTGTCGCTTAATCCCTAA‑AAAACAGGCAT‑GACCC‑AAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   <  2:1172536/150‑1 (MQ=255)                                                    tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCA‑AGAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTGA‑AAAACATGCAT‑GAGCC‑AAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                   <  1:767664/150‑1 (MQ=255)                                                      ttttttTAATACATTCTTCTTGCGCAGCAGGCTTTTACA‑AAAGCAGGTAAAATAATTTTTTTTTCGCTAAATGCCTAA‑AAATCATGCAG‑GCGCC‑AAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                 <  2:227745/150‑1 (MQ=255)                                                        tATTTCATCCAATCTTCTTCCTGCACAAGCTTTTTCA‑AGAAAAGTTTAGATAATTTTTTTGTCTTTTACTCCTTAA‑AAAACATGAAG‑GAGCC‑AAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                <  2:207764/149‑1 (MQ=255)                                                         aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                             >  2:465045/1‑151 (MQ=255)                                                         aTTTCATGCAATCTTCCTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                             >  1:1090086/1‑151 (MQ=255)                                                          tttCATGCATTTTTCTTGCTGCGAAAGCTTTTCCC‑AGAAAAGGTTAGATGTTTTTTTTGTCGTTAAATCCCTGA‑AAAACATGCAT‑GACCA‑AAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                             <  1:568732/150‑1 (MQ=255)                                                            tCATGCAATCTGCTTGCTGCGCAAGCGTTTTGC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                          >  2:994551/1‑151 (MQ=255)                                                              aTGCATTTTTCTTGCTGCGCAAGCTTTTCCC‑AAAACAGTTTAGATGTTTTTTTTGTCGTTTAATGCCTGA‑AAAACATGCAT‑GAGCC‑AAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                           <  1:1062071/148‑1 (MQ=255)                                                                    tcttctTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                  >  1:439972/1‑151 (MQ=255)                                                                    tcttctTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                  >  1:1237224/1‑151 (MQ=255)                                                                      ttctgGCTGCGGAAGGTTTTTCG‑ATAAATGTTTGGATGCTTTTTTTGCCG‑TTAAGGCCTTTAAAAACGTGCGG‑AACCC‑AAAAAATAAAAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                 <  1:926741/150‑1 (MQ=255)                                                                      ctcatgccgcGCAAGCGTTTTCG‑AGACCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAAAAGGCAT‑GAGCC‑ACAAACTAATATAAAAAATCCCGCCATTAAGTTGTCTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  2:1238447/5‑151 (MQ=255)                                                                       tcttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  1:339745/1‑150 (MQ=255)                                                                        cttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  1:726423/1‑149 (MQ=255)                                                                        cttCCTGCGCAAGCGTTATCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                >  2:192462/1‑149 (MQ=255)                                                                         ttGCTGCGAAAGCTTTTCCA‑AAAACAGGTTAGATGATTTTTTTGCCGCTTATTGCCTGA‑AAAACATGCAG‑GACCC‑AAAAAAAAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                  <  1:99181/146‑1 (MQ=255)                                                                           gctgcGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                            >  2:452813/1‑150 (MQ=255)                                                                               cgcAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                       >  1:946143/1‑151 (MQ=255)                                                                                  aggctttttCC‑AAAAGAGTTTAGATAATTTTTTTGTCGTTAAATCCCTGA‑AAACCAGGCAG‑GACCC‑AAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                     <  1:472801/145‑1 (MQ=255)                                                                                  aaGCTTTTTCC‑AAAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGA‑AAAACATGCAG‑GACCC‑AAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                          <  1:78544/145‑1 (MQ=255)                                                                                      ttttcCC‑AAAACAGTTTAGATAATTTTTTTGCCGCTTAATCCCTAA‑AAAACATGAAT‑GAGCC‑AAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                        <  1:155156/142‑1 (MQ=255)                                                                                        tttCC‑AGAAAAGGTTAGATTTTTTTTTTG‑AGGTTAAAGGCGTA‑AAAAAATGCATCGCCCC‑AAAAAAAAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                      <  1:268917/143‑1 (MQ=255)                                                                                         ttCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                              >  2:909680/1‑150 (MQ=255)                                                                                             aaaaacggtttAGATGTTCTTTTTGCCGTTTAAGGCGTGA‑AAAACATGGAT‑‑AGCCAAAAAAATAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAt                                                          <  1:113013/142‑1 (MQ=255)                                                                                              aaaaCAGGTTAGATGTTTTTTTTGTCGTTAAATCCTTAA‑AAAAAATGAAT‑GACCA‑AAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                           <  1:357071/147‑1 (MQ=255)                                                                                              aGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                         >  2:786533/1‑151 (MQ=255)                                                                                               gAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                        >  1:607516/1‑151 (MQ=255)                                                                                                aaCAGGTTAGATTTTTTTTTTGTCGTTAAATGCTTAA‑AAAACATGCAT‑GAGCC‑AAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                       <  1:515995/151‑1 (MQ=255)                                                                                                aaCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCA‑ACAAAATAATATAAAAAATCACGCCATTAAGTTGACTTTTAGCGCACATATCTCCAGAATGCCGCCGTTTGCCAGAAGTTCGTCGGTAAGCAGATTTg                                                       >  2:1017402/1‑151 (MQ=255)                                                                                                   aGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                    >  2:732020/1‑151 (MQ=255)                                                                                                        aGATGTTTTTTTTGTCATTAAATGCTTTA‑AAAACAGGAAT‑GAGCC‑AAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                     <  1:155164/145‑1 (MQ=255)                                                                                                           tgatCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                            >  2:653350/1‑151 (MQ=255)                                                                                                              tttttttGTCGTTAAATGCTTAA‑AAAACATGCAG‑GACCC‑AAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                           <  1:1021123/147‑1 (MQ=255)                                                                                                                tttttGTCTTTAAGTGCTTGA‑AAAAAAGGAAT‑GACCA‑AAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                          <  1:70822/148‑1 (MQ=255)                                                                                                                 ttttCTCGTTAAATCCCTAA‑AAAACTTGTAT‑GGCCC‑AAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                       <  2:262313/150‑1 (MQ=255)                                                                                                                  tttGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                      >  2:785494/1‑150 (MQ=255)                                                                                                                   ttGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                    >  1:241969/1‑151 (MQ=255)                                                                                                                            aaTGCCTGA‑AAAACATGCAT‑GAGCC‑AAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                             <  1:954033/149‑1 (MQ=255)                                                                                                                              tGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATGGATTTACGTCATCATTGTGAATTAat                         >  1:1008768/1‑151 (MQ=255)                                                                                                                               gCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGCCTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                         >  1:1281956/1‑150 (MQ=255)                                                                                                                                 cctaaaaaaaaTGCAT‑GACCA‑AAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                       <  2:185720/146‑1 (MQ=255)                                                                                                                                 cTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                      >  2:414572/1‑151 (MQ=255)                                                                                                                                  tGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                     >  1:866756/1‑151 (MQ=255)                                                                                                                                        aaCATGCAT‑GACCC‑AAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                <  1:228197/151‑1 (MQ=255)                                                                                                                                         aCAGGCAG‑GAGCA‑AAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa               <  1:973372/151‑1 (MQ=255)                                                                                                                                           atgcat‑gAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                >  1:129476/1‑148 (MQ=255)                                                                                                                                            tgcat‑gAGCC‑AAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag             <  1:380520/150‑1 (MQ=255)                                                                                                                                               at‑gAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg           >  2:661822/1‑149 (MQ=255)                                                                                                                                               at‑gAGCC‑AAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag         <  1:171158/151‑1 (MQ=255)                                                                                                                                                acgaccc‑AAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag         <  2:430833/146‑1 (MQ=255)                                                                                                                                                    gCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga      >  1:401482/1‑150 (MQ=255)                                                                                                                                                     cc‑AAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt            <  1:54091/143‑1 (MQ=255)                                                                                                                                                        acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat  >  2:311336/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                      TTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGT‑AAAACATGCAT‑GAGCC‑ACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29338‑29635

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