Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,553,252	A→C	6.2%	E84A (GAG→GCG)	PP_1363 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,553,252	0	A	C	6.2%	257.6 / 9.4	98	E84A (GAG→GCG)	PP_1363	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/52);  new base C (4/2);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

CATGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGG  >  NC_002947/1553103‑1553392                                                                                                                                                      |                                                                                                                                             cATGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGa                                                                                                                                              >  1:620652/1‑150 (MQ=255)        caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAg                                                                                                                                       <  2:1251041/150‑1 (MQ=255)        caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAg                                                                                                                                       >  1:2067176/1‑150 (MQ=255)        caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCCAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAg                                                                                                                                       <  1:837699/150‑1 (MQ=255)          ccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGc                                                                                                                                     <  2:2060908/150‑1 (MQ=255)          ccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGc                                                                                                                                     >  1:928024/1‑150 (MQ=255)            aTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGt                                                                                                                                   <  2:1712991/150‑1 (MQ=255)             tCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGtt                                                                                                                                  >  1:1836335/1‑150 (MQ=255)             tCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGtt                                                                                                                                  <  2:187166/150‑1 (MQ=255)             tCTGAACAATCCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCTGGCCGTTCTTCTACCCTACCGCGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGtt                                                                                                                                  <  1:1566621/150‑1 (MQ=255)               tGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAgcgc                                                                                                                                          <  2:183737/140‑1 (MQ=255)               tGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAgcgc                                                                                                                                          >  1:183737/1‑140 (MQ=255)               tGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGt                                                                                                                                >  2:1184039/1‑150 (MQ=255)                gAACAATTCCGACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTa                                                                                                                               <  1:492009/150‑1 (MQ=255)                gAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTa                                                                                                                               >  2:2155703/1‑150 (MQ=255)                      tGCCTACCTCGCCTACGCCCGCGTCGTTTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGa                                                                                                                         <  1:492460/150‑1 (MQ=255)                         cTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGACCTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAATTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCACGCTTTTTACGTGGAcgc                                                                                                                      >  1:1075550/1‑150 (MQ=255)                          tACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAcgcg                                                                                                                     >  2:1217628/1‑150 (MQ=255)                          tACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAcgcg                                                                                                                     >  2:1750196/1‑150 (MQ=255)                          tACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAcgcg                                                                                                                     >  1:983075/1‑150 (MQ=255)                            ccGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTTATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACgcggc                                                                                                                   <  2:1639082/150‑1 (MQ=255)                              gcgcCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCCgcgc                                                                                                                 >  2:816355/1‑150 (MQ=255)                                 ccTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCaaa                                                                                                              <  1:1434589/150‑1 (MQ=255)                                 ccTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCaaa                                                                                                              <  1:1184039/150‑1 (MQ=255)                                  cTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGCGCTCCAGGcactg                                                                                                                                 <  2:2121992/130‑5 (MQ=255)                                  cTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACAAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGCGCTCCAGGcactg                                                                                                                                 >  1:2121992/1‑126 (MQ=255)                                      gACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGCAGTGTCACTGCTGCGCGGCCCGCCGCTACAGCCCCAGGAGTTGGACGGTGACCCTGCGCAAACCTag                                                                                                         <  2:1594125/150‑2 (MQ=255)                                      gACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTAt                                                                                                         <  2:983075/150‑1 (MQ=255)                                       aCCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATc                                                                                                        >  1:389462/1‑150 (MQ=255)                                        cctcggcgATTGCGGTAGGGTAATGGTAGACATGTATCCGGTTGATAGCGCCAGCTGTGCCGGGTCGTTCATTCATTTTTTTGTGGCACTGCTGCGCGGCCAGTGCCTGGAGGGCCAGGCGTTGTACGTGGACGCGTCGCGAACCTATCa                                                                                                       <  2:1084727/147‑1 (MQ=255)                                        ccGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCa                                                                                                       <  2:928024/150‑1 (MQ=255)                                            ggcgATTGCGCACGGGTCATGCTCGACCTGTCGGAGCTACACCGCTCCAGCCGTGCTCGGCCGTTCAGCCGCCCCGAAGTGTCACTCCTGCGCGG‑CAGCGCCTGGAGCGCCCGGCGTTGTACGTGTACGCGGCGCAAACCTATCAGTAcc                                                                                                  <  2:1104820/150‑1 (MQ=255)                                                aTTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTa                                                                                               <  2:842034/150‑1 (MQ=255)                                                aTTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTa                                                                                               >  1:842034/1‑150 (MQ=255)                                                   gcgcACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATcc                                                                                            <  2:1787346/150‑1 (MQ=255)                                                    cgcACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGGCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCCGGCGTTGTACGTGGACCCAGAGCACACGTGACAGTGCCTTATAca                                                                                           >  2:178598/1‑150 (MQ=255)                                                      cACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCcagc                                                                                         >  2:203243/1‑150 (MQ=255)                                                       aCGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCcagca                                                                                        <  2:187814/150‑1 (MQ=255)                                                              aTGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCg                                                                                 <  2:299613/150‑1 (MQ=255)                                                               tGCTCGCACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCGGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGCCTCGGCGCAAACCTATCAGCACCTTATCCAGCAGTACCCgg                                                                                <  1:1750196/150‑1 (MQ=255)                                                                  tCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGt                                                                                                                            <  1:1666039/103‑1 (MQ=255)                                                                  tCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGt                                                                                                                            >  2:1666039/1‑103 (MQ=255)                                                                  tCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGa                                                                                                                         <  2:1667564/106‑1 (MQ=255)                                                                  tCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGa                                                                                                                         >  1:1667564/1‑106 (MQ=255)                                                                  tCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCa                                                                             >  2:191686/1‑150 (MQ=255)                                                                   cGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCaa                                                                            >  2:2250171/1‑150 (MQ=255)                                                                     aaCTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACg                                                                          >  1:806936/1‑150 (MQ=255)                                                                        tGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAAGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAat                                                                       <  2:2170896/150‑1 (MQ=255)                                                                          tCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCa                                                                                                                                        >  1:1347107/1‑83 (MQ=255)                                                                          tCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCa                                                                                                                                        <  2:1347107/83‑1 (MQ=255)                                                                          tCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGt                                                                                                                            <  2:233481/95‑1 (MQ=255)                                                                          tCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGt                                                                                                                            >  1:233481/1‑95 (MQ=255)                                                                          tCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACCTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAatat                                                                     <  1:1489088/150‑1 (MQ=255)                                                                            gCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                   <  2:1571775/150‑1 (MQ=255)                                                                                gTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTAt                                                               >  1:1163510/1‑150 (MQ=255)                                                                                 ttGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATc                                                              <  2:2067176/150‑1 (MQ=255)                                                                                     tCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCgcgc                                                          >  2:1817939/1‑150 (MQ=255)                                                                                         aCCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATAAGTCAATTAGTGATCAGTACCCGGGCAACGAATATGCCTATCGCGCCCgg                                                      >  1:1740258/1‑150 (MQ=255)                                                                                           cAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGTGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                   <  2:1951291/135‑1 (MQ=255)                                                                                           cAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGTGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                   >  1:1951291/1‑135 (MQ=255)                                                                                             gCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAATGAATATGCCTATCGCGCCCGGgcgc                                                  <  1:2155703/150‑1 (MQ=255)                                                                                               cGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgc                                                <  1:201415/150‑1 (MQ=255)                                                                                                 tGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgata                                              <  1:1587441/150‑4 (MQ=255)                                                                                                    ccGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCa                                           <  1:191686/150‑1 (MQ=255)                                                                                                      ggCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGa                                         <  2:1163510/150‑1 (MQ=255)                                                                                                         cGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGAccc                                      <  1:203243/150‑1 (MQ=255)                                                                                                         cGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGAccc                                      >  1:2448370/1‑150 (MQ=255)                                                                                                           ttCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTg                                    >  1:1939544/1‑150 (MQ=255)                                                                                                              aTCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGaa                                 <  2:620652/150‑1 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGTAAACCTATCAGTACCTTATc                                                                                             <  2:1422612/89‑1 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATc                                                                                             >  1:1422612/1‑89 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                   >  1:1056184/1‑115 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                   <  2:1056184/115‑1 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgatag                                             <  2:2385520/137‑5 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgatag                                             <  2:2292152/137‑5 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgatag                                             >  1:2292152/1‑133 (MQ=255)                                                                                                               tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgatag                                             >  1:2385520/1‑133 (MQ=255)                                                                                                                ccACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGaaaa                               >  2:1625264/1‑150 (MQ=255)                                                                                                                  aCCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAAc                             <  1:1623007/150‑1 (MQ=255)                                                                                                                    cccGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTc                           <  2:1836335/150‑1 (MQ=255)                                                                                                                     ccGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCAATGACGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCg                          <  1:2409634/150‑1 (MQ=255)                                                                                                                     ccGAAGGGTCACTGCTGCGCGGCCAGTGCCTGGGGCGCCAGGCGTTGCACGTGGACGCGGCGGAATCCTATTAGTACCCTATCCAGCCGTACCCGGGCAACGAATACGCCTAGCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAACCTCg                          <  1:2250171/150‑1 (MQ=255)                                                                                                                          gtgtCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTAt                                                                                                         <  2:933411/66‑1 (MQ=255)                                                                                                                          gtgtCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTAt                                                                                                         >  1:933411/1‑66 (MQ=255)                                                                                                                            gtCACTGCTGCGCGGCCAGTGCCTGGAGGGCTAGGCGTTGTACGTGGACGCGGCGAAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCAttt                   <  2:1075550/150‑1 (MQ=255)                                                                                                                            gtCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCAttt                   <  1:1217628/150‑1 (MQ=255)                                                                                                                               aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTg                                    <  2:1996047/130‑1 (MQ=255)                                                                                                                               aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTg                                    >  1:1996047/1‑130 (MQ=255)                                                                                                                               aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCg                          >  1:736108/1‑140 (MQ=255)                                                                                                                               aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCg                          <  2:736108/140‑1 (MQ=255)                                                                                                                               aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTgcg                >  1:373072/1‑150 (MQ=255)                                                                                                                                 tgctgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAat                                                                       <  1:1709732/93‑1 (MQ=255)                                                                                                                                 tgctgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAat                                                                       >  2:1709732/1‑93 (MQ=255)                                                                                                                                    tgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTg           >  2:1802805/1‑150 (MQ=255)                                                                                                                                           ccAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTCCCATATCCGGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAAc    <  2:2157391/150‑1 (MQ=255)                                                                                                                                             aacgCCTGGCGCTCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACgg  >  2:1959755/4‑150 (MQ=255)                                                                                                                                             aacgCCTGGCGCTCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACgg  >  2:177779/4‑150 (MQ=255)                                                                                                                                              acgCCTGGCGCTCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACGTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTc                           >  1:1357631/3‑124 (MQ=255)                                                                                                                                              acgCCTGGCGCTCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTc                           <  2:1357631/122‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             CATGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGG  >  NC_002947/1553103‑1553392

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