Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,950,342	Δ1 bp	coding (21/99 nt)	ilvL →	ilvG operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,950,342	0	G	.	98.4%	272.7 / ‑3.4	64	coding (21/99 nt)	ilvL	ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (26/37);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAC  >  NC_000913/3950205‑3950479                                                                                                                                          |                                                                                                                                           ttATATTGAAATGATTAAACGCATCATATGAATCCGCCAACCAAGATCTTTTACTATTTACAAAACCTATGGTAACTCTTTAGGCATACCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGc                                                                                                                                      <  1:244840/143‑1 (MQ=255) ttATATTGAAAGTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGc                                                                                                                                      <  2:155206/143‑1 (MQ=255)    aatTGAAAGTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTg                                                                                                                                   <  1:56351/142‑1 (MQ=255)     atTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTgg                                                                                                                                  >  2:702723/1‑143 (MQ=255)       tGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTc                                                                                                                                <  1:361875/143‑1 (MQ=255)           attattAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACACGACGCAAGACCAGACGATATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGa                                                                                                                            <  1:1038327/143‑1 (MQ=255)              attaAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTCGCCTGGTCGTGATTa                                                                                                                         >  1:131869/1‑143 (MQ=255)              attaAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTa                                                                                                                         <  1:260095/143‑1 (MQ=255)                taAACGGATCAGAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACGCTTTAGACATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGc                                                                                                                       <  2:635504/143‑1 (MQ=255)                        tcatAAAAGTCGGCATAAAGATAGCTTGTGCTATTTACAAAGCGTATGGTAACTCTTTAGGCATTCCTGCGAGCAAGATGCAAGGAAAGACAAAATGACGGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCgtggtggt                                                                                                               <  2:496802/143‑1 (MQ=255)                           tAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAACCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGat                                                                                                            >  1:847499/1‑143 (MQ=255)                            aaaaaTCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTACCTCTTTAGGCATTCCTTCGAACAAGATGCAGGCGAAGACTAAATGTTAGTCCTGTTACGACT‑ATTACCCTGGTCGTGATTAGCGTGGTGGAGata                                                                                                           >  1:681392/1‑142 (MQ=255)                             aaaaTCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGatta                                                                                                          <  1:483257/143‑1 (MQ=255)                              aacttctccAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCTAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGattat                                                                                                         <  2:935381/136‑1 (MQ=255)                                 tCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACATGATGCAAGAAAAGACAAAATGACAGACCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATtattat                                                                                                      <  2:1043812/143‑1 (MQ=255)                                 tCGGCCAAAAAAGATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATtattat                                                                                                      <  2:573841/143‑1 (MQ=255)                                  cGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATc                                                                                                     <  1:712365/143‑1 (MQ=255)                                    gCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCCAGAAAAAACATGAAGACTGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATGAGccc                                                                                                   >  1:1084669/1‑143 (MQ=255)                                       aaaaaaTATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCAcc                                                                                                >  1:936473/1‑143 (MQ=255)                                        aaaaaTATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGCACCTCACTCAAAAAAATACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCg                                                                                               <  1:505208/143‑1 (MQ=255)                                             tatCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCgg                                                                                          <  1:253278/143‑1 (MQ=255)                                               tCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATTACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCgggg                                                                                        <  2:886698/143‑1 (MQ=255)                                                 ttGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCt                                                                                      >  1:160131/1‑143 (MQ=255)                                                 ttGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCt                                                                                      >  1:698482/1‑143 (MQ=255)                                                    tACTATTTACAAAACCTATGGTAACTCTTTAGGGGTTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCCACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGCTTATTATCCCACCGTGCGGGGCTGCa                                                                                   >  2:1012909/1‑143 (MQ=255)                                                    tACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCGCTTCTACGACT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGCGGCTGCa                                                                                   >  2:884793/1‑143 (MQ=255)                                                       tATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACtt                                                                                >  2:537862/1‑143 (MQ=255)                                                        aTTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTg                                                                               >  2:696693/1‑143 (MQ=255)                                                          ttACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑a                                                                            <  1:1063513/143‑1 (MQ=255)                                                          ttACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑a                                                                            <  2:698482/143‑1 (MQ=255)                                                          ttACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑a                                                                            <  2:160131/143‑1 (MQ=255)                                                             cAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGGGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGa                                                                         <  2:1006768/143‑1 (MQ=255)                                                             cAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGa                                                                         <  2:235185/143‑1 (MQ=255)                                                              aaaaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAg                                                                        >  1:827681/1‑143 (MQ=255)                                                              aaaaCCTATGGGAACGGTGGAGGCATCCCGGCGAACAAGATGCAAGAAAAGACAAAATGACAGCTCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGGGGTGATTATTATCCCATCGTGCGGGGCTGCACTTGG‑ACGAg                                                                        <  1:666618/143‑1 (MQ=255)                                                                aaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGAGGTGGTGCTTACTATCCCACCGTGCGGGGCCGCACTTGG‑ACGAGGa                                                                      >  2:482756/1‑143 (MQ=255)                                                                   cTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTAAGAGT‑ATTAGCCTGGTCGTAATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGaagg                                                                   >  2:634809/1‑141 (MQ=255)                                                                       ggTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCtt                                                               >  1:742916/1‑143 (MQ=255)                                                                       ggTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCtt                                                               <  2:742916/143‑1 (MQ=255)                                                                           aCTCTGTAGGCAGTCCTTCGAACACGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTATCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGC‑ACGAAGAAAGGCTTagag                                                           <  2:603710/143‑1 (MQ=255)                                                                            ctctGTAGGCACGCCTTAGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAgaga                                                          <  2:562319/143‑1 (MQ=255)                                                                                    gCATTCCCTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCt                                                  <  2:1084669/143‑1 (MQ=255)                                                                                       ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACCGCCCTTCTACGAGT‑ATTAGCATGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCCGGGATGCACGTGG‑CCGAGGAAAGGCTCAGCGATCCAGCCGTaa                                               >  2:468989/1‑143 (MQ=255)                                                                                        tCCTTCGAACAAGATGCAAGAAAAGACAGAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGAGATCCGCATGGCGGTGCTTATTATCACACCGTGCGCGGCTGCACTTGG‑GCGAGGACAGGCTCAGAGCTCAAGCATTAAc                                              >  2:886521/1‑143 (MQ=255)                                                                                           ttCGAACAAGATGCAAGAAAAGACAACATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGAGTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGaa                                           <  1:299611/143‑1 (MQ=255)                                                                                           ttCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGaa                                           >  2:815799/1‑143 (MQ=255)                                                                                                  aaGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                    >  1:607220/1‑143 (MQ=255)                                                                                                     tgccAAGAAGCGACGAAAGGACAGCCCTTCTACGCGT‑AGTAGCCTGGTCGGCATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAcccc                                 <  2:889729/140‑1 (MQ=255)                                                                                                     aTGCAAGAAAAGACAAAATGACATCCCTTCTACGAGT‑ATTAGCCGGGTCGTGATTAGCGTGGTGGTGATTATTAGCCCACCGTGCGGGGCCACCATAAA‑AATATTAAAGGTTTATAGATCAAGCCTTAACGAACTAAGAcccc                                 <  1:207976/143‑1 (MQ=255)                                                                                                      tacAAGAAAAGACAAAATGACGGCCCTGCTACGTGT‑ATTAGTCTGTTCGTGATTAGCGTTGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAATCCTTAACGAACTAAGAccccc                                <  2:405208/141‑1 (MQ=255)                                                                                                        cAAGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGc                              <  1:4004/143‑1 (MQ=255)                                                                                                         aaGAAAAGACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGACTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATGAAGCCTTAACGAACGAAGACCCCCGCa                             >  2:963686/1‑143 (MQ=255)                                                                                                            aacaGACAAAATGACCGCCCTTCCACGAGT‑AGTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCg                          <  1:124578/140‑1 (MQ=255)                                                                                                                gACAAAATGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAg                      <  2:151909/143‑1 (MQ=255)                                                                                                                   aaaaTGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCACGCCTTAACGAACTAACACCCCCGCACCGCCAGGTc                   >  1:885230/1‑143 (MQ=255)                                                                                                                   aaaaTGACAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTc                   <  2:59767/143‑1 (MQ=255)                                                                                                                         acggctcTTGTACGAGTG‑TTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCggggg             <  2:448900/137‑1 (MQ=255)                                                                                                                         aCAGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGC‑AACAGGGTGATGTTTGTGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCggggg             >  1:1082545/1‑143 (MQ=255)                                                                                                                           aGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGGCCCGAAAGGACCAGCGGtg           >  2:317063/1‑142 (MQ=255)                                                                                                                           aGCCCTTCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGtt           >  1:138279/1‑143 (MQ=255)                                                                                                                             cacgtCTAGGAGT‑ATTGGCCTGGTCGTGATGAGGGGGGTGGTGATTATTATCCCACCGTGGGGGGTTG‑GCTCGGTACAACGGTGGGCTTATCCATCAACCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGtttt         <  1:1088091/139‑1 (MQ=255)                                                                                                                                 tCTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGtttttttt     <  1:531164/143‑1 (MQ=255)                                                                                                                                  cTACGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTg    >  2:113381/1‑143 (MQ=255)                                                                                                                                    aCGAGT‑ATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc  <  1:1030856/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           TTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGG‑ACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAC  >  NC_000913/3950205‑3950479

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