Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,630,703	C→A	6.3%	A427A (GCC→GCA)	PP_3200 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,630,703	0	C	A	6.3%	256.6 / 5.1	79	A427A (GCC→GCA)	PP_3200	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (3/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CGGTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGAT  >  NC_002947/3630562‑3630852                                                                                                                                              |                                                                                                                                                      cGGTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGa                                                                                                                                               <  1:1349223/150‑1 (MQ=255) cGGTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGGCTCCAGGCTGGCAAGCCAGGa                                                                                                                                               <  2:661802/150‑1 (MQ=255)        ttGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCAGGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAAc                                                                                                                                        >  1:519744/1‑150 (MQ=255)        ttGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAAc                                                                                                                                        <  2:2056710/150‑1 (MQ=255)        ttGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGACACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAAc                                                                                                                                        <  1:1939614/150‑1 (MQ=255)           tCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTc                                                                                                                                     >  1:2104322/1‑150 (MQ=255)             aaCCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCAt                                                                                                                                   >  2:1650184/1‑150 (MQ=255)                ccGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGc                                                                                                                                                    >  1:232681/1‑130 (MQ=255)                ccGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGc                                                                                                                                                    <  2:232681/130‑1 (MQ=255)                  gCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTg                                                                                                                              >  1:1861642/1‑150 (MQ=255)                   cTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGc                                                                                                                             >  1:488222/1‑150 (MQ=255)                         ccAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTc                                                                                                                       <  2:488222/150‑1 (MQ=255)                          cAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCACTGCCGAGTCg                                                                                                                      >  1:1098188/1‑150 (MQ=255)                           aGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGa                                                                                                                     >  2:2097971/1‑150 (MQ=255)                            gCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCAGGCTGGCAAGCCAGGAgacgttg                                                                                                                                        >  2:951278/1‑124 (MQ=255)                            gCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCAGGCTGGCAAGCCAGGAgacgttg                                                                                                                                        <  1:951278/130‑7 (MQ=255)                               gTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCtt                                                                                                                 >  1:1106415/1‑150 (MQ=255)                                ttGGCCTCGACCAACCAGCATGCGTTGGCCGTGATCTGGTCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCCGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTc                                                                                                                <  2:1398635/150‑1 (MQ=255)                                    ccacGAGCAACCGCCATTCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCg                                                                                                            >  2:2401724/4‑150 (MQ=255)                                      tCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGc                                                                                                          <  2:136277/150‑1 (MQ=255)                                      tCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGc                                                                                                          >  1:1343112/1‑150 (MQ=255)                                          ccaaccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGaaa                                                                                                      <  2:1861642/150‑1 (MQ=255)                                          ccaaccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACCTCATCGCTGCCGAGTCGAGCTTCCTCGGCGaaa                                                                                                      >  1:1697261/1‑150 (MQ=255)                                           caaccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAAt                                                                                                     >  1:908653/1‑150 (MQ=255)                                             accaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCg                                                                                                   >  1:2161010/1‑150 (MQ=255)                                             accaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCGTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGATCTTCCTCGGCGAAATCg                                                                                                   >  1:398827/1‑150 (MQ=255)                                                agaattcGTTGGCGGTGGCTGAGGCAATCTATAGTGATGATAACCCGGTTTATGCGACCGGCCGCACGCTGCAGGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCGGAGTCGAGCTTCCTCGGCGAAATCGCCt                                                                                                <  2:141415/144‑1 (MQ=255)                                                 gCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTg                                                                                               <  1:2097971/150‑1 (MQ=255)                                                        ttGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                        <  1:2401724/150‑1 (MQ=255)                                                        ttGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                        <  2:418174/150‑1 (MQ=255)                                                        ttGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                        <  1:131422/150‑1 (MQ=255)                                                        ttGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACCTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                        <  2:1697261/150‑1 (MQ=255)                                                         tGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGTTCGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCggg                                                                                       <  1:637035/150‑1 (MQ=255)                                                         tGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCggg                                                                                       <  2:2058019/150‑1 (MQ=255)                                                                 ggcgaggcgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGc                                                                               <  2:1343112/150‑1 (MQ=255)                                                                  gcgaggcgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCg                                                                              >  2:2066397/1‑150 (MQ=255)                                                                     aggcgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGa                                                                           <  2:2124319/150‑1 (MQ=255)                                                                       gcgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTTGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGc                                                                         >  2:241154/1‑150 (MQ=255)                                                                         gaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCaa                                                                       <  1:2066397/150‑1 (MQ=255)                                                                           tCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTg                                                                                                                                                 <  1:281779/74‑1 (MQ=255)                                                                           tCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTg                                                                                                                                                 >  2:281779/1‑74 (MQ=255)                                                                           tCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTTGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGa                                                                     <  1:241154/150‑1 (MQ=255)                                                                           tCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGa                                                                     <  1:1650184/150‑1 (MQ=255)                                                                             gACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAAc                                                                   <  2:908653/150‑1 (MQ=255)                                                                                  tGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGc                                                                               <  2:275400/133‑1 (MQ=255)                                                                                  tGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGc                                                                               >  1:275400/1‑133 (MQ=255)                                                                                  tGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGa                                                              >  1:2076723/1‑150 (MQ=255)                                                                                     cGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAATGCCGATGc                                                           <  1:2054673/150‑1 (MQ=255)                                                                                      gATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGcc                                                          >  2:2172986/1‑150 (MQ=255)                                                                                      gATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACCCCGATGcc                                                          >  2:2037674/1‑150 (MQ=255)                                                                                            ccGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGAc                                                    >  1:1786102/1‑150 (MQ=255)                                                                                             cGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGTCGATCTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGAcc                                                   >  2:2133218/1‑150 (MQ=255)                                                                                             cGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGAcc                                                   >  1:497245/1‑150 (MQ=255)                                                                                               gTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGCTTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGTGCAAGAACGCCGATGCCGTTGActct                                                 >  2:839457/1‑150 (MQ=255)                                                                                                     gCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGcc                                           >  2:2205941/1‑150 (MQ=255)                                                                                                      ccACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGt                                                                                      >  1:231053/1‑106 (MQ=255)                                                                                                      ccACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGt                                                                                      <  2:231053/106‑1 (MQ=255)                                                                                                       cACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGa                                         >  1:2006188/1‑150 (MQ=255)                                                                                                                         cGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGt                                                                             <  2:1651746/96‑1 (MQ=255)                                                                                                                         cGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGt                                                                             >  1:1651746/1‑96 (MQ=255)                                                                                                                            gttcggctcCAGGCTGGCAAGCCAGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTgcg                    >  1:678863/7‑150 (MQ=255)                                                                                                                                gTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTg                <  2:113972/150‑1 (MQ=255)                                                                                                                                gTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTg                <  2:2104322/150‑1 (MQ=255)                                                                                                                                gTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCCTCCTCGGCGCAATCGCCTGGACCCGGGTGATGAGCGTAAGCAAGTACGCCGATGCCGTTTTCCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTg                <  2:2080693/150‑1 (MQ=255)                                                                                                                                 gctcCAGGCTGGCAAGCCAGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTgg               >  1:1274160/2‑150 (MQ=255)                                                                                                                                   tcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                        <  1:891151/75‑1 (MQ=255)                                                                                                                                   tcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                        >  2:891151/1‑75 (MQ=255)                                                                                                                                   tcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGc                                <  2:2140928/131‑1 (MQ=255)                                                                                                                                   tcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGc                                >  1:2140928/1‑131 (MQ=255)                                                                                                                                   tcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGgtg             >  1:1337321/1‑150 (MQ=255)                                                                                                                                   tcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGgtg             >  1:1322547/1‑150 (MQ=255)                                                                                                                                     cTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGtt           >  1:1769650/1‑150 (MQ=255)                                                                                                                                      tGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc          <  2:1322547/150‑1 (MQ=255)                                                                                                                                      tGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc          <  2:1337321/150‑1 (MQ=255)                                                                                                                                      tGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc          >  1:1658767/1‑150 (MQ=255)                                                                                                                                      tGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc          >  1:2078900/1‑150 (MQ=255)                                                                                                                                      tGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc          <  1:2259162/150‑1 (MQ=255)                                                                                                                                      tGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc          >  1:254397/1‑150 (MQ=255)                                                                                                                                            gCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCg    <  2:1098188/150‑1 (MQ=255)                                                                                                                                             ccAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGa   <  1:788798/150‑1 (MQ=255)                                                                                                                                             ccAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCAACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGa   >  2:1535090/1‑150 (MQ=255)                                                                                                                                              cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGGTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGc                                                           <  2:1565231/93‑1 (MQ=255)                                                                                                                                              cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGc                                                           >  1:1565231/1‑93 (MQ=255)                                                                                                                                              cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGt                                     >  1:966874/1‑115 (MQ=255)                                                                                                                                              cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGt                                     <  2:966874/115‑1 (MQ=255)                                                                                                                                              cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGAt  >  1:717092/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                      CGGTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGAGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGAT  >  NC_002947/3630562‑3630852

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