Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,446,482	C→T	L791L (CTG→TTG)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,446,482	0	C	T	97.6%	124.8 / ‑5.3	41	L791L (CTG→TTG)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/23);  total (17/24)
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.70e-01

AACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCAACTTGCGTCTG  >  NC_000913/4446246‑4446727                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      aaCTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGttt                                                                                                                                                                                                                                         <  1:139113/251‑1 (MQ=255)    tgtgCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTc                                                                                                                                                                                                                                        >  1:206638/1‑249 (MQ=255)               cAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCt                                                                                                                                                                                                                           >  2:83182/1‑251 (MQ=255)                               gggCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGt                                                                                                                                                                                                           <  2:186793/251‑1 (MQ=255)                                       gggCGTGCGGCGGTGAATCTCAACCGCTTCGACCCCGCCATGCCGAAACCGTTTATGCCAGAACCCACTCAGGCCCGCCGTATCTTCACCGGTACAGCAGCCGTTGCCTGGGACACCACGAAAGAGGCGCTGCCGCAGGGCAGTATCACCCTTTCGGGCCGTAACGTGCAGGTAACGCAAACCGTCAACGCTGCGGAGTCGCCCGTGGCGTCTCAGACACTCACTATGCCGGCGGAATTGCGTAAAAACCg                                                                                                                                                                                                   >  2:179829/1‑251 (MQ=255)                                          cGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTCTGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACAGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTCACAACCGTGc                                                                                                                                                                                                >  1:165670/1‑251 (MQ=255)                                                                           gCCATTCTGAAACCGTTTATGCCAGAAACCACTCAGGCCATCGGTATCTTCAAGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGTGCAGTATCACCCTTTCGGGGCGTAACGTGCATGTAACGCTACCCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGTATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTTGGCTGGACCATCCGCCTGACCAATa                                                                                                                                                                <  1:149666/250‑1 (MQ=255)                                                                           gCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCGGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATaa                                                                                                                                                               <  2:103513/251‑1 (MQ=255)                                                                             cATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACg                                                                                                                                                             <  1:130360/251‑1 (MQ=255)                                                                                       ccggttATGCCAGAAACCACTCAGGCCAGCGGTATCTGCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTGTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACGATGCGCCTGACCAATAACGGCGAGTTTGa                                                                                                                                                   <  2:37604/247‑1 (MQ=255)                                                                                            tATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGa                                                                                                                                               <  1:105560/250‑1 (MQ=255)                                                                                             aTGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAc                                                                                                                                              <  1:137174/250‑1 (MQ=255)                                                                                                       cgacgcAGGCCAGCGGGATCTTGACGGGTAAAGCGGTCGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCCCCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTGCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCGGGACCATCCGCCTGACCACTCATGGCCAGTCTGATGGAcaggtgcaggtg                                                                                                                                   <  2:110368/246‑1 (MQ=255)                                                                                                         aCTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCCGTTTGATGGACAGGTGCAGGTGAc                                                                                                                                 >  1:268620/1‑251 (MQ=255)                                                                                                              ggCCAGCGGTATCTTTACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGTTGCAGGTGACCGATc                                                                                                                            <  2:40542/251‑1 (MQ=255)                                                                                                                  aGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTTCAGGTGACCGATCCGCa                                                                                                                        <  1:84699/251‑1 (MQ=255)                                                                                                                  aGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGACGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATACCGGCCAGTTTGATGGGCAGGTGCAGGTGACCGATCCGcc                                                                                                                        >  2:108177/1‑250 (MQ=255)                                                                                                                        aTCTTCACGGGTAAAGCGGACGTTGCCTGGGACCCCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgcc                                                                                                                   <  2:147606/250‑1 (MQ=255)                                                                                                                          cTTCACGGGTCAAGCTGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgccgcc                                                                                                                >  1:206712/1‑251 (MQ=255)                                                                                                                                      aGCGGACGTGGCCTGGAACCCCACGAAAGAGGGGCGGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAAGGTGCAGGTAACGCAAACCGTCAACGATGGGGCGTTGCCCGTGGCGTTTCAGACCCGGAATCTGACGGCGGAATGGCGTAACACCCGTGCCGAATTGGGCTGGACCGTCGGCCTGACCAATAACGGCCAGTTGTATGGACAGGTGCAGGTGACCGAGCCGTGAGTCCGCCGTAATCTtggtg                                                                                                    <  2:134799/251‑1 (MQ=255)                                                                                                                                           aCGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAAc                                                                                               >  2:59712/1‑251 (MQ=255)                                                                                                                                           aCGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAAc                                                                                               <  1:120384/251‑1 (MQ=255)                                                                                                                                                gCCTGGGACACCCCGAAAGAGGGGTTGCCGCAGGGCAGTATCCCCTTTTCGGGGCGTAACGGGCAGGTAACGCAAACCGTAAACGATGCGGCGTTGCCCGGGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGGGCCGAATTGGGCGGGACCATCCGCCTGACCAATAACGGCCAGTTGGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCa                                                                                           <  2:222394/250‑1 (MQ=255)                                                                                                                                                   tGGGACACCAAGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAatat                                                                                       >  1:191847/1‑251 (MQ=255)                                                                                                                                                    gggACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCCCCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATc                                                                                      <  1:194602/251‑1 (MQ=255)                                                                                                                                                    gggACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATc                                                                                      >  1:19453/1‑251 (MQ=255)                                                                                                                                                     ggACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATcc                                                                                     >  1:159268/1‑251 (MQ=255)                                                                                                                                                      gACACAACGAAAGAGGGGCTGCCGCAGGCCAGAATGACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAAAGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGAGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGACTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCACGGCCGCCGCAATCTTGGTGGCAACGTCTATATCCg                                                                                    >  2:13860/1‑251 (MQ=255)                                                                                                                                                            acGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACtt                                                                              >  1:57939/1‑251 (MQ=255)                                                                                                                                                             cGAAAGCGGGGCTGCCGCAGGGCAGTATCACCCTTGCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTc                                                                             <  2:8533/251‑1 (MQ=255)                                                                                                                                                             cGAAAGAGGGGCTGCCGCAGGGCAGTCTCACCCTTTGGGGGCGAAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATGGGGCTGGACCACCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGGGACCGACCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTc                                                                             <  2:57939/251‑1 (MQ=255)                                                                                                                                                                      ggCTGCCGCAGGGCAGTATCCCCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCg                                                                    <  1:82467/251‑1 (MQ=255)                                                                                                                                                                      ggCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAACCCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATGGCGTAACAACCGTGCCGAATTGGGCTGGACCACCCGCCTGCCCAATAACGGCGAGTTTGAGGGCCATTTGCAGGTGACCGATCCGCAAGGCCGCCGTAAGCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCg                                                                    <  1:174944/251‑1 (MQ=255)                                                                                                                                                                      ggCTGACGCAGGGCAGTATAACACTTTCGGGGCGTAACGTGCCGGTCACGCAAACCGTCAACGATGCGGCGTTGCCCGCGGCGTTTCAGACACTGCATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGAAAACTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCg                                                                    >  2:176328/1‑251 (MQ=255)                                                                                                                                                                               aGGGCAGTATCACCATTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAAc                                                           >  2:167344/1‑251 (MQ=255)                                                                                                                                                                                         cACCCTTGCGGGGCGAAAGGTGCAGGTAACGCAAGACGTCTACGGTGCGGCGCTGCCCGTTGCCTTTCAGCCACTGAATCGGACGGCGGAATTGCGGAACAACCGGGCCGAATTGGGCTGGCCCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCTTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTa                                                 <  2:206653/251‑1 (MQ=255)                                                                                                                                                                                                                       gCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGCCACTGAATCTGACGGCGGAATTGCGTAACCACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACAGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTg                    >  2:127442/1‑250 (MQ=255)                                                                                                                                                                                                                            gcgtgaacggtggggcgTTGCCCGTGGCGTTTCCGACGCTGAATCTGACGGCGGAATTGCATAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAAGATCGGTCAGTTTGAGGGACAGGGGCAGGGGACCGAGCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAAGATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGcc              <  1:179829/238‑1 (MQ=255)                                                                                                                                                                                                                             cGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGGGAGTGcc              >  2:62452/1‑250 (MQ=255)                                                                                                                                                                                                                             cGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACGTAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATGGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCa             <  1:153655/251‑1 (MQ=255)                                                                                                                                                                                                                                        cggcgTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCAACTTGCGTCTg  <  1:83182/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      AACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCAACTTGCGTCTG  >  NC_000913/4446246‑4446727

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