Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,602	G→C	7.2%	P479A (CCC→GCC)	PP_3486 ←	cytochrome c

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,949,602	0	G	C	7.2%	255.8 / 7.7	80	P479A (CCC→GCC)	PP_3486	cytochrome c
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base C (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGC  >  NC_002947/3949455‑3949750                                                                                                                                                     |                                                                                                                                                     cATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCggg                                                                                                                                                      >  2:8583/1‑148 (MQ=255)  aTTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCggg                                                                                                                                                      >  6:233393/1‑147 (MQ=255)     gACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGTCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAgg                                                                                                                                                  <  5:341774/148‑1 (MQ=255)     gACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAgg                                                                                                                                                  <  1:271742/148‑1 (MQ=255)     gACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAgg                                                                                                                                                  <  4:321351/148‑1 (MQ=255)                 cgGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAggg                                                                                                                                                 <  8:93245/137‑1 (MQ=255)                 cgGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAggg                                                                                                                                                 >  7:93245/1‑137 (MQ=255)                 cgGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGccacca                                                                                                                                     <  7:39578/149‑1 (MQ=255)                  ggCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCcaccac                                                                                                                                    <  5:385334/149‑1 (MQ=255)                       aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATggg                                                                                                                               >  5:120998/1‑149 (MQ=255)                       aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATggg                                                                                                                               <  3:222728/149‑1 (MQ=255)                       aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAAGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATggg                                                                                                                               <  6:207667/149‑1 (MQ=255)                            cccAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCt                                                                                                                          >  5:378375/1‑149 (MQ=255)                             ccAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTg                                                                                                                         >  1:212196/1‑149 (MQ=255)                             ccAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTg                                                                                                                         <  5:70463/149‑1 (MQ=255)                                  tCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATg                                                                                                                    >  3:36654/1‑149 (MQ=255)                                   cACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGc                                                                                                                   >  3:138118/1‑149 (MQ=255)                                    aCGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGc                                                                                                                   <  5:233393/148‑1 (MQ=255)                                        tGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTgg                                                                                                              >  7:81287/1‑149 (MQ=255)                                          gTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggcg                                                                                                            >  8:17570/1‑149 (MQ=255)                                          gTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggcg                                                                                                            >  3:213441/1‑149 (MQ=255)                                          gTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTgg                                                                                                              <  7:37776/147‑1 (MQ=255)                                              gCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGa                                                                                                        >  3:192297/1‑149 (MQ=255)                                               cACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGAc                                                                                                       >  4:243264/1‑149 (MQ=255)                                                  cAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTAGGCCGCCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGAc                                                                                                       >  1:99557/1‑147 (MQ=255)                                                  cAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAg                                                                                                    <  3:176517/149‑1 (MQ=255)                                                   aGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAgg                                                                                                   <  2:212196/149‑1 (MQ=255)                                                       aCAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCTAAGTAAAGCGTGGCCCGGCCAGGCCCACCGCAACCACATGGGCACCTGCAATGCCCTGGCGGCGACCAGGGcc                                                                                                >  5:201645/1‑148 (MQ=255)                                                          ggCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAg                                                                                            <  3:243264/149‑1 (MQ=255)                                                             acacCCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGCTCCAAAGcc                                                                                          >  8:106985/1‑148 (MQ=255)                                                             acaACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGccc                                                                                         <  4:213441/149‑1 (MQ=255)                                                              caACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGcc                                                                                          >  5:413052/1‑147 (MQ=255)                                                              caACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCa                                                                                        <  8:81287/149‑1 (MQ=255)                                                               aaCCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAg                                                                                       >  1:204058/1‑149 (MQ=255)                                                                 ccACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTg                                                                                     >  5:75205/1‑149 (MQ=255)                                                                 ccACACGATCCGCATGACAAACTCCGCTCAAGGGCT‑GGCCGGCACTGGACGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAACGCCCAGTg                                                                                     >  8:131405/1‑149 (MQ=255)                                                                   acacGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGcc                                                                                   >  5:47323/1‑149 (MQ=255)                                                                     acGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCgg                                                                                 <  4:138118/149‑1 (MQ=255)                                                                      ctatacGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCggg                                                                                <  6:260701/144‑1 (MQ=255)                                                                           cGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAg                                                                                       <  7:257344/137‑1 (MQ=255)                                                                           cGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAg                                                                                       >  8:257344/1‑137 (MQ=255)                                                                            gCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAgg                                                                          >  3:149493/1‑149 (MQ=255)                                                                            gCATGACAACCTCCGCTCAAGGGCT‑GGCCGGAACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAgg                                                                          >  5:293868/1‑149 (MQ=255)                                                                             cATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGc                                                                         <  6:75205/149‑1 (MQ=255)                                                                                gACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCgg                                                                      >  4:141006/1‑149 (MQ=255)                                                                                      cTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGca                                                                >  7:165875/1‑149 (MQ=255)                                                                                          gCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATc                                                            >  4:94133/1‑149 (MQ=255)                                                                                              aaGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTg                                                         >  3:93123/1‑148 (MQ=255)                                                                                                    t‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCt                                                  >  7:169592/1‑149 (MQ=255)                                                                                                    t‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCt                                                  >  2:378244/1‑149 (MQ=255)                                                                                                      ggCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCtt                                                 <  6:120998/149‑1 (MQ=255)                                                                                                       gCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCt                                                  <  8:169592/147‑1 (MQ=255)                                                                                                       gCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTc                                                <  3:141006/149‑1 (MQ=255)                                                                                                        ccGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCTCCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTc                                                >  1:309226/1‑148 (MQ=255)                                                                                                            cACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAgcg                                           >  5:108383/1‑149 (MQ=255)                                                                                                            cACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAgcg                                           >  7:234094/1‑149 (MQ=255)                                                                                                              cTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         <  4:93123/149‑1 (MQ=255)                                                                                                              cTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         <  2:99557/149‑1 (MQ=255)                                                                                                              cTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         <  6:378375/149‑1 (MQ=255)                                                                                                               tGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTgg                                                                   <  1:59516/122‑1 (MQ=255)                                                                                                               tGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTgg                                                                   >  2:59516/1‑122 (MQ=255)                                                                                                                   cGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAAt                                    <  1:378244/149‑1 (MQ=255)                                                                                                                     gTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCa                                  <  6:413052/149‑1 (MQ=255)                                                                                                                                cctatgTAATGCGTGGCCGGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTcgcg                        <  2:204058/143‑1 (MQ=255)                                                                                                                                 cGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGc                      >  6:313934/1‑149 (MQ=255)                                                                                                                                  gAATTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCt                     <  6:201645/149‑1 (MQ=255)                                                                                                                                   aaGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTg                    >  4:100693/1‑149 (MQ=255)                                                                                                                                    aGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTgg                   <  1:8583/149‑1 (MQ=255)                                                                                                                                      tAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGc                             <  4:217324/137‑1 (MQ=255)                                                                                                                                      tAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGc                             >  3:217324/1‑137 (MQ=255)                                                                                                                                           ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         >  7:133849/2‑120 (MQ=255)                                                                                                                                           ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         <  8:139597/119‑1 (MQ=255)                                                                                                                                           ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         <  8:133849/119‑1 (MQ=255)                                                                                                                                           ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                         >  7:139597/2‑120 (MQ=255)                                                                                                                                           ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTAGCGGCTGGCATCCCg            >  2:409090/2‑149 (MQ=255)                                                                                                                                            gTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCa                 >  3:188300/1‑143 (MQ=255)                                                                                                                                            gTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCa                 <  4:188300/143‑1 (MQ=255)                                                                                                                                             tGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTc          <  2:372622/149‑1 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAAt      <  7:17570/149‑1 (MQ=255)                                                                                                                                                  cGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCTGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATc     <  7:131405/149‑1 (MQ=255)                                                                                                                                                  cGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATc     <  4:192297/149‑1 (MQ=255)                                                                                                                                                  cGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATc     >  4:165616/1‑149 (MQ=255)                                                                                                                                                     gCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGc  <  1:355616/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGC  >  NC_002947/3949455‑3949750

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