Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,583,969	G→A	W98* (TGG→TAG)	insB6 →	IS1 family protein InsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,583,969	0	G	A	97.1%	97.9 / ‑3.7	34	W98* (TGG→TAG)	insB6	IS1 family protein InsB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/4);  total (29/5)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGCCGGGCAGTGACGTCATCGTCTGCGCGGAAATGGACGAACAGTGGGGCTATGTCGGGGCTAAATCGCGCCAGCGCTGGCTGTTTTACGCGTATGACAGTCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTA  >  NC_000913/3583771‑3584262                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      aGCCGGGCAGTGACGTCATCGTCTGCGCGGAAATGGACGAACAGTGGGGCTATGTCGGGGCTAAATCGCGCCAGCGCTGGCTGTTTTACGCGTATGACAGTCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATACCGCCTGAAGGGAAAGCTGCACGGAAGCAGCAAGAGGTATGCGCAGTGGATGGAGCgggagggggggagtc                                                                                                                                                                                                  <  1:651428/300‑14 (MQ=1)                                                                                                                                         aCGCACTATGGCTACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAATTTGGAGTCATTACCGACATTAAAAg                                                          <  1:620258/300‑1 (MQ=18)                                                                                                                                         aCGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                          <  1:620266/300‑1 (MQ=18)                                                                                                                                           gCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTa                                                        >  1:9566/1‑300 (MQ=21)                                                                                                                                            cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                                       >  2:637092/1‑300 (MQ=21)                                                                                                                                            cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                                       >  1:569305/1‑300 (MQ=21)                                                                                                                                            cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATCAACGACATTAAAAGt                                                         >  2:380707/1‑298 (MQ=17)                                                                                                                                            cACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCAACCTTTGACGTGGTGATATAGATGACGGATGGCTGGCAGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGTGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATCAACGACATTAAAAGTaa                                                       >  2:380611/1‑300 (MQ=255)                                                                                                                                              cTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCTGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTAATTACCGCCATTAAAAGTAACt                                                     >  2:167382/1‑300 (MQ=19)                                                                                                                                              cTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATCATCTGAACATAAAACACTATCAATAAGTTGGAGTCGTTACCGACATTAAAAGTAACt                                                     >  2:193667/1‑300 (MQ=25)                                                                                                                                                      aCGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAAACTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaa                                             >  1:111974/1‑300 (MQ=255)                                                                                                                                                       cGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                                            >  1:16660/1‑300 (MQ=39)                                                                                                                                                           ggggCGTCTTATGAGCCTGCTGTCACCCTTTGACGTTGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTAACGACATTAAAAGTAACTTTTTTGAAATTa                                         >  1:583592/1‑299 (MQ=255)                                                                                                                                                           ggggCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCANAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                         >  2:514173/1‑299 (MQ=255)                                                                                                                                                               cGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat                                    >  2:579823/1‑300 (MQ=255)                                                                                                                                                                gTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata                                   >  1:167838/1‑300 (MQ=255)                                                                                                                                                                 tCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAANCGATATACGCAGCGAATTGAGCGGCATAACCTNAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCCTGAAAAAGTCATCGGGCATTATCTGAACATAACACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                                  >  2:551793/1‑300 (MQ=255)                                                                                                                                                                   ttATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGgtaata                                >  2:671599/1‑300 (MQ=255)                                                                                                                                                                   ttATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata                                >  1:711331/1‑300 (MQ=255)                                                                                                                                                                   ttATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATNATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata                                >  1:506797/1‑300 (MQ=255)                                                                                                                                                                   ttATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCAGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata                                >  2:290793/1‑300 (MQ=255)                                                                                                                                                                    tATGAGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat                               <  1:106355/300‑1 (MQ=255)                                                                                                                                                                        aGCCTGCTGTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                           >  2:608646/1‑300 (MQ=255)                                                                                                                                                                             gctgTCACCCTTTGACNTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACACAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAANTAACTTTTTTGAAATTAGNGATAATATTGGCACAGa                      >  1:611487/1‑300 (MQ=255)                                                                                                                                                                             gctgTCACCCTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGa                      >  2:477186/1‑300 (MQ=255)                                                                                                                                                                                    cccTTTGACGTGGTGATCTAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGCACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAACAGTAACTTTTTTGCAATTAGAGATAATATTGGCACAGAAAatata               >  1:22217/1‑300 (MQ=255)                                                                                                                                                                                    cccTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACANTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAatata               >  1:512417/1‑300 (MQ=255)                                                                                                                                                                                     ccTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAatatat              >  1:36325/1‑300 (MQ=255)                                                                                                                                                                                      cTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAACCACTATCAATAAGTTGGAGTCAGTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt             >  2:386185/1‑300 (MQ=255)                                                                                                                                                                                        ttGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGTGGCATAACCTGAATCTTAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTATAGATAATATTGGCACAGAAAATATATTc            <  2:343762/299‑1 (MQ=255)                                                                                                                                                                                                gtgATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTTAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGttt   >  1:50991/1‑300 (MQ=255)                                                                                                                                                                                                gtgATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGttt   >  2:4936/1‑300 (MQ=255)                                                                                                                                                                                                gtgATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGttt   >  1:538080/1‑300 (MQ=255)                                                                                                                                                                                                  gATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAATAACTTTTTTGAAATTAGAGATAATATTGGAACAGAAAATATATTCAGTTCGtcta  >  2:66/1‑296 (MQ=255)                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      AGCCGGGCAGTGACGTCATCGTCTGCGCGGAAATGGACGAACAGTGGGGCTATGTCGGGGCTAAATCGCGCCAGCGCTGGCTGTTTTACGCGTATGACAGTCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTA  >  NC_000913/3583771‑3584262

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