Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	582,597	C→T	intergenic (+500/‑278)	tfaD → / → tfaX	putative tail fiber assembly protein TfaD/protein TfaX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	582,597	0	C	T	97.6%	117.4 / ‑4.3	41	intergenic (+500/‑278)	tfaD/tfaX	putative tail fiber assembly protein TfaD/protein TfaX
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/22);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GAGTTCCATCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAA  >  NC_000913/582358‑582843                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       gagTTCCATCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGt                                                                                                                                                                                                                                             >  1:64950/1‑251 (MQ=255)              aGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGt                                                                                                                                                                                                                                >  1:294263/1‑251 (MQ=255)                    ggCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAAcc                                                                                                                                                                                                                          <  1:190848/251‑1 (MQ=255)                          cgccgcATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTc                                                                                                                                                                                                                    >  1:256474/1‑251 (MQ=255)                            caggaTTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCACATCCATTAACGATATCAGATTACTCCCAGGTACGTGAGTAAGGTAAAACCCAGTTcc                                                                                                                                                                                                                   >  2:271540/5‑250 (MQ=255)                               cagtGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAgg                                                                                                                                                                                                                                              <  2:207353/217‑1 (MQ=255)                               cATTGCGTATTACACAGCTTTTTTTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAgg                                                                                                                                                                                                                                              >  1:207353/1‑220 (MQ=255)                                       aTTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATc                                                                                                                                                                                                       >  1:396056/1‑251 (MQ=255)                                               gCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGgt                                                                                                                                                                                               <  2:359984/251‑1 (MQ=255)                                                            atatTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATg                                                                                                                                                                                  <  1:40432/251‑1 (MQ=255)                                                                   caggccggCAGTGAAGAAGGTAAGACATGGACATGGAATCAAATTGACCATGTCACGCCGCGGGATAAGGGTCAAAAACATCATGGCTACTTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTACAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACCCCAATATCCAGATGGTTGTTACCTAAATGTTCCaa                                                                                                                                                                            >  2:113595/1‑249 (MQ=255)                                                                             gTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCAGGGCTAATGTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGgaa                                                                                                                                                                 <  2:154638/251‑1 (MQ=255)                                                                                              tggaaatggaaTCAAATTGCCCATGAAACGCCGCGGGATAAGGTTCAAAAACATAATGGCTAATTTTATGTTTAATTTTTGATTACCCAGCCCTTGTAGATGCCGCGTTAAGTCTAGCTTCGTTAAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTGAACTAAATGTTCCAGAAAGTTTGGAAGAAGATGTTCCATTGt                                                                                                                                                 >  1:335143/1‑250 (MQ=255)                                                                                               ggaaatggaaTCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAAGATCCAGATGTTTGGTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTa                                                                                                                                                <  2:132900/250‑1 (MQ=255)                                                                                                               ggaccatggagagccgcGGGATAAAGTTCAGAAACATCATGTCTAGTTTTATGATTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCTGTCAAATCCATTAAAGTTATCTGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGTTGAATGTTAACATATTATTCCAGAAAGTGTGGAATAAGGTTTTCCTTTGTAGGACAtacacattca                                                                                                                                 >  2:334576/13‑240 (MQ=255)                                                                                                                        aaCGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTTTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACCTCCCCATGCAAGCCGATTTg                                                                                                                       <  1:340411/250‑1 (MQ=255)                                                                                                                               cgGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGAGGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTccc                                                                                                               <  1:208570/251‑1 (MQ=255)                                                                                                                                     aaGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTTCGTGAGTAAGGTAAAACCCAGTTCCAACACCAACATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTGCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAAcc                                                                                                         <  1:436002/251‑1 (MQ=255)                                                                                                                                            aaaaaCATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGt                                                                                                  >  1:433064/1‑251 (MQ=255)                                                                                                                                                    catGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGc                                                                                          >  1:33629/1‑251 (MQ=255)                                                                                                                                                     atGGCTAATTTTATGTTTACTTTTTAATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGCTACTCCCAAAACCCACCAGTCATaa                                                                                             >  2:79198/1‑247 (MQ=255)                                                                                                                                                       ggCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTACAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCCACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTATGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCttt                                                                                       >  2:418353/1‑251 (MQ=255)                                                                                                                                                           aaTTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAggg                                                                                   <  2:64950/251‑1 (MQ=255)                                                                                                                                                               ttATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAgggga                                                                                 >  2:68240/1‑247 (MQ=255)                                                                                                                                                                tATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGa                                                                                                                                                                                                                                       <  1:444138/98‑1 (MQ=255)                                                                                                                                                                tATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGa                                                                                                                                                                                                                                       >  2:444138/1‑98 (MQ=255)                                                                                                                                                                     ttAATTTTTGATTCCCCAACCCTTGTAGATGCCGGGGTCAGGCTAGCTTCTTTTAAATCCTGTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAATCCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCGTTGGAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAgtggtg                                                                          <  2:329416/250‑1 (MQ=255)                                                                                                                                                                          ttttGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCACAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTaa                                                                     <  1:14591/250‑1 (MQ=255)                                                                                                                                                                            ttGATTCCCCAGCCCTTGGAAAGGCCGCGTTCAGGCTAGCTTCGGTGAAATCCATTAAAGAAAACAAATTACACTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATAGCCAGATGGTTGTTACCTAAATGTCACATAAAGTGTGGAAGAAGTTGTGCCTATGTAGGACATCCCCAGGCAAGCCGATTTGATACGCACAAAACCCGCCAGGCTTAAAGCTTTAGGGTAAGTTGTGTGTAAAt                                                                   <  2:499584/250‑1 (MQ=255)                                                                                                                                                                               tttCCCCAGCACGTGTAGTTGCCGCTTTCAGGCTAGCTTCGTACAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGGTCCAGAAAGTGTGGAAGAAGCTGTTCCTGTGTAGTACACCCCCATGCAAGCCGATTTGATACTCCCAATACCCACCAGACATAAAGCTTTAGGGTAAGGGGTGTGTAAATTCt                                                                <  2:470891/249‑1 (MQ=255)                                                                                                                                                                                   cccAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAACCCATTAAAGATATCAGATTACTCGCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAAGTCTAGccc                                                           <  2:302144/251‑1 (MQ=255)                                                                                                                                                                                    ccAGCCCGTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAATTTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTTTGGAAGAAGGTGTGCCTTTGTAGTACATCCCCATGCAAGCCGATTTGATTCCCCCAAAACCCACCAGTCATAAAGCTTTATGGTAAGTGGTGTGTAAATTCTAGcccc                                                          <  1:439257/251‑1 (MQ=255)                                                                                                                                                                                        cccTTGTAGATGCCGCGTGCAGGCTATCTTCGTTAAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCAAGTTCCAACACCAATATCCAGATGTTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTTTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCGTTAGGGTAAGTGGTGTGTAAATTCTAGCCCcatca                                                      <  2:184320/251‑1 (MQ=255)                                                                                                                                                                                                               cTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAACTTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTc                                >  2:160678/1‑250 (MQ=255)                                                                                                                                                                                                                          cAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTGCCAACACAAAAATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCAGGCAAGCCGATTTGATACTCCCAAAACCCCCCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCTCCATGTTATAg                    <  2:163115/251‑1 (MQ=255)                                                                                                                                                                                                                            aaaccATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATcc                                                                                                                                        <  2:122317/130‑1 (MQ=255)                                                                                                                                                                                                                            aaTCCATTAAAGATCTCAGATTACTCTCAGGCCCGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATcc                                                                                                                                        >  1:122317/1‑133 (MQ=255)                                                                                                                                                                                                                               ccaatAAAGATATCATATTACTCTCAGGTACTTGAGTAAGGTAAAACCCAGGACCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTTTTCCCTTGTAGGACATCCCCATGCAAGGCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTATGTGGTGTGTAAATTCTAGCCCCATCATCGGTGTTGTTTTATTAATTTCACCATGTTATAGTTTTa               <  2:208358/247‑1 (MQ=255)                                                                                                                                                                                                                                    aaaGATATCAGATTACTCTCAGGTACATGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTtgt          <  2:33629/251‑1 (MQ=255)                                                                                                                                                                                                                                        atatCAGATTACTCTAAGGTACGTGAGTAAGGTAAAACCCCGATCCAACCCCAATATCCAGCAGGTTGTTACCTAAATGTCCCAGAACGTGTGGAAGAAGGTGATCATTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAACACCCACCAGTCATAACACTTTAGGGTAAGTGGTGTGTAAATACTATCCCCGTAATCTGTTTTTTTTTCTTAATTTCACCTTGTTCTAGTTTTATTTGTGAAt      >  2:294871/1‑251 (MQ=255)                                                                                                                                                                                                                                           acAGATTACTCTCAGTTACGTGAGTAAGGTAAA‑CCCAGTACCAACACCAAGATACAGATGGTTGGTACCTAAATGTTCTAGAAAGTGTGGAAGAAGGTGTTCCGTCGTAGGACAACCCCATGCAAGCCGAGTTGATACTACCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCTCATCATCTGTGTTTATTTATTAATTTCACCATGTTAGAGTTTCATTTGTGAATTaaa  <  2:415582/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GAGTTCCATCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTATTTGTGAATTAAA  >  NC_000913/582358‑582843

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