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

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29348‑29646                                                                                                                                                        |                                                                                                                                                       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:660046/1‑150 (MQ=255)    tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                                    >  2:366063/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:587968/1‑149 (MQ=255)                             ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                            >  1:556116/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:19383/1‑151 (MQ=255)                                                        aaTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                               >  2:1075882/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)                                                                                                                                   ggcaggaccaaaaaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                        <  2:335849/138‑1 (MQ=255)                                                                                                                                    gaatcaccaataaaaTAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                        <  1:39052/138‑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)                                                                                                                                               aaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt            <  1:340113/150‑1 (MQ=255)                                                                                                                                               aaaaaaATAAAAAATATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg                  <  2:317048/144‑1 (MQ=255)                                                                                                                                               aaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc           <  2:554133/151‑1 (MQ=255)                                                                                                                                                aaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct          <  2:964389/151‑1 (MQ=255)                                                                                                                                                  ataataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCAATATCTCAAAAATGCCCCCGTTTCCCAGAAATTCGTGGGTAAGCAGATTTGCATTGATTTACGTCATCTTTGGGATTTAATATGCAATAAAGGTGAGTGattatt            >  1:468077/1‑147 (MQ=255)                                                                                                                                                    aataAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      <  1:407299/151‑1 (MQ=255)                                                                                                                                                     atatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  1:462335/1‑150 (MQ=255)                                                                                                                                                      tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  2:386306/1‑149 (MQ=255)                                                                                                                                                        taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  2:476637/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29348‑29646

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