Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,496	0	T	A	42.2%	91.4 / 40.1	71	intergenic (+301/‑155)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (29/12);  new base A (0/30);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.20e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTG  >  NC_000913/29348‑29641                                                                                                                                                           |                                                                                                                                                  tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAAtata                                                                                                                                                  <  2:194793/150‑1 (MQ=255)  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATaaa                                                                                                                                                <  2:339005/151‑1 (MQ=255)   aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATaaa                                                                                                                                                >  2:600307/1‑150 (MQ=255)      tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATaaaaaa                                                                                                                                             <  1:39954/150‑1 (MQ=255)       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAAt                                                                                                                                            <  1:437955/150‑1 (MQ=255)       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATc                                                                                                                                           <  2:1142935/151‑1 (MQ=255)       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATc                                                                                                                                           <  2:1047748/151‑1 (MQ=255)       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATc                                                                                                                                           <  1:32045/151‑1 (MQ=255)        aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATcc                                                                                                                                          <  1:389348/151‑1 (MQ=255)        aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATc                                                                                                                                           >  1:1238031/1‑150 (MQ=255)          ttttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCg                                                                                                                                        >  2:1009328/1‑151 (MQ=255)                         ccttttttttCTGCTCTTTTTTATTTCATGAATTCTTCTTGCTGCGCAAGCTTTT‑TCCAA‑AA‑CGGT‑TT‑AGATGATTTTTTTGTCTCTTAATGCCTGAAAAACATGCATGACC‑CAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACtt                                                                                                                         <  2:452908/147‑1 (MQ=255)                          cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                        >  2:96297/1‑151 (MQ=255)                             tttttccgtctctTTTTTTTTTCACCCATTCTTTTTCCTGCGAATGCTTTT‑CCAAA‑AA‑CAGT‑TT‑AGATAATTTTTTTGTCGATATGGTCCTTTAAAACATGCATGAGC‑AAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                      <  1:1155432/141‑1 (MQ=255)                               ttttcttttcTTTTTTATTGCATGCATTCTTTTTGCT‑CGCAGGCGTTTCGCCAG‑AG‑CGGG‑TT‑ATATTTTTTTTTTATCGATTACCCGCTG‑AAAACCTTTTTGGCCGCAGAAAAAAATAAAAAAAGTCCCGCCATTAAGTTGACTTTTAg                                                                                                                     <  1:55618/149‑1 (MQ=255)                               ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                                   >  2:88485/1‑151 (MQ=255)                                tttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                                  >  1:703727/1‑151 (MQ=255)                                   cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                >  2:91555/1‑150 (MQ=255)                                    tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                              >  2:556789/1‑151 (MQ=255)                                          ttttttttTCATGAATTTTTTTTGCTGCGAAGGCGTTT‑CCAAA‑AA‑AAGG‑TA‑AAATTTTTTTTTTTTCGTTTATTCCTTAAAAAACATACTTGCGC‑AAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                          <  1:1281414/149‑1 (MQ=255)                                           ttttttttCATGCATTTTTTTTGCTGCGAAAG‑TTTT‑TCCAA‑AA‑AAAG‑TTAAAATAATTTTTTTGCCGCTTAATGCCTAAAAAACATGCTTGACC‑CAAAAAATAAAAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                        <  1:131180/150‑1 (MQ=255)                                            tttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGTTTAATGCCTGTAAAACATCCATGAGC‑CACAAATTAATTTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                         >  2:35014/1‑148 (MQ=255)                                              tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                     >  2:369876/1‑150 (MQ=255)                                               aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                   >  2:93338/1‑151 (MQ=255)                                                     tGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                             >  1:724993/1‑151 (MQ=255)                                                          tcttctTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                        >  2:154336/1‑151 (MQ=255)                                                            ttctgGCTGCGCAGGC‑TTT‑TCCAA‑AA‑AAGGTTT‑AGATGTTTTTTTTGTCTTTTAATGCTTGAAAAACAGGAGGGAGC‑CAAAAAAAAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                           <  2:14223/146‑1 (MQ=255)                                                                 gctgcGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                  >  1:1063668/1‑150 (MQ=255)                                                                   ctgcgaacgcttt‑tCCAA‑AACCAGT‑TG‑AGATAATTTTTTGGTCGCTAAAGCCCTGAAAACCTTCCTTGCCC‑CAAAAAATAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                  <  1:485658/139‑1 (MQ=255)                                                                    gcgaaggttttt‑tCCAA‑AA‑CAGT‑TT‑AGATGATTTTTTTGTCCATAAACGCCTGAAAAACTTGCATGACC‑CAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                <  2:914576/140‑1 (MQ=255)                                                                        aaGTGTTT‑TCCAA‑AA‑CGGA‑AA‑AGATTATTTTTTTGGCCTTTAATGCCAGTAAAACTGGAATGAGC‑AAAAAAAAAATAAAAAAAAGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                             <  2:36132/148‑1 (MQ=255)                                                                        aaGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATTTTTTTGTCTCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                           <  2:1195942/150‑1 (MQ=255)                                                                           gcttt‑tCCAAAAA‑CAGT‑TT‑AAATTTTTTTTTTTTCGATTAATCCCTAAAAAACTTGAATGAAC‑CAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                        <  1:356567/149‑1 (MQ=255)                                                                             ttt‑tCCAG‑AC‑AAGG‑TT‑TGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGACAGAAATTCGTCGGt                                                                        >  1:921124/1‑148 (MQ=255)                                                                                 tCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                    >  2:1014954/1‑149 (MQ=255)                                                                                       aa‑CAGT‑TT‑AAATGTTCTTTTTGTCCTTTAACCCCTAAAAAACTTGATTGATC‑AACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                              <  1:939558/150‑1 (MQ=255)                                                                                         aacggttt‑AGATGATTTTTTTTTCGTTAAATCCCTAAAAAACATGCATGACC‑CAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCTAGAAATTCGTCGGTAAGCAGATTTg                                                             <  2:362034/145‑1 (MQ=255)                                                                                                t‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                                      >  2:726831/1‑151 (MQ=255)                                                                                                  aaattatCTTTTTGTCCTTTATTGCTTAAAAAACGTCCATGAGC‑CAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                        <  2:941423/143‑1 (MQ=255)                                                                                                     ttttttttttGTCTTTTAATCCCTATAAAACAGGCATGACC‑CAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                     <  2:1064747/143‑1 (MQ=255)                                                                                                     tgatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                  >  2:1246987/1‑151 (MQ=255)                                                                                                        tttttttCTCGTTTAATCCTTAAAAAACATGAAGGACC‑CAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                     <  2:932958/143‑1 (MQ=255)                                                                                                          ttttttTCGATTAATGCCTAAAAAACATGCAGGAGC‑AAAAAAAAAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                              <  1:385353/150‑1 (MQ=255)                                                                                                           ttttGTCGTTTATTGCCTAAAAAACTTGCTGGACC‑CAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                             <  1:825156/150‑1 (MQ=255)                                                                                                             ttGTCCCTAAATGCCTAAAAAACGTGCATGACC‑CAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                           <  2:973869/150‑1 (MQ=255)                                                                                                               gtcctttATTGCCAGAAAAACAGACAGGACC‑AAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                         <  1:1208463/145‑1 (MQ=255)                                                                                                                 cGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                      >  2:1255628/1‑151 (MQ=255)                                                                                                                  gctcagtgccgaaaaaaTATTCATGAGC‑CAAAAAATAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                      <  1:490868/137‑1 (MQ=255)                                                                                                                  gCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                      >  2:919027/1‑150 (MQ=255)                                                                                                                     tAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                  >  1:1255138/1‑151 (MQ=255)                                                                                                                      aaTGCCTAAAAAACATGCAGGAGC‑AAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                  <  2:812343/150‑1 (MQ=255)                                                                                                                          cttaaaaaaCAGGCATGACC‑AAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                                <  1:1030031/143‑1 (MQ=255)                                                                                                                          cctaaaaaaCATGTAGGACC‑CAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                              <  1:496102/145‑1 (MQ=255)                                                                                                                           cTGAAAAACATGCATGAGC‑CAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                               <  2:179381/148‑1 (MQ=255)                                                                                                                            tgaaaaaCAGGCATGAGC‑AAAAAAAAAATATAAAAAGCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                   <  2:387444/140‑1 (MQ=255)                                                                                                                              aaaaacggacatgAGC‑CAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                             <  1:458090/138‑1 (MQ=255)                                                                                                                               aaaaCATGAATGAGC‑CAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                          <  2:593666/149‑1 (MQ=255)                                                                                                                                  aCATGCTGGACC‑CAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                     <  1:422254/151‑1 (MQ=255)                                                                                                                                  aCATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                      >  1:431650/1‑150 (MQ=255)                                                                                                                                      gcatgAGC‑CACAAAATAATATAAAAAATCCTGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                 >  1:999288/1‑151 (MQ=255)                                                                                                                                        atgAGC‑CAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga                <  1:369862/150‑1 (MQ=255)                                                                                                                                           aaccaaaaaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag               <  2:159279/142‑1 (MQ=255)                                                                                                                                           aGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg             >  2:239590/1‑150 (MQ=255)                                                                                                                                            cc‑cAAAAAATAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt              <  1:232188/147‑1 (MQ=255)                                                                                                                                             c‑cACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat          >  2:964666/1‑151 (MQ=255)                                                                                                                                             c‑cAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa           <  2:558173/150‑1 (MQ=255)                                                                                                                                                acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat        >  2:181741/1‑151 (MQ=255)                                                                                                                                                aaaaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg             <  2:733622/144‑1 (MQ=255)                                                                                                                                                 aaaaaTAAAAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt       <  1:11371/150‑1 (MQ=255)                                                                                                                                                  aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc      >  2:228837/1‑151 (MQ=255)                                                                                                                                                    aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc    <  2:771392/151‑1 (MQ=255)                                                                                                                                                    aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc      >  1:1218289/1‑149 (MQ=255)                                                                                                                                                     aaaataAAAAAAAGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc      <  2:1218356/146‑1 (MQ=255)                                                                                                                                                      taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg  >  2:522144/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                  TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTT‑TCCAG‑AA‑CAGG‑TT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGC‑CACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTG  >  NC_000913/29348‑29641

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