Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,827,131	T→G	25.4%	D185A (GAC→GCC)	coaBC ←	bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,827,131	0	T	G	25.4%	67.4 / 5.0	45	D185A (GAC→GCC)	coaBC	bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC
Reads supporting (aligned to +/- strand):  ref base T (26/7);  new base G (10/2);  total (36/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCATTCACGGCGGCTTCCATTTCCAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGAGCACGGCTACAGGCGCAGGTGTAGCCAGACAAAT  >  NC_007779/3826895‑3827357                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     aaGCATTCACGGCGGCTTCCATTTCCAGCGCGGTCATCACATCAACACGGTTAACAAACGGTGGCGTCGGTAGGGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGGCAA‑AATAATGTTCAGATGTTTCAGGGCg                                                                                                                                                                                                                                    >  1:2143920/1‑239 (MQ=255)      ttCACGGCGGCTTCCATTTCCAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGGTTCAGGTTG‑TTGAc                                                                                                                                                                                                                              >  1:1618114/1‑239 (MQ=255)           ggcggcTTCCATTTCCAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGAGGGTTTAGGTCG‑GTTGCGGGCg                                                                                                                                                                                                                         >  1:232673/1‑239 (MQ=255)              ggcTTCCATTTCCAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGaa                                                                                                                                                                                                                       <  2:532597/238‑1 (MQ=255)                         ccAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGGAATCCATAATTGTC‑GCTGGTTCTAGTCGTTTGAC‑GGCGAAAAAA‑gcgccaa                                                                                                                                                                                                           >  2:2242000/1‑235 (MQ=255)                          cAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCCGATGGTTCAAGGCG‑GTGACGGGGGGAAAAT‑GCGCTaac                                                                                                                                                                                                          >  1:1951728/1‑236 (MQ=255)                             cgcgGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGGCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCa                                                                                                                                                                                                      >  1:1363150/1‑239 (MQ=255)                                           caACACTTTTAACAGGCGGGGGCGGCGGTAGTGACACCGGCCCTGATCCCGGGGTGCCGGTCGGCCCCCGAGGGGCGGGGGCGGCGGCGAGAGCAAACCCCATCTGGCCGGC‑GTTGTGATTAGAGATATACCGCCCGGGATGGGGCGGTTCGCGGGGGGGGCGGGCGGTAAT‑CATAATGTTCATATGTTTCAGGTCG‑TTGGGGGGCAAAAAAT‑GCGCTACCGCCATATCCACAATGGtt                                                                                                                                                                                        <  2:1327699/239‑1 (MQ=255)                                            aaCACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CAAAATGTTCAGATGGTTCAAGGCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACCATGGTTa                                                                                                                                                                                       >  1:848121/1‑239 (MQ=255)                                                           cGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGGCG‑TTGACCGGGGGAAAAT‑GCGCCTCCGCCCTAATCCCAAAGGTTAAAGGATctgtgtctt                                                                                                                                                                        >  1:44953/1‑231 (MQ=255)                                                               ggCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGGAAAAT‑GCGCTACCGCCCTATCCACAAAGGTTaa                                                                                                                                                                                      >  1:1528592/1‑221 (MQ=255)                                                                            aaaCCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑TATAAAGTTCAGATGGTTTAGGGCG‑TTGACGGGGGGAAAAA‑GCGCCACCGCCATAACCACAAAGGTTAAAGGGGCGAGGAATAGCCCACGGCCGGTATc                                                                                                                                                        >  2:785125/1‑238 (MQ=255)                                                                                               cGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGa                                                                                                                                                                >  1:1639894/1‑211 (MQ=255)                                                                                                   aCGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGGTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGat                                                                                                                                                           >  1:1885957/1‑212 (MQ=255)                                                                                                        cgcgCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGGCG‑TTGACGGGCGGAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCCCCCAGGACCCATATCACCACAAACCTGACTGCCCCTGTCTGGc                                                                                                                            >  1:948086/1‑238 (MQ=255)                                                                                                        cgcgCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGGCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGGGCATTCGCCCAGGACCGATATCACCACCAGCCTGACCGCCACTGCCTGGc                                                                                                                            >  1:1128711/1‑238 (MQ=255)                                                                                                        cgcgCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGGCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAAAGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCCCAAGCCTGACTGCCACTGTCTGGc                                                                                                                            >  1:444060/1‑238 (MQ=255)                                                                                                                                    aaaaaaaacCCCATCTGGCGGGG‑GGTGTAATTAGAAATAAACCCCCCCGGTTGGAGTGGTTCGCGGGGCGGG‑CGGGGGAAATACATAATGTTCATATGTTTCAGGTGT‑TGGACGGGAAAAAAAC‑GCGCTCCCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGaagcaag                                                                                               <  2:444060/234‑1 (MQ=255)                                                                                                                                     tAGCAAAACCCATCTTGCGGGG‑GGTGTGATTAGAGATATAACGCCCCGGATCGGGCGGTTCCCGCGTGGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGAGGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGaagcaag                                                                                               <  1:2242000/238‑1 (MQ=255)                                                                                                                                            accccttgtcccGGAGGCTGTGTTTAGAATTATAACGCACCTGTTGCGGCGTTCCGCGGTCCGGGCCGCCTGTATT‑ATTATTTTTCATTTTTTTCAGGTCG‑TGGGCGGACAAAAAAC‑CCCCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACctct                                                                                        <  2:981870/229‑1 (MQ=255)                                                                                                                                            aCCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑AATAATGGTCAGATGTTTCCAGTCG‑TTGACGGGCGGAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGCGCAAACCACGCGAAGCAAGCACCCCTa                                                                                       >  1:2219453/1‑239 (MQ=255)                                                                                                                                                aTCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACTGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAA‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGCGCATTCCCCCAGGACCGATAACACCACAAGCCTGACAGCCACTGTCCCGCCCCCCGAAGAGCCAACCACGGGGAGAAAACCCCCCTaaaat                                                                                   >  2:589526/1‑237 (MQ=255)                                                                                                                                                aTCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTTAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCAAGGGAAGCAAGCACCTCCaaa                                                                                     >  1:485886/1‑237 (MQ=255)                                                                                                                                                 tCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTa                                                                                  >  1:1354695/1‑239 (MQ=255)                                                                                                                                                 tCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGAGGTAAA‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGCAGCAAGCCACCCTAAATTa                                                                                  >  1:807248/1‑239 (MQ=255)                                                                                                                                                  cTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCCAGCACCTCTAAATTAt                                                                                 >  1:758373/1‑239 (MQ=255)                                                                                                                                                  cTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGCTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑TCGCTACCGCCATATCCACAATGGTTAACGGATCGAGGATTCGCCCAGGACCCATATCCCCACAAACCTAACTGCCACTGTCTGGc                                                                                                                            >  1:2080811/1‑196 (MQ=255)                                                                                                                                                     gCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAAGTCG‑TTGACGGGCGAAAAAA‑GCGCTACCGCCCTATCCACAATGGTTAACGGATCGGGCATTCGCCCAGGACCGATATcacca                                                                                                                                                    >  1:1957918/1‑169 (MQ=255)                                                                                                                                                          a‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCCCAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCCCCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGGGCAAACCACGGGAACCAAGCACCTCTCAATTAtccgtctct                                                                         >  1:1841105/1‑231 (MQ=255)                                                                                                                                                                     agagATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAA‑AATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAA‑GAGCTACCGCCATATCCACAATGGTTAACGGAACGGGCATTCGCCCCGGACCGCAATCACCACAAACCCGACCGCCACCGGCTGGCCCCCAGAGGAGGAAAACACGGGAAGCAAGCACCCCTAAAATTAGCCGCGTGGCAGCGGCa                                                                >  2:454928/1‑238 (MQ=255)                                                                                                                                                                            taACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTAtg                                                                                >  1:87007/1‑215 (MQ=255)                                                                                                                                                                            taACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACa                                                         >  1:1440709/1‑238 (MQ=255)                                                                                                                                                                                    cGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGAagag                                                                                                                 >  1:141004/1‑174 (MQ=255)                                                                                                                                                                                                ttCACGCGTCGGGCCGGCGGTAAT‑CCTAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTAACGCCATATCCAAAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGc                                                                                                                                               >  1:1954769/1‑132 (MQ=255)                                                                                                                                                                                                ttCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GGGCTACCGCCATATCCACAATGGTTAACGGATCGAGCCTTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGcc                                                                                                                           >  1:1867850/1‑152 (MQ=255)                                                                                                                                                                                                 tCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTTAGGTCG‑TTGACGGGCGAAAAAT‑GCGATACCGCCATATCCACAATGGTTAACGGATTGAGCATTCGCCCAGGACCGATATCACCACAAACCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAACACCTCTaaa                                                                                     >  1:1994825/1‑189 (MQ=255)                                                                                                                                                                                                          gggCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGGCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCCATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGGGCACgggtta                          >  1:1232150/1‑235 (MQ=255)                                                                                                                                                                                                              gtgggctatt‑CATATTGTGAAGATTTTCAAGGTGT‑TGGAGGGACGAAAAAT‑GCCCTCCCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGAGCACGGCTACAGgc                      <  1:589526/230‑1 (MQ=255)                                                                                                                                                                                                              cggcagTAAA‑CAAAATGTTCAGATGTTACGGTGTG‑TTGAGGGGAGAAAAAT‑GCCCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCGGCGTGGCAGCGGCACGGTACAGCTGCTGGTTCATGGCGGGGAGCACGGCTACAgg                       <  1:1284654/238‑1 (MQ=255)                                                                                                                                                                                                               ggcgaTAAT‑CATAATGTTAAATTTTTTGAGGGCG‑TTGGCGGGAAAAAAAG‑CCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGAGCACGGCTACAGgcg                     <  1:721510/239‑1 (MQ=255)                                                                                                                                                                                                                                  ccgATGTTTTAGGTGA‑TTTGCCGGCGAAAAATCGCCTCCCCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGAGCACGGCTACAGGCGCAGGTGTAGCCAGACaa    <  2:807248/237‑1 (MQ=255)                                                                                                                                                                                                                                   aGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGGGCAAGGcg                            >  1:382520/1‑212 (MQ=255)                                                                                                                                                                                                                                   aGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGAACATCTGCTGGTTCATGGCGGGGAGCCGCGCTACAGGCGCCGGGGGAGccccgtccct  >  1:1282005/1‑231 (MQ=255)                                                                                                                                                                                                                                          tCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTaaa                                                                                     >  1:2209683/1‑149 (MQ=255)                                                                                                                                                                                                                                          tCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTa                                                                                       >  1:10913/1‑147 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     AAGCATTCACGGCGGCTTCCATTTCCAGCGCGGTCATCACATCAACACGTTTAACAAACGGTGGCGTCGGTAGTGAAACCGGACCTGATACCAGCGTGACGTTCGCGCCACGACGGGCAGCGGCGGCGGCGATAGCAAAACCCATCTTGCCGGA‑GCTGTGATTAGAGATATAACGCACCGGATCGAGCGGTTCACGCGTCGGGCCGGCGGTAAT‑CATAATGTTCAGATGTTTCAGGTCG‑TTGACGGGCGAAAAAT‑GCGCTACCGCCATATCCACAATGGTTAACGGATCGAGCATTCGCCCAGGACCGATATCACCACAAGCCTGACTGCCACTGTCTGGCCCCCAGATGAGCAAACCACGGGAAGCAAGCACCTCTAAATTATGCTGCGTGGCAGCGGCACGGTACATCTGCTGGTTCATGGCGGGGAGCACGGCTACAGGCGCAGGTGTAGCCAGACAAAT  >  NC_007779/3826895‑3827357

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