Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,799,370	Δ1 bp	7.8%	intergenic (‑100/+23)	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,370	0	A	.	7.8%	179.9 / 10.7	64	intergenic (‑100/+23)	pheM/rplT	pheST‑ihfA operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base . (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-02

GGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGAT  >  NC_000913/1799229‑1799512                                                                                                                                              |                                                                                                                                                gttGCAAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGAc                                                                                                                                                <  2:326883/141‑1 (MQ=255)   tGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTGTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTa                                                                                                                                              >  1:162492/1‑143 (MQ=255)       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)                                                                                                                                              |                                                                                                                                                GGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGAT  >  NC_000913/1799229‑1799512

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