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

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29348‑29642                                                                                                                                                      |                                                                                                                                                   tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                   <  1:180989/150‑1 (MQ=255)    tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATATTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                                >  1:257207/1‑150 (MQ=255)     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               >  2:887772/1‑150 (MQ=255)          ttttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                         >  2:362452/1‑151 (MQ=255)            ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:500622/1‑151 (MQ=255)                     tGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGAAATCTTCTTGCTGCGCAAGCGTTTTTC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAATCATGCATGAGCCAGACATTAATATAAAATATCCCGCCATTAAGTTg                                                                                                                              >  1:233743/1‑151 (MQ=255)                      ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                              >  2:624375/1‑150 (MQ=255)                          cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                         >  2:348176/1‑151 (MQ=255)                             ttttttCTGCTTGTTTTTATTTAATGAATTTTTCTTGCTGCACAGGCTTTTCCA‑AAAAAAGGATAAATAATTTTTTTGTCATTAAATGCTTGTAAATAATAGTTGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  1:800141/150‑1 (MQ=255)                             ttttttCTGCTCTTTTTTTTTTCATGCAATCTTCTTGCTGGGCAGGTTTTTTCC‑AGAACAGTTTAGATTTTCTTTTTGTCGTTTATTCCCTAAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  2:630021/150‑1 (MQ=255)                              tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                     >  2:526332/1‑151 (MQ=255)                                    tGCTCTTTTTTTTTTCTTGCATTCTTCTTGCTGCGCAAGCTTTTCCC‑AGAACAGTTAAGATGTTTTTTTTGTCGTTAAATGCGTGAAAAACAGGAAGGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                                    <  2:290914/146‑1 (MQ=255)                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                           >  2:731889/1‑151 (MQ=255)                                         gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                          >  1:256512/1‑151 (MQ=255)                                          tttttATTCCATCCATTCTTCTTGCTGCGCAAGCTTTTTCA‑AGAACAGTTTAGATATTTTTTTTGTCGTTTAATGCCTGAAAAACATGCGTGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                           <  2:196054/149‑1 (MQ=255)                                           ttttttttCATGCATTCTTCTTGCTGCACATGTTTTTCCA‑AAAAAAGGTTATATGTTTTTTTTGTCGCTTAATGCCTGTAAAACATGAGTGACCCACAAAATAAAAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                          <  1:785767/149‑1 (MQ=255)                                             ttATTTCATGCATTCTTCTTGCTGCGCAAGCTTTTCCC‑AGAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                       <  2:895404/150‑1 (MQ=255)                                              tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                      >  1:120060/1‑150 (MQ=255)                                                tttCATACAATTTTCTTGCTGCGCAAGCTTTTCCC‑AAAACAGTTAAAATGATTTTTTTTTCTTTTAATGCCTTAAAAACATGCATGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                      <  1:933406/148‑1 (MQ=255)                                                       cAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                            >  2:450400/1‑151 (MQ=255)                                                       cAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                            >  1:634036/1‑151 (MQ=255)                                                        aaTCTTCTTGCTGCGAAAGCTTTTCCC‑AAAACGGTTTAGATGTTTTTTTTGCCGCTTAATGCTTAAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                            <  2:449227/150‑1 (MQ=255)                                                            tttttGCTGCCCAAG‑GTTTTCCAAAAACGGTTTAAATAATTTTTTTTTCGCTTAATCCCTTAAAAACATCCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                              <  1:845352/141‑1 (MQ=255)                                                            ttcttGTGGCGAATGCTTTTCCA‑AAAACGGTTTAGATTTTTTTTTTGTCGTTAAATCCTTAAAAAACATGAAGGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                           <  2:541019/147‑1 (MQ=255)                                                            ttcttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  2:59171/1‑151 (MQ=255)                                                            ttcttGCAGCGAAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                        >  2:862992/1‑150 (MQ=255)                                                              cttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  2:551450/1‑149 (MQ=255)                                                               ttGCTGCCGAAGCTTTTCCA‑AAAAAAGGTTAGATGTTTTTTTTTTCGTTTAATCCTTAAAAAACATGCAGGGCCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                          <  1:497315/145‑1 (MQ=255)                                                                             ttttCC‑AGAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTAAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                             <  2:519828/144‑1 (MQ=255)                                                                                  caaaaacattttAGTTAATTTTTTTGTCTCTTAATCCTTTTAAAACATGTTTGCACCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                   <  1:388501/140‑1 (MQ=255)                                                                                    aaaaCAGGTTAGATTTTTTTTTTTTCTTTTAATGCCTAAAAAACATGCATGACCCAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                <  2:852126/149‑1 (MQ=255)                                                                                      aaCGGGTTAGATGATTTTTTTGTCGTTTAATCCTTGAAAAACATGAATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAt                                                                 <  2:757502/148‑1 (MQ=255)                                                                                         aGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                           >  2:188809/1‑151 (MQ=255)                                                                                            ttAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                         >  2:243602/1‑150 (MQ=255)                                                                                                atgatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     >  1:584450/1‑150 (MQ=255)                                                                                                 ttttttttttGTCGCTAAATCCCATAAAACCATGTATGCCCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     <  1:265415/144‑1 (MQ=255)                                                                                                    tttttttGTCTCTTAATGCTTTAAAAAAATGCAGGACCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                   <  1:412056/146‑1 (MQ=255)                                                                                                    tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                >  2:628503/1‑151 (MQ=255)                                                                                                     ttttttGTCTTTTAATGCCTAAAAAACATGCATGCGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                <  1:767192/149‑1 (MQ=255)                                                                                                     cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                >  1:655137/1‑150 (MQ=255)                                                                                                      tttttGTCCTTTAATGCCTAAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                 <  1:264184/148‑1 (MQ=255)                                                                                                       ttttGTCGATAATTGCTTGAAAAACAGGCATGCACCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                 <  2:98570/147‑1 (MQ=255)                                                                                                           gTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                         >  2:410648/1‑151 (MQ=255)                                                                                                                ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                    >  1:359169/1‑151 (MQ=255)                                                                                                                ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                    >  2:727935/1‑151 (MQ=255)                                                                                                                  aaggcttgaaaaaCATGAAGGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                   <  2:58165/141‑1 (MQ=255)                                                                                                                         gTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  2:493044/1‑150 (MQ=255)                                                                                                                           aaaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                          >  2:660235/1‑150 (MQ=255)                                                                                                                            aaaCATGTATGAGCCAAAAAAAAATATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        <  1:345980/151‑1 (MQ=255)                                                                                                                            aaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAGTTAATATGCAAAt                        >  1:337260/1‑151 (MQ=255)                                                                                                                             aaCATGAATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                       <  2:339036/151‑1 (MQ=255)                                                                                                                                    atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  >  2:199221/1‑149 (MQ=255)                                                                                                                                    atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  >  2:806645/1‑149 (MQ=255)                                                                                                                                      gacccAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                    <  2:631313/142‑1 (MQ=255)                                                                                                                                           aaaaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           <  1:551417/147‑1 (MQ=255)                                                                                                                                           aaaaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           <  1:592111/147‑1 (MQ=255)                                                                                                                                            aaaaaTAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt        <  2:34438/150‑1 (MQ=255)                                                                                                                                             aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       >  2:753477/1‑151 (MQ=255)                                                                                                                                             aaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       <  2:792041/151‑1 (MQ=255)                                                                                                                                                ataataTAAAAAATCCCGCCATTAAGATGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    >  2:82246/1‑151 (MQ=255)                                                                                                                                                  aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg   <  2:86571/150‑1 (MQ=255)                                                                                                                                                   ataaAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  <  1:527768/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29348‑29642

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