Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,027,754	G→A	A340V (GCG→GTG)	fadA ←	3‑ketoacyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,027,754	0	G	A	94.3%	161.6 / ‑6.4	54	A340V (GCG→GTG)	fadA	3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (27/23);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGCACACGAGCGCGCGGCTTAAGACCT  >  NC_000913/4027508‑4027996                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ggTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGAt                                                                                                                                                                                                                                                <  1:140500/251‑1 (MQ=255) ggTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGAt                                                                                                                                                                                                                                                <  2:87828/251‑1 (MQ=255) ggTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGAt                                                                                                                                                                                                                                                >  2:854/1‑251 (MQ=255)        gatgccataTGCGACGCGTGCGTGTCTTATCAGGCCTACATTGGTGCCGGATCGGGAGACCGGGTAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATTCCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGGGGATGACCCAGCGCGATCGCGcc                                                                                                                                                                                                                                          <  2:242965/242‑1 (MQ=255)           gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCg                                                                                                                                                                                                                                      >  1:131847/1‑251 (MQ=255)                  gCGACGCTGGCGCGTCTTAGCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGt                                                                                                                                                                                                                               <  1:324898/251‑1 (MQ=255)                    gACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTg                                                                                                                                                                                                                             >  1:140856/1‑251 (MQ=255)                       gCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGAtc                                                                                                                                                                                                                          <  1:346902/251‑1 (MQ=255)                                    aTCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCt                                                                                                                                                                                                             >  1:278277/1‑251 (MQ=255)                                      cAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCTTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTg                                                                                                                                                                                                            <  2:138163/250‑1 (MQ=255)                                      cAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGc                                                                                                                                                                                                           >  2:402095/1‑251 (MQ=255)                                        ggCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTc                                                                                                                                                                                                         >  2:86212/1‑251 (MQ=255)                                           cTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAAt                                                                                                                                                                                                      >  2:354825/1‑251 (MQ=255)                                             aCATTGGTGCCGGATCGGTAGACCGGATAGGGCCTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATCCCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTa                                                                                                                                                                                                    <  2:9200/251‑1 (MQ=255)                                                     gCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGa                                                                                                                                                                                            <  2:293938/251‑1 (MQ=255)                                                      ccGGATCGGAAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATTCGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGAACTTCTCGTCAATCTGCTCAATTAGTCCCAGAt                                                                                                                                                                                           <  2:131071/251‑1 (MQ=255)                                                       cGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATc                                                                                                                                                                                          <  1:329026/251‑1 (MQ=255)                                                        ggATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCt                                                                                                                                                                                         >  1:207717/1‑251 (MQ=255)                                                        ggATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCt                                                                                                                                                                                         >  1:92405/1‑251 (MQ=255)                                                            cGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTaa                                                                                                                                                                                     <  1:16373/251‑1 (MQ=255)                                                             ggTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAAt                                                                                                                                                                                    <  1:365284/251‑1 (MQ=255)                                                                atacaggATAAGGCTTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACGCCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGATCCGCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGATGATACGCGCACCCGAACAACCCAGCGGATGTCCCAGCACGATCGCGCCGCCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGGTCTTTACTACa                                                                                                                                                                                 <  2:350000/246‑1 (MQ=255)                                                                             cGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAAGTAGTCCCAGATCTTTAATACATGTCAGGATCTgc                                                                                                                                                                    <  1:351177/251‑1 (MQ=255)                                                                                 cACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGcgcg                                                                                                                                                                 <  1:178624/250‑1 (MQ=255)                                                                                   cgccgcATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGc                                                                                                                                                               >  1:321266/1‑250 (MQ=255)                                                                                   cgccgcATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGGCAATCGGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCa                                                                                                                                                              <  1:416359/251‑1 (MQ=255)                                                                                                aaGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGCTCTGCGCGGCAACGGCTTCGTTc                                                                                                                                                  >  1:52175/1‑250 (MQ=255)                                                                                                    ggTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTc                                                                                                                                             <  1:211931/251‑1 (MQ=255)                                                                                                         aacccgcaccaccacCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCCGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCACTCTGCTCAATCAGCCCCAGACCTTTAATACATGGAAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAac                                                                                                                                         >  2:406457/13‑250 (MQ=255)                                                                                                             cGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGATGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCACAGGCTTCGTTCAtttaaaccanggc                                                                                                                                    >  2:382420/1‑241 (MQ=255)                                                                                                                     cacCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGc                                                                                                                            >  1:76252/1‑251 (MQ=255)                                                                                                                           cGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGc                                                                                                                        >  2:32860/1‑249 (MQ=255)                                                                                                                                  cccTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACTCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGccc                                                                                                               >  2:178654/1‑251 (MQ=255)                                                                                                                                        cccAGACCGTTACACATCGTCGCCAGCCCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGTTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACTATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTTAATTAGTCCCAGACCTTAATTACCTGGCAGTATCTTCTCGCCAAAGGCTTCGTCCATTTAAACCACGCCGATATCGCTGGCAGAAAGCCCCGCttt                                                                                                         <  1:150721/251‑1 (MQ=255)                                                                                                                                         ccagtccgtaacacacCGCCGCCAGACCAAACTAAACGTCTTTCCGTTCCCTCAGATTCGCCAGCGGGGTGCTGATCCGCGCACCGGAACAACCCAGCGGAGGACCCCGCCGGATCGCGCCCCCTTTGAGTTGGATCTTCTCGTCAATCTGCTCAATTATCCCCAGATCTTTAATACATGGCAGGATCTGCGCGCCAAAGGCTTCTTTCATTTCAAACACGCCGATATCGCTGGCAGAAATCCCCGCtttt                                                                                                        <  1:382420/241‑1 (MQ=255)                                                                                                                                          cAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCttttt                                                                                                       >  2:65316/1‑251 (MQ=255)                                                                                                                                          cAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGATCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCttttt                                                                                                       <  1:127702/251‑1 (MQ=255)                                                                                                                                               cGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGCTACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCATTTTTCAgc                                                                                                  >  2:231040/1‑251 (MQ=255)                                                                                                                                                 aTACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAgcgc                                                                                                >  1:57666/1‑251 (MQ=255)                                                                                                                                                 aTACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGAGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCCCtttttcaagcgc                                                                                                >  1:206446/1‑244 (MQ=255)                                                                                                                                                   acacATCGTCGCCAGACCAAACTGACCGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCCCCGGAACAACCCAGCGGATGCCCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCa                                                                                              <  1:32860/251‑1 (MQ=255)                                                                                                                                                     acaTCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGt                                                                                            >  1:146262/1‑251 (MQ=255)                                                                                                                                                         cgtcgCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATCCGCGCACCGGAACAACCCAGCGGATGCCCCAGCGCGATCGCGCCGCCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCTTTCATTTCAACCACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCg                                                                                        <  2:6313/251‑1 (MQ=255)                                                                                                                                                              ccAGACCAAACTGAACGTTTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCGGATCCGCGCCCCGGACAAACCCAGCGGATGCCCCAGCACGATCGCGCCACCGTTGAGGTTGATTTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGcc                                                                                   <  2:214139/251‑1 (MQ=255)                                                                                                                                                                                        tCCATCAGATTCAGCCGGAGGGTGCTGATGCCCTCACCGGAACAACCCCGTGCACGACCCGGCGCGCTCGCGACCCGGGTGAGGTTGACCTTCTCGTCAATCTGCACAATTGGTCCCAGATGTTTAATTCATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCTTGTCAATCGCGCCGATATCGCTGGCAGAAAGCCCCGCATTTTTCAGCGCCAGTTTCGAGGCCGGACCCGGGCCGTAACCCATAATCg                                                          <  2:86413/250‑1 (MQ=255)                                                                                                                                                                                                   cagcagCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaac                                              <  1:305268/251‑1 (MQ=255)                                                                                                                                                                                                       agcGTGGTGCTGATACGCGCACCGGAACACCCCAGCGGATGACCAAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCAcaaccaac                                          <  1:354825/251‑1 (MQ=255)                                                                                                                                                                                                               gCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGCTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGcc                                   >  1:22733/1‑250 (MQ=255)                                                                                                                                                                                                                 tGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATc                                >  2:236385/1‑251 (MQ=255)                                                                                                                                                                                                                       gcgcACCGGAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGAAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGc                          >  1:411056/1‑251 (MQ=255)                                                                                                                                                                                                                               gAACAACCCAGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGCACACGAgc                  >  2:125569/1‑251 (MQ=255)                                                                                                                                                                                                                                        aGCGGATGACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGCACACGAGCGCGCGGCtt         >  1:330757/1‑251 (MQ=255)                                                                                                                                                                                                                                             atgccccAGCACTATCGCGCCACCTTTGATGTTGATCTTCTCGTCAATCTCCTCAATTAGTCCCAGATCTGAAATACATGGCAGGATCTGCGCGCCAAAGGCTTGTTTCATTTCAACCACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGCACACGAGCGCGCGGCCTAAGAc    <  1:154124/247‑1 (MQ=255)                                                                                                                                                                                                                                               gACCCAGCACGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGCACACGAGCGCGCGGCTTAAGNCCt  >  2:379700/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATCTTCTCGTCAATCTGCTCAATTAGTCCCAGATCTTTAATACATGGCAGGATCTGCGCGGCAAAGGCTTCGTTCATTTCAAACACGCCGATATCGCTGGCAGAAAGCCCCGCTTTTTTCAGCGCCAGTTTCGAGGCCGGAACCGGGCCGTAACCCATAATCGATGGGTCACAACCAACGACCGCCATCGAACGCACACGAGCGCGCGGCTTAAGACCT  >  NC_000913/4027508‑4027996

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