Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	895,016	A→G	6.3%	Q9R (CAG→CGG)	potG →	putrescine ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	895,016	0	A	G	6.3%	197.2 / 2.2	63	Q9R (CAG→CGG)	potG	putrescine ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base G (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCTGTGGCAAG  >  NC_000913/894869‑895164                                                                                                                                                    |                                                                                                                                                     gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAg                                                                                                                                                     <  4:48107/149‑1 (MQ=255) gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAg                                                                                                                                                     <  5:160504/149‑1 (MQ=255)  tGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAgg                                                                                                                                                    <  6:365796/149‑1 (MQ=255)    gACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGc                                                                                                                                                   <  8:43376/148‑1 (MQ=255)     aCCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCTGCGGGc                                                                                                                                                   >  5:263954/1‑147 (MQ=255)     aCCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCTGCGGGc                                                                                                                                                   >  5:263956/1‑147 (MQ=255)         aaGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAAc                                                                                                                                             >  1:297753/1‑149 (MQ=255)           gTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCTGCGGGCGagggata                                                                                                                                           <  6:263954/149‑7 (MQ=255)           gTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCTGCGGGCGagggata                                                                                                                                           <  6:263956/149‑7 (MQ=255)            tGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATCCTGCGCTAGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCACGCCAGCATGCGAAAACCCg                                                                                                                                          >  6:361112/1‑149 (MQ=255)              aaGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTa                                                                                                                                        <  8:261228/149‑1 (MQ=255)                     gAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCt                                                                                                                                 <  2:297753/149‑1 (MQ=255)                        aataatCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAc                                                                                                                              <  5:312311/149‑1 (MQ=255)                                   gTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTAt                                                                                                                      >  3:286339/1‑146 (MQ=255)                                      gTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                                                                                    <  2:175231/145‑1 (MQ=255)                                      gTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                                                                                    >  1:175231/1‑145 (MQ=255)                                        aGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAg                                                                                                                   <  6:362784/144‑1 (MQ=255)                                        aGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAg                                                                                                                   >  5:362784/1‑144 (MQ=255)                                         gATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGATAACCCGTAAGGCGCTGACGCCGCTATTAGAAATc                                                                                                              <  1:222766/148‑1 (MQ=255)                                           tGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGc                                                                                                           <  3:89056/149‑1 (MQ=255)                                              cGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAAc                                                                                                        <  1:24144/149‑1 (MQ=255)                                                   gggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGAc                                                                                                   <  8:226378/149‑1 (MQ=255)                                                     gcgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGAcc                                                                                                  >  8:80906/1‑148 (MQ=255)                                                       gcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaaa                                                                                               <  7:363089/149‑1 (MQ=255)                                                                ccgccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGAt                                                                                      >  5:179440/1‑149 (MQ=255)                                                                 cgccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATg                                                                                     <  1:201214/149‑1 (MQ=255)                                                                   ccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCACATCCTACGATGGt                                                                                   >  1:46227/1‑149 (MQ=255)                                                                   ccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCCACGATGGt                                                                                   <  7:80906/149‑1 (MQ=255)                                                                                 cATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCTGTGga                                                                     <  8:209511/149‑1 (MQ=255)                                                                                 cATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga                                                                     >  8:156463/1‑149 (MQ=255)                                                                                  aTTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACACGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat                                                                    <  5:361112/149‑1 (MQ=255)                                                                                   ttATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat                                                                    >  4:100258/1‑148 (MQ=255)                                                                                      tGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgatgat                                                                 <  5:311729/148‑1 (MQ=255)                                                                                        gTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga                                                                     <  4:2172/142‑1 (MQ=255)                                                                                        gTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga                                                                     >  3:2172/1‑142 (MQ=255)                                                                                        gTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGt                                                               >  2:70127/1‑148 (MQ=255)                                                                                         tGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGt                                                               >  7:203820/1‑147 (MQ=255)                                                                                           cgcTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGc                                                           >  3:32074/1‑149 (MQ=255)                                                                                            gcTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGcc                                                          <  1:70127/149‑1 (MQ=255)                                                                                              ttGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACACTATCCCTAGCCCTCAGGCGAAAACCCGTAAGGCGCTGAAGCCCCTATTAGAAATCCGCAACCTGAACAAATCCTACGAAGGTCAACATGATGTGGATGATCTCAGCCt                                                         >  3:109117/1‑148 (MQ=255)                                                                                                    acaTTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGa                                                                                                    <  3:162357/99‑1 (MQ=255)                                                                                                    acaTTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGa                                                                                                    >  4:162357/1‑99 (MQ=255)                                                                                                      aTTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGCCCATCTa                                                >  6:312603/1‑149 (MQ=255)                                                                                                      aTTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTa                                                >  1:169490/1‑149 (MQ=255)                                                                                                       ttCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTa                                                >  7:312946/1‑148 (MQ=255)                                                                                                               cGGAGAGCAGCCGTGAATGGCGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTg                                        >  2:268523/1‑148 (MQ=255)                                                                                                                 gagaGCAGCCGTGAATGGCGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGaaa                                     <  1:268523/149‑1 (MQ=255)                                                                                                                  agagCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAAt                                    <  7:156463/149‑1 (MQ=255)                                                                                                                      cagcCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTc                                >  2:127845/1‑149 (MQ=255)                                                                                                                             gAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGAAGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCgctgc                         >  2:373078/1‑149 (MQ=255)                                                                                                                                tGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTgg                      <  6:179440/149‑1 (MQ=255)                                                                                                                                tGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTgg                      >  4:317887/1‑149 (MQ=255)                                                                                                                                tGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTgg                      >  7:69044/1‑149 (MQ=255)                                                                                                                                     cTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCCGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCa                 >  5:220417/1‑149 (MQ=255)                                                                                                                                     cTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCa                 >  5:368995/1‑149 (MQ=255)                                                                                                                                        tCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATcc              >  8:44204/1‑149 (MQ=255)                                                                                                                                          ccTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCgg            <  8:312946/149‑1 (MQ=255)                                                                                                                                           cTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTg                                        <  3:10057/120‑1 (MQ=255)                                                                                                                                           cTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTg                                        >  4:10057/1‑120 (MQ=255)                                                                                                                                            tCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCt          >  5:14400/1‑149 (MQ=255)                                                                                                                                                 cGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCTGTGGc     <  8:203820/149‑1 (MQ=255)                                                                                                                                                   cAGGCGAAAACCCTTAAGTCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCATGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCTGTGGCaa   <  4:109117/149‑1 (MQ=255)                                                                                                                                                   cAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCTGTGGCaa   <  4:286339/149‑1 (MQ=255)                                                                                                                                                    aGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCTGTGGCAAg  >  5:314337/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAAAGGTGAAATCTTCGCGCTGCTGGGCGCATCCGGCTGTGGCAAG  >  NC_000913/894869‑895164

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