Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,296,419	G→A	G202G (GGC→GGT)	ccmB ←	cytochrome c maturation protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,296,419	0	G	A	97.8%	140.5 / NA	46	G202G (GGC→GGT)	ccmB	cytochrome c maturation protein B
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/17);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AATCAGCCGGAGCAAAGCCGAATCCCCAGATCCAGCCGACGGTAAGCACGACCACACTGGCAATTGCCAGCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCACTTGCCAGCCATAAACATCCATT  >  NC_000913/2296179‑2296665                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       aaTCAGCCGGAGCAAAGCCGAATCCCCAGATCCAGCCGACGGTAAGCACGACCCCACTGGCACTTGCCAGCCACGGTATAACCCAGCCACAAATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAACATCAGTAATCCCGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGc                                                                                                                                                                                                                                               >  1:252393/1‑250 (MQ=255)  aTCAGCCGGAGCAAAGCCGAATCCCCAGATCCAGCCGCCGGTAATCACGACCACACTGGCAATTGCCCCCCACCGTATAAACCAGCCACAGACTTGATACAGCCGTGGTGGCATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCCAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCTCCGTCGCACAACGCCTTAATGTGGCGGTACCTGCCAGCa                                                                                                                                                                                                                                              >  2:25175/1‑250 (MQ=255)             cAAAGCCGAATCCCCAGATCCAGCCGACGGTAAGCACGACCACACTGGCAATTGCCAGCCACGGTATAAACCAGCCACAGAGTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTaaa                                                                                                                                                                                                                                   >  2:9208/1‑250 (MQ=255)                     aaTCCCCAGATCCAGCCGACGGTAAGCACGACCACACTGGCAATTGCCAGCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGcc                                                                                                                                                                                                                             >  1:57680/1‑248 (MQ=255)                      acccccAGTTCCAGCCTACGGTAAGCCCGTCCCAACTGGCAATTGCCAGCCACGGTATAAACCCGCCACAGATTTGATACAGCTGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATATTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTTTTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTGCCTGCCAGCAACGCGCCTAAAATTGCCAGAt                                                                                                                                                                                                                         <  1:163091/249‑1 (MQ=255)                        ccccAGATCCAGCCGACGGTAAGCACGACCACACTGGCAATTGCCAGCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCAAGATAc                                                                                                                                                                                                                       >  2:204223/1‑251 (MQ=255)                          ccAGATCCAGCCGACGGTAAGCACGACCACACTGGCAATTGCCATCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCATACTCGAACAAAAATCAGTAATCCAGGGTTATTGAATGCTGAGTCGTAACGCTGCCTCCGTCGCAAAATGACTTAATGTGGCGGTACTTGCCAGCAACGCGCCTAAAATTGCCAGATAccc                                                                                                                                                                                                                     <  1:354449/251‑1 (MQ=255)                                               aacgacgtcacTGGCAATTGCAAGCCAAGGTATAAACCCGCCACAGATTTAATACAGCAGTGGAGGCATCACCAGTTGCTGCAGTGGGTTACACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCCACGTTATTGCATGCTGATTCGAAACACTGCCGCCGTAGCAAAAGAACTTAATGTCGCGGAACCTGACAGCACCGCGCCTAAAATTACCAGATACCCGTCAAAGGGCAACTGCATaa                                                                                                                                                                                                 >  2:335086/9‑249 (MQ=255)                                                      ccacgggCAATTCCCGGCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCGGCCGCCGGCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCcgcg                                                                                                                                                                                         <  2:174930/246‑1 (MQ=255)                                                             cAATTGCCAGCCACGTTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATcg                                                                                                                                                                                  <  1:166633/251‑1 (MQ=255)                                                              aaTTGCCAGCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATcgc                                                                                                                                                                                 >  1:2307/1‑251 (MQ=255)                                                                              ataAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATg                                                                                                                                                                 >  2:139169/1‑251 (MQ=255)                                                                                      gCCACAGATTTGGTACAGCCGTGGTGGGATCGCCAGTGGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATgg                                                                                                                                                         <  2:249439/251‑1 (MQ=255)                                                                                       ccACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATggg                                                                                                                                                        >  1:289888/1‑251 (MQ=255)                                                                                        cacaGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATggg                                                                                                                                                        >  1:324495/1‑250 (MQ=255)                                                                                                tGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGc                                                                                                                                               >  2:106465/1‑251 (MQ=255)                                                                                                tGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGc                                                                                                                                               >  1:207451/1‑251 (MQ=255)                                                                                                   tACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGt                                                                                                                                            >  1:154458/1‑251 (MQ=255)                                                                                                                             aTGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCAcacc                                                                                                                  >  2:61355/1‑251 (MQ=255)                                                                                                                                   tgtTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACcgcg                                                                                                            >  1:159735/1‑251 (MQ=255)                                                                                                                                    gtTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGCGCAGCACACCACCgcgc                                                                                                           >  2:200020/1‑251 (MQ=255)                                                                                                                                         ccACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTaa                                                                                                       <  1:104709/250‑1 (MQ=255)                                                                                                                                               aGTTTCGATACCAGACTCGAACAAAAAGCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCa                                                                                                 <  2:57137/250‑1 (MQ=255)                                                                                                                                                gTTTCGATACCAGACTCAAACAAAAATCAGTAATCCAGCTTTATTGAATTCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACTCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCTTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACt                                                                                               <  1:169311/251‑1 (MQ=255)                                                                                                                                                gTTTCGAGACCAGACGCGAACAAAAATCAGTAATCCAGCGTTATTGAAGGCGGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAAGGTCGGGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACt                                                                                               <  2:376847/251‑1 (MQ=255)                                                                                                                                                 tttCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGAATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTg                                                                                              >  1:339192/1‑251 (MQ=255)                                                                                                                                                                 gAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCgg                                                                              >  1:211056/1‑251 (MQ=255)                                                                                                                                                                  aaaaaaaaTCAGTAATCCAGCGTTATTGAATGCTGATTCGTATCGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAgctcccctccctgt                                                                             <  1:369792/248‑13 (MQ=255)                                                                                                                                                                  aaCAAAAATCATTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCTCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAgctcctccccctgt                                                                             <  1:369774/251‑13 (MQ=255)                                                                                                                                                                    cAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGGCAGCGCCACGCCCGGTGc                                                                           <  1:294683/251‑1 (MQ=255)                                                                                                                                                                             gTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCTGTGCACCGAGaaa                                                                  <  1:108615/251‑1 (MQ=255)                                                                                                                                                                              tAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTACCACTGTCAGCGCCACGCCCGGTGCACCGAGaa                                                                   >  1:47325/1‑249 (MQ=255)                                                                                                                                                                                        ttATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGt                                                         >  2:67591/1‑249 (MQ=255)                                                                                                                                                                                               aggCTGATCCGTAACTGTGCCGCCTTCGCAAAAGCGCTTAATGGCGCTGTACCTGCCAGCAACGCGCCTAAAATTGCCCAATACCCGTCAACGGGCAAATTCATAGAAGCCGGGTGCATCGCGGGTGTGTCAAAGTTGAGGAATGGGATGGTCAGCGGTAACCCCAGTATGTTGAGGTGCAGACGACCGCGCTTAAGTCCTACTGTCATCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTcc                                                <  1:51169/249‑1 (MQ=255)                                                                                                                                                                                                  cTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTccc                                               >  1:255127/1‑249 (MQ=255)                                                                                                                                                                                                        cGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCcagcagc                                         >  1:43968/1‑249 (MQ=255)                                                                                                                                                                                                          tAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCg                                     >  1:341423/1‑251 (MQ=255)                                                                                                                                                                                                           aaCGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGCAACACCCGTATGCTGAGCAGCACCCCACCGCGCTTAAGTCCCACTGCCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGCTCCcagcagccgccg                                    >  2:308310/1‑246 (MQ=255)                                                                                                                                                                                                            aCGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTc                                   <  1:61355/251‑1 (MQ=255)                                                                                                                                                                                                                      cgTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGAAACAACAGTATg                                                                                                                                >  2:211157/1‑148 (MQ=255)                                                                                                                                                                                                                      cgTCGCAAAAGGACTTAATGGCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATg                                                                                                                                <  1:211157/148‑1 (MQ=255)                                                                                                                                                                                                                       gTCGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCAc                        >  1:125286/1‑251 (MQ=255)                                                                                                                                                                                                                         cGCAAAAGGACTTAATGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCACtt                      >  2:67148/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaaGGACTTAAGGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCACTTgcc                   <  2:289888/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaTGTCGCGGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCACTTGCCAGCCATAAAc         <  2:287783/251‑1 (MQ=255)                                                                                                                                                                                                                                             cgGTACCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCATCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCACTTGCCAGCCATAAACATCCAtt  >  1:321/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       AATCAGCCGGAGCAAAGCCGAATCCCCAGATCCAGCCGACGGTAAGCACGACCACACTGGCAATTGCCAGCCACGGTATAAACCAGCCACAGATTTGATACAGCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAGCAGCACACCACCGCGCTTAAGTCCCACTGTCAGCGCCACGCCCGGTGCACCGAGAAAGCCAAGCGTAGGCGTTCCCAGCAGCAGCGTCAGCGCCATCACTTGCCAGCCATAAACATCCATT  >  NC_000913/2296179‑2296665

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: ‑