Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,299,162	G→A	Y139Y (TAC→TAT)	napH ←	ferredoxin‑type protein NapH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,299,162	0	G	A	97.7%	152.0 / NA	44	Y139Y (TAC→TAT)	napH	ferredoxin‑type protein NapH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/21);  minor base C (0/1);  total (22/22)
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.99e-01

GGCCGCAACGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACGTAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAT  >  NC_000913/2298919‑2299402                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 ggCCGCAACGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGa                                                                                                                                                                                                                                            >  1:351533/1‑250 (MQ=255)       aaCGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGc                                                                                                                                                                                                                                     <  1:305701/251‑1 (MQ=255)        aCGGTAATCACACCTTTGCTACCCAGCACTCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGAGCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCg                                                                                                                                                                                                                                    <  2:293926/251‑1 (MQ=255)             aaTCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCCTTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCACATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAggg                                                                                                                                                                                                                                >  2:315732/1‑250 (MQ=255)                            aCCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGGCCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTg                                                                                                                                                                                                                  <  1:197047/249‑1 (MQ=255)                                                   cAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGCTCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAAc                                                                                                                                                                                         >  2:193658/1‑251 (MQ=255)                                                     aCCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCa                                                                                                                                                                                       <  2:337096/251‑1 (MQ=255)                                                       cGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGt                                                                                                                                                                                     <  1:353145/251‑1 (MQ=255)                                                              aTGTGACCGCACCAGACGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATAAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTa                                                                                                                                                                               >  2:28817/1‑250 (MQ=255)                                                                  gCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTc                                                                                                                                                                           <  1:32891/250‑1 (MQ=255)                                                                        aCCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCGAAGTCGGTAATc                                                                                                                                                                    <  1:315732/251‑1 (MQ=255)                                                                           aGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGAGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCCACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTACTCggg                                                                                                                                                                 >  1:401810/1‑251 (MQ=255)                                                                                 gtTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGc                                                                                                                                                           >  2:332571/1‑251 (MQ=255)                                                                                    cAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCggg                                                                                                                                                        <  2:351533/251‑1 (MQ=255)                                                                                           cAGTAGATCAAATAACAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGCTTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGAccc                                                                                                                                                 >  1:48006/1‑251 (MQ=255)                                                                                                gATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGGTCAGCGGGCAGACCCAGCTg                                                                                                                                            <  1:155314/251‑1 (MQ=255)                                                                                                gATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGGAATCGGGTTCAGCGGGCAGACCCAGCTg                                                                                                                                            >  1:183556/1‑251 (MQ=255)                                                                                                                 gcgAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTcc                                                                                                                           >  1:287750/1‑251 (MQ=255)                                                                                                                     gaATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGcc                                                                                                                       >  1:216903/1‑251 (MQ=255)                                                                                                                          agcagcGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGc                                                                                                                  >  1:343615/1‑251 (MQ=255)                                                                                                                              gcgcgcCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGt                                                                                                                >  1:94482/1‑249 (MQ=255)                                                                                                                              gcgcccCGCTGCCGAAGCCCATTACCAGGCTACGGCCCATCAGAGAAACGGGGGTGATCCATTCCCAAAGGAGCGTGCCAGTGAAGGCTGATCCCACCAGAAGGACCACCAACAGCACATAGCGAAGGTGGCGAGGGAGCGGCGCAGACTGAGTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCGAGGGCTGTAATCGGGTTCAGCGGGCAGGCCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTag                                                                                                              <  2:196869/251‑1 (MQ=255)                                                                                                                                 cgccgcTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTgg                                                                                                                                                                                                                                      <  1:11082/128‑1 (MQ=255)                                                                                                                                 cgccgcTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTgg                                                                                                                                                                                                                                      >  2:11082/1‑128 (MQ=255)                                                                                                                                    cgcTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACg                                                                                                        >  1:377476/1‑251 (MQ=255)                                                                                                                                        gCCGAAGCCCATTTCCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTg                                                                                                     <  2:401856/250‑1 (MQ=255)                                                                                                                                            aaGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAAt                                                                                                >  1:78730/1‑251 (MQ=255)                                                                                                                                                    tACCAGGCTACGCCCCATCAGAGAAAAGGGGTTGATCCATTCCCAAATGAGCGAGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGAACTCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATcaccgcac                                                                                        >  2:286353/1‑251 (MQ=255)                                                                                                                                                       cAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCgg                                                                                     >  1:242113/1‑251 (MQ=255)                                                                                                                                                        aGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGt                                                                                    <  2:287750/251‑1 (MQ=255)                                                                                                                                                             aCGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAAc                                                                                <  1:28817/250‑1 (MQ=255)                                                                                                                                                                    aTCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCATTTAAGGCTTATCCCACCAGAATGACCACCAACAGCACATAGCGAATGGGGCGAGGGATCGTCGCAGCCTGATTGAGGTCAAAGCGCCTGCGTAACCATTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGt                                                                        <  1:205668/251‑1 (MQ=255)                                                                                                                                                                             aCGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACAGAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGgcaggcag                                                               <  1:349987/251‑1 (MQ=255)                                                                                                                                                                                gggTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGGGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCCGATg                                                            >  1:416508/1‑251 (MQ=255)                                                                                                                                                                                gggTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTGCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCTCCGGTCAACGCCACCGTGGCAGGCAGATg                                                            <  2:147205/251‑1 (MQ=255)                                                                                                                                                                                    tGATCCATTCCCAAAGGAGCTGGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCGAAGTCGGGAATCGGGTTCAGCTGGCAGACCCAGCTGCAAAATAATCGCTGTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGTACCGTTCAACGCCACGGTGTCAGGCAGGTGAcc                                                         <  2:214785/250‑1 (MQ=255)                                                                                                                                                                                                             aGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCCGATGACCAATGGCAAGACTTTGCAGCGgcatcc                               >  2:151492/1‑250 (MQ=255)                                                                                                                                                                                                              gTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAg                              <  1:71806/251‑1 (MQ=255)                                                                                                                                                                                                                    gCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATc                        <  1:17655/251‑1 (MQ=255)                                                                                                                                                                                                                         tCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGCTCGGTc                    >  1:387412/1‑250 (MQ=255)                                                                                                                                                                                                                           ccaccaGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGATGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAAGCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGc                 <  1:256224/251‑1 (MQ=255)                                                                                                                                                                                                                               caGAATGACCACCAACAGCACCTAGCGAATGTGGCGAGGGATCGGCGCAGACTGATTGAGGTCAAAGCGCCGGCGGAACCAGTTCGCTAAGGCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTCCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGaa             <  2:39056/251‑1 (MQ=255)                                                                                                                                                                                                                                  aaTGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGaaagg          >  2:218019/1‑248 (MQ=255)                                                                                                                                                                                                                                          ccaACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTAAACGCCACCGGGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGAGTCgcat  >  1:4359/1‑248 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GGCCGCAACGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACGTAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAAT  >  NC_000913/2298919‑2299402

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