Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,604	A→C	6.1%	L478R (CTG→CGG)	PP_3486 ←	cytochrome c

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,949,604	0	A	C	6.1%	258.0 / 5.2	81	L478R (CTG→CGG)	PP_3486	cytochrome c
Reads supporting (aligned to +/- strand):  ref base A (37/39);  new base C (3/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTT  >  NC_002947/3949459‑3949752                                                                                                                                                   |                                                                                                                                                     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                               >  3:217324/1‑137 (MQ=255)                                                                                                                                  tAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGc                               <  4:217324/137‑1 (MQ=255)                                                                                                                                       ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                           >  7:139597/2‑120 (MQ=255)                                                                                                                                       ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                           >  7:133849/2‑120 (MQ=255)                                                                                                                                       ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                           <  8:133849/119‑1 (MQ=255)                                                                                                                                       ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCgca                                           <  8:139597/119‑1 (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)                                                                                                                                                   aGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCtt  <  8:165875/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCT‑GGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTCCAATCTGCTT  >  NC_002947/3949459‑3949752

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