Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,799,374:1	+T	7.7%	intergenic (‑104/+19)	pheM ← / ← rplT	pheST‑ihfA operon leader peptide/50S ribosomal subunit protein L20

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,799,374	1	.	T	7.7%	242.0 / 15.2	65	intergenic (‑104/+19)	pheM/rplT	pheST‑ihfA operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base . (30/30);  new base T (2/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

AAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAT  >  NC_000913/1799235‑1799516                                                                                                                                             |                                                                                                                                               aaaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTc                                                                                                                                             <  2:638840/143‑1 (MQ=255)  gaGTAAAACAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTcc                                                                                                                                            <  2:1107801/142‑1 (MQ=255)    gTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCct                                                                                                                                          <  1:642674/143‑1 (MQ=255)     tAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctc                                                                                                                                         >  1:675394/1‑143 (MQ=255)       aaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCtt                                                                                                                                       >  2:274757/1‑143 (MQ=255)        aaacaGAAGCGGAAACTAGCAGCATTCATTGCTCGCGTTACCTTTTGGTACTCTCCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCATTGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttt                                                                                                                                      <  1:1057694/143‑1 (MQ=255)        aaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttt                                                                                                                                      <  1:409595/143‑1 (MQ=255)          aaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACACTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCa                                                                                                                                    >  1:901054/1‑143 (MQ=255)             gaagCGGAAAATAGCAGCATACATCGCTTGCGGTACCTTTAGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttttcgtt                                                                                                                                 <  1:151192/143‑6 (MQ=255)              aagCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTg                                                                                                                                <  2:126490/143‑1 (MQ=255)                gCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGc                                                                                                                              <  2:34642/143‑1 (MQ=255)                    acaagagcagcAGTCATGGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTAt                                                                                                                          <  1:394187/138‑1 (MQ=255)                     aaaTAGCAGCATTCAGTGCTTGCGTTACCTTTTAGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATg                                                                                                                         <  1:839132/143‑1 (MQ=255)                           cagcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATCGCGTTGAAAAGAAAAGAGGGCGACTAG‑CTCTGTGTTTCGCCTGGCTTACGccgaag                                                                                                                   >  2:680555/1‑139 (MQ=255)                              cattcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAgctg                                                                                                                >  2:394639/1‑143 (MQ=255)                                   attGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTc                                                                                                           <  1:630852/143‑1 (MQ=255)                                     tgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGc                                                                                                         <  1:407260/143‑1 (MQ=255)                                      gcgtgggTTGCCTTTTGGGACTCTTCAAAGGACCGTTCTTGAAAAGGCTACGGCGATAAAGGTCAATGTTTTGATGGCGTTGAAACGAGAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCt                                                                                                        <  2:27812/137‑1 (MQ=255)                                       cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                       >  2:524446/1‑143 (MQ=255)                                       cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAAGATTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                       >  1:607385/1‑143 (MQ=255)                                        ttgcGTTACCTTTTGGTACTCTTCAAACGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      <  1:419269/143‑1 (MQ=255)                                        ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAACAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      >  1:723673/1‑143 (MQ=255)                                        ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      >  2:307848/1‑143 (MQ=255)                                                 cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCCTGTTCAACCag                                                                                             >  2:8012/1‑143 (MQ=255)                                                  ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCaga                                                                                            >  1:133707/1‑143 (MQ=255)                                                   tttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCagag                                                                                           >  1:822059/1‑143 (MQ=255)                                                    ttGGTACTATTCAAAAGATCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGc                                                                                          <  1:69401/143‑1 (MQ=255)                                                          ctctTCAAAAGACCTTTATGGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCCTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGCGCGGTGat                                                                                    <  2:1091478/143‑2 (MQ=255)                                                             ttCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGCACGc                                                                                 >  2:1006768/1‑143 (MQ=255)                                                                aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTAc                                                                              >  2:454938/1‑143 (MQ=255)                                                                    gACCTGTATTGAAAGGGCTACGGCGATAACAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTGCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTg                                                                          <  2:102735/143‑1 (MQ=255)                                                                      tcttgagtGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                        <  2:480996/136‑1 (MQ=255)                                                                      cgtttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                        <  2:539658/141‑1 (MQ=255)                                                                      ccTTTATCGAAAAGGCCACGGCGATAAAAGTCAATGGTTTGATGGCGTTGAAACGAAAAGAGGGAGACCAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                        <  1:208374/143‑1 (MQ=255)                                                                       cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGACAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTtcgtcg                                                                       >  2:75943/1‑143 (MQ=255)                                                                        tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGa                                                                      >  1:794701/1‑143 (MQ=255)                                                                         ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGaa                                                                     >  1:255666/1‑143 (MQ=255)                                                                          tATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGaac                                                                    >  1:795230/1‑142 (MQ=255)                                                                              gAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                                <  1:582259/143‑1 (MQ=255)                                                                               aaaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                               <  1:464451/143‑1 (MQ=255)                                                                                 gaGGCTACGGCGATAAAAGTCAAGGTTTTGATGTCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGa                                                             >  2:1056548/2‑143 (MQ=255)                                                                                         ggCGATAAAAGTCAATGTGTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACAGGCTTATGACAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGGTACTTTGTCGAATACTGCGATATCAGcc                                                     >  2:671452/1‑143 (MQ=255)                                                                                            gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAgg                                                  <  2:616088/143‑1 (MQ=255)                                                                                             aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTCCTTTGTCGAATACTGCGATATCAGCCAGGa                                                 >  2:872810/1‑143 (MQ=255)                                                                                             aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGCATACTGCGCTCTCAGCCAGGa                                                 >  2:809515/1‑143 (MQ=255)                                                                                                  aGTCAATGTTTTGATGGCGCTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATCCTGCGATATCAGCCAGGATCTTa                                            >  1:236087/1‑143 (MQ=255)                                                                                                             tGATGGCGTTGAGACGAAAAGAGGGAGACTAG‑CTCCCCCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTc                                 <  2:920431/143‑1 (MQ=255)                                                                                                             tGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTc                                 <  1:329684/143‑1 (MQ=255)                                                                                                               aTGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCaa                               <  1:871694/143‑1 (MQ=255)                                                                                                                tGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAc                              >  2:246842/1‑143 (MQ=255)                                                                                                                gggCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAc                              <  2:1052333/142‑1 (MQ=255)                                                                                                                 ggCGTTGAAACGAAAAGAGGGAGACTCG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTAGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCCATATCAGCCAGCATCTTACGGGCGCGTTCAACa                             >  2:656163/1‑143 (MQ=255)                                                                                                                   cGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                           <  1:353870/143‑1 (MQ=255)                                                                                                                     ttGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTGCGGTCGATTGCAACAGAgg                         >  1:741278/1‑143 (MQ=255)                                                                                                                         aaCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGGTTTCAACAGAGGCttt                     >  1:375073/1‑143 (MQ=255)                                                                                                                             aaaaGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGCGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTc                 >  2:165806/1‑143 (MQ=255)                                                                                                                             aaaaGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTAAGCCGGAGCTGCTTTCGCTATATCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGAGCTTACGGGCGGTCTAAAAAGAGGATTTTTTc                 >  2:21388/1‑143 (MQ=255)                                                                                                                              aaaGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCa                <  1:274757/143‑1 (MQ=255)                                                                                                                               agggggggGACTAG‑CTCCCTCGTTCAACTGGCTTATGCCAGAGCTGCTTTTGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCCGCCAGGATCTTACGGTCGATTTCAACGGAGTCTTTTTTCAg               <  2:82543/139‑1 (MQ=255)                                                                                                                               aaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg               >  1:489540/1‑143 (MQ=255)                                                                                                                                agagGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAgg              >  2:1100427/1‑143 (MQ=255)                                                                                                                                    ggAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAAAGAGGCTTTTTTCAGGCCAt          >  1:458775/1‑143 (MQ=255)                                                                                                                                     gagaCTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGCGCTGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt         <  1:630005/143‑1 (MQ=255)                                                                                                                                      agaCTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg        <  1:762639/143‑1 (MQ=255)                                                                                                                                        aCTAG‑CTCCCTCTTTCAACTGGCTTATGCGAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgat      <  2:566556/143‑1 (MQ=255)                                                                                                                                            g‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGCTCCTACGGTCGATTCCAACCGAGGCTTTTTTCAGGCTATTGCTGaaa  >  2:812065/1‑142 (MQ=255)                                                                                                                                             |                                                                                                                                               AAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAT  >  NC_000913/1799235‑1799516

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