Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,155,197	G→A	Q394Q (CAG→CAA)	mdtA →	multidrug efflux pump membrane fusion protein MdtA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,155,197	0	G	A	85.7%	92.7 / 4.6	42	Q394Q (CAG→CAA)	mdtA	multidrug efflux pump membrane fusion protein MdtA
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (21/15);  total (25/17)
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.69e-01

CATGTTAGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGGTGGTCACGCTCT  >  NC_000913/2154951‑2155440                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    cATGTTAGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAgcg                                                                                                                                                                                                                                                 <  1:73335/251‑1 (MQ=255)        ttctatacctatcataaCGCCTTAGTGATTCCAACAGCTGCCATGAAAAAGGGCAAGGAAGGCCAATTGGTCTGGGTGATGAATAGCGAAAAAAATGTCAGCAAACATCCGGGGACGCCGGGCACTCAGGACAGTCAGAAAGTTGGGAACCGTTCAGGTATTTATTCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAATAGTCGAAAGTTTAAGTGGTCGAATCCCAAAGCGCCactc                                                                                                                                                                                                                                           <  2:163744/235‑2 (MQ=255)            acacCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGCAAACAAGGTCAGCAAACATCTGGTGCCGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGATGGGGCGAAAGTGGATGTGGTGGAAGCCCAGAGCGCCACTACTCCg                                                                                                                                                                                                                                      >  2:115982/1‑251 (MQ=255)                               aGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGcc                                                                                                                                                                                                                   <  2:117557/251‑1 (MQ=255)                                   aTCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGa                                                                                                                                                                                                               <  1:204927/251‑1 (MQ=255)                                        aaCAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATAcg                                                                                                                                                                                                          <  1:266506/251‑1 (MQ=255)                                                   ggtaaaTTGGCAATGAAGGCCATTGTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGAGCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACGCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGa                                                                                                                                                                                               <  2:279870/248‑1 (MQ=255)                                                     cAAATGGGCAGTGAAGGCCATTTGTTCTGGGTGGTGAATATGGAAAACAAGGTCAGCAAACATCTGGTGACGCCGTGCAGTCAGGACAGTCAGAAAGTTGTGAGCCTTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACATACGGCATTGAGCGCCTTACCGAAGGGGCGAAAGTGGAAGTTGTGGAAGCCCAAAGCGCCACTACGCCGGAAGAGAAAGCCACCAGCCGCGAAGACGCGTGAAAAGGAg                                                                                                                                                                                              <  2:169180/250‑1 (MQ=255)                                                             cAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCt                                                                                                                                                                                     >  1:248836/1‑251 (MQ=255)                                                              aatgaatgCCATTTTTGCTGGGTGCTGAATATCGAAAACAAGGTCAGCAAACGTCGGGTGACGCCGGGAATTCAGGACAGGCAGAAAGTGGAGATCCGTGCAGGTATTTCGGCGGGCGATCGCGTGGTGACAGACGGCAGTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAGGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTg                                                                                                                                                                                    <  2:181816/251‑1 (MQ=255)                                                               aTGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGAAAGCTCCTgg                                                                                                                                                                                   >  1:78958/1‑250 (MQ=255)                                                                      ccATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGACAGAACGTGGTGATCCGAGCCGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGGTCGCCTGACCGCAGGGGCGAAAGGGGAAGTGGTGGAAGCACAGGGCGCCACAACTCCGGAAGATAAAGACACCAGCCTCGCACAAGCACAAAAAGGAGCAAGCTCCGGAAGCAggg                                                                                                                                                                            >  2:95484/1‑250 (MQ=255)                                                                                   gTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGagc                                                                                                                                                               <  1:206199/251‑1 (MQ=255)                                                                                      cTGAATAGCGACAACAAGGTCAGCAAACATCTGGAGCCGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGGGGTGCCAGACGGCATTGATCGCCTGACCGAAGGGGCGCAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAACCACCAGCCGCGAATACGCGCAAAAAGGAGCCAGCTTCTGATGCCAGTGTTAACCCCGggcatc                                                                                                                                                            >  2:203142/1‑249 (MQ=255)                                                                                       tGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGagcagc                                                                                                                                                            >  1:157440/1‑250 (MQ=255)                                                                                          aTAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGcac                                                                                                                                                          >  1:38523/1‑249 (MQ=255)                                                                                              cGAAAACAAGGTCAGCAAACATCTTGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGCTCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCTAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcgg                                                                                                                                                    <  1:281929/251‑1 (MQ=255)                                                                                               gAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGAGTTACCCCCGAGCAGCACAggcggc                                                                                                                                                   >  1:197485/1‑251 (MQ=255)                                                                                                aaaaCAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCCCCACTACTCCGGAAGAGAAAGCCACCAGACGCGAATACGCGAAAAAAAGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGcc                                                                                                                                                  >  1:276341/1‑251 (MQ=255)                                                                                                    cAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTCTTTCTGCGGGCGATCGCGGGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGCGGGGGAAGCCCAAAGCGCCACTACTACGGAAGAGAAAGCCACCAGCCGCGAATACTCGAACACAGGAGCACACTCCTGCTGCCGGGGTTACCCCCGAGCAGCACCGGCGGCCCggc                                                                                                                                              >  2:67448/1‑249 (MQ=255)                                                                                                        gTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCCGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAACAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGCGCAGCACAGGCGGcccagcgcgc                                                                                                                                          >  2:107400/1‑244 (MQ=255)                                                                                                        gTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATTCGTGCAGGTATTTCTGCGGGCGATCGCGTGGGGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGAGGCGGTGGGTGACGCCCAAAGCGCCACTACTCCGGATGGGAACGCCACCAGCCGCGAATACGCGAACAACGGGGCAGGCTCCTGAAGCAGGAGGTGCACCCGAGCAGCCCAGGCGGCCCGCCgcgc                                                                                                                                          >  2:9775/1‑251 (MQ=255)                                                                                                         tCAGCAAACATCTGGTGACGCCGGGCATTCCGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCCCGCTCCTGATGCAGGTGTTACCCCCGCGCAGCACAGGCGGCACGTCGCGcc                                                                                                                                         >  1:143504/1‑251 (MQ=255)                                                                                                                      ggTGACGCCGGGAATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatta                                                                                                                                <  1:153664/247‑1 (MQ=255)                                                                                                                               gggCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGCTGCAGGTGTTACCCCCGCGCAGCCCAGGCGGCCCGTCGCGCCTATTTATTGTGCGTCCTGTGGc                                                                                                                   >  2:92808/1‑251 (MQ=255)                                                                                                                                      cAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACg                                                                                                            <  1:230060/251‑1 (MQ=255)                                                                                                                                       aGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGcc                                                                                                                                                                                                                   >  1:124205/1‑147 (MQ=255)                                                                                                                                       aGGACAGTCAGAAAGTGGTGATCCGTGAAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCCCCAGcc                                                                                                                                                                                                                   <  2:124205/147‑1 (MQ=255)                                                                                                                                                   aaGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGAtggtgg                                                                                               <  1:65421/251‑1 (MQ=255)                                                                                                                                                             acGTGTAGGTATTTCTGCTGGCGTTCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCTCTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTAc                                                                                      <  2:161188/249‑1 (MQ=255)                                                                                                                                                              cGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGCGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCCCGCTGCTGATGGGGGCGAACTTaatc                                                                                    >  2:164797/1‑248 (MQ=255)                                                                                                                                                                    ggTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCAACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGACACCACGCTGCTGATGGTGGCGATCTTACTCGCCggg                                                                              >  1:194642/1‑251 (MQ=255)                                                                                                                                                                          tCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGACAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATc                                                                        >  1:109939/1‑251 (MQ=255)                                                                                                                                                                                   gATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATcgc                                                               <  1:209654/251‑1 (MQ=255)                                                                                                                                                                                   gATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATcgc                                                               >  1:4788/1‑251 (MQ=255)                                                                                                                                                                                   gATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCAAGCTGCTGATGGTGGCGATCTTCCTCGCCGGGATTATCGGTTATcgc                                                               >  1:191569/1‑251 (MQ=255)                                                                                                                                                                                                 aGACGGCATTGATCGCCTGCCCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCACCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTCTTATGCGTCCTGTGGCACCCACGCGGCTGATGGTGGCGATATTACTCGCCGGGAGTATCGGTTCTCGCACCCTGCCCCTTTc                                                 >  2:88311/1‑251 (MQ=255)                                                                                                                                                                                                    cGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGCTCTTACTCGCCGGGATTATCGGTTATAGCGCCCTGCCCGTTTCGgc                                              >  1:271174/1‑251 (MQ=255)                                                                                                                                                                                                     ggCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGgcg                                             >  1:247853/1‑251 (MQ=255)                                                                                                                                                                                                     ggCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGGGGAAGCCTAAAGCGCCACAACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGgcg                                             <  2:58848/251‑1 (MQ=255)                                                                                                                                                                                                                        aaGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACAAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTGCTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGAGGACTATc                          >  2:272775/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGGGGACTATCCGACCATTCAGGTGGTCACGctct  >  1:90072/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CATGTTAGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGGTGGTCACGCTCT  >  NC_000913/2154951‑2155440

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