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

TTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGA  >  NC_000913/29339‑29631                                                                                                                                                       |                                                                                                                                                  ttgttgTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACaa                                                                                                                                                 <  1:1117146/150‑1 (MQ=255)  tgttgtTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACaaa                                                                                                                                                <  2:9915/150‑1 (MQ=255)   gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACaaaa                                                                                                                                               <  2:900057/150‑1 (MQ=255)    ttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGG‑TAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAaata                                                                                                                                             <  2:1107891/150‑1 (MQ=255)          tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                        >  2:1011189/1‑150 (MQ=255)            aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                     <  2:210510/151‑1 (MQ=255)             tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                     >  2:366063/1‑150 (MQ=255)             tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                     >  2:660046/1‑150 (MQ=255)              gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                    <  2:908435/150‑1 (MQ=255)               tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                  >  1:196397/1‑151 (MQ=255)                  aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                               <  2:195585/151‑1 (MQ=255)                    tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                              <  2:388655/150‑1 (MQ=255)                    tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                             >  2:672880/1‑151 (MQ=255)                         ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTTC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                         >  1:47853/1‑150 (MQ=255)                            tttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGt                                                                                                                     >  1:249953/1‑151 (MQ=255)                                 tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                >  1:421210/1‑151 (MQ=255)                                      ttttttttGCTCTTTTTTATTTCATAAATTCTTCTTGCGGCACAAGGTTTTTCA‑AGAACGGGTTAGATGTTTTTTTTGTCGTTAAATGCC‑TG‑AAAAACAGGAAGGGGCCAAAAAATAAAATAAAAAACCCCGCCATTAAGTTGACtttt                                                                                                             <  2:1030354/149‑1 (MQ=255)                                      ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                             >  1:556116/1‑149 (MQ=255)                                      ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                             >  1:587968/1‑149 (MQ=255)                                       tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                          >  1:835061/1‑151 (MQ=255)                                        ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                         >  1:664848/1‑151 (MQ=255)                                         tttCTGCTCTTTTTTATTTCATGCATTCTTCTTGCGGCGCAAGCTTTTTCC‑AGAACAGTTTAGATGATTTTTTTGTCGTTTAATGCC‑TG‑AAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                         <  2:630108/150‑1 (MQ=255)                                          ttCTGCTTTTTTTTATTTCATGAATTCTTCTTGCTGCGAAAGCTTTTTCA‑AAAACGGTTAAGATGATTTTTTTGTCGTTTAAGCCC‑TG‑AAAAACATGCAGGGGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                        <  2:537504/150‑1 (MQ=255)                                               cTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                  >  2:712890/1‑151 (MQ=255)                                                ttttttttATTACATGCAATTTTTTTGCTGCGCAAGCGTTTCCC‑AAAACAGGTTAGATGATTTTTTTGGCGTTAAATGCT‑TG‑AAAAACATGCAGGGCCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCata                                                                                                   <  2:181362/146‑1 (MQ=255)                                                   ttttttttTCATCAATTCTTTTTGCGGCGCAAG‑GTTTTCCAAAAACAGGTAAGATATTTTTTTTGTCGATAAAT‑CC‑TGAAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCata                                                                                                   <  1:99087/146‑1 (MQ=255)                                                        aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                         >  1:709073/1‑151 (MQ=255)                                                           tCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                       >  1:864203/1‑150 (MQ=255)                                                             aTGCAATCTTTTTGCTGCGAAAGCTTTTCCA‑AAAACGGTTTAGATAATTTTTTTGTCGTTAAATGCC‑TG‑AAAAACATGAAGGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                     <  2:790185/150‑1 (MQ=255)                                                             aTGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                    >  1:1523/1‑151 (MQ=255)                                                                cAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                 >  1:4267/1‑151 (MQ=255)                                                                cAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                 >  2:867173/1‑151 (MQ=255)                                                                 aaTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                >  2:1075882/1‑151 (MQ=255)                                                                 aaTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                >  2:19383/1‑151 (MQ=255)                                                                    cttcttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                             >  1:417608/1‑151 (MQ=255)                                                                     tctttGCTGCGAAGGCTTTTCCA‑AAAACATGTTAAATAATTTTTTTGCCCCTTAACCCT‑TT‑AAAAACAGTTATGAGCAAAAAAATAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                   <  1:65481/141‑1 (MQ=255)                                                                      tcttGCTGCGCAAGCGTTTTCC‑AAAACGGTTGAGATGATCTTTTTGTCGCTTAATCCC‑TA‑AAAAACATGCATGACCCACAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                               <  2:217086/147‑1 (MQ=255)                                                                       cttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                            >  2:539327/1‑149 (MQ=255)                                                                          gctgcGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                        >  1:934848/1‑150 (MQ=255)                                                                                 aaGCTTTTTCC‑AGAACGGGTTAGATAATTTTTTTGTCGTTTAATCCT‑TA‑AAAAACATGAATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                    <  2:589253/147‑1 (MQ=255)                                                                                 aaGCTTTTTCC‑AGAACGGGTAAGATGTTTTTTTTGTCGTTAAATGCT‑TG‑AAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                 <  2:200304/150‑1 (MQ=255)                                                                                 aaGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                >  2:632909/1‑151 (MQ=255)                                                                                  aGCTTTTCCC‑AAAACAGTTTAGATGATTTTTTTG‑CCTTTAATTCCTTA‑AAAAACTTGCATGACCCAAAAAATAAAAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                               <  1:672315/151‑1 (MQ=255)                                                                                  aGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                               >  1:209345/1‑151 (MQ=255)                                                                                      tttcCC‑AGAAAAGGTTAGATGTTTTTTTTGTCGCTTAATCCC‑TG‑TAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                             <  2:907189/149‑1 (MQ=255)                                                                                         tCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                          >  2:1073217/1‑149 (MQ=255)                                                                                             aaaaCAGTTTAGATGTTTTTTTTGTCGTTTAATCCC‑TA‑AAAAACATGCAGGACCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                     <  2:412035/149‑1 (MQ=255)                                                                                                aCAGGTTAGATGATTTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                  <  2:979556/151‑1 (MQ=255)                                                                                                     ttATATGATCTTTTGGTCCTTAAATGCC‑TT‑AAAAACATGCATGGGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                              <  2:47858/150‑1 (MQ=255)                                                                                                     ttAGTTGTTTTTTTTGTCGTTAAATCCC‑TG‑AAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                               <  2:823408/149‑1 (MQ=255)                                                                                                       aGATAATTTTTTTGCCGCTAAATCCT‑TA‑AAAAACACACATGACCCACAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                             <  2:556216/149‑1 (MQ=255)                                                                                                           gttCTTTTTGTCGCTTAATGCC‑TT‑AAAAACAGGCATGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                        <  2:498249/148‑1 (MQ=255)                                                                                                            atCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                      >  2:202307/1‑151 (MQ=255)                                                                                                              tttttttCCGCTTAATGCC‑TA‑AAAAACATTCAAGCCCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                      <  2:806416/148‑1 (MQ=255)                                                                                                                  ttGTCTTTTAATGCC‑TG‑AAAAACATGCAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                 <  2:396925/150‑1 (MQ=255)                                                                                                                  tgttCGCTTAATGCT‑TG‑AAAAACATGCATGACCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                 <  1:582489/147‑1 (MQ=255)                                                                                                                    gTCGTTAAATGCT‑TG‑AAAAACATGCATGACCCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                <  2:504305/149‑1 (MQ=255)                                                                                                                      cgtttAATGCC‑TA‑AAAAACTTGCATGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                            <  1:1015545/148‑1 (MQ=255)                                                                                                                             aatccttaaaaaaCATGAATGCCCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGaa                          <  1:229539/139‑1 (MQ=255)                                                                                                                               cc‑TG‑AAAAACATGAATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                   <  2:377719/151‑1 (MQ=255)                                                                                                                                c‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                  >  1:423396/1‑151 (MQ=255)                                                                                                                                   tatAAAACATGCCTGCCCCAAAAAATAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                     <  2:75400/145‑1 (MQ=255)                                                                                                                                   taaaaaaCTTGAATGCGCCAAAAAATAAAAAAAAATATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                   <  1:369325/146‑1 (MQ=255)                                                                                                                                     aaaaaCATGCATGCGCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                <  1:59474/149‑1 (MQ=255)                                                                                                                                       aaaCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt             <  2:880462/151‑1 (MQ=255)                                                                                                                                              catgcaccAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg       <  1:588397/144‑1 (MQ=255)                                                                                                                                              catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg       >  2:875017/1‑150 (MQ=255)                                                                                                                                               atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg       >  2:314440/1‑149 (MQ=255)                                                                                                                                                tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt    >  2:358644/1‑151 (MQ=255)                                                                                                                                                tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt    >  1:1053821/1‑151 (MQ=255)                                                                                                                                                    ccAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga  <  2:556686/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  TTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGA  >  NC_000913/29339‑29631

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