Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,227,377	T→C	8.7%	R72R (CGT→CGC)	PP_2828 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,227,377	0	T	C	8.7%	192.1 / 7.2	70	R72R (CGT→CGC)	PP_2828	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/4);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

ACATGTCCCAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAGGATGACCCGCG  >  NC_002947/3227229‑3227520                                                                                                                                                     |                                                                                                                                                 acaTGTCCCAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGt                                                                                                                                                  <  3:124776/149‑1 (MQ=255)        ccAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGTCCATGTTGCGTGCcgcg                                                                                                                                            >  8:116996/1‑148 (MQ=255)        ccAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCcgcg                                                                                                                                            >  4:161292/1‑148 (MQ=255)        ccAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCcgcg                                                                                                                                            >  7:303295/1‑148 (MQ=255)        ccAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCgcgc                                                                                                                                           <  5:291710/149‑1 (MQ=255)         cAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCgcgc                                                                                                                                           >  7:277651/1‑148 (MQ=255)              gcAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAcatg                                                                                                                                     <  8:250631/149‑4 (MQ=255)              gcAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAcatg                                                                                                                                     <  7:178259/149‑4 (MQ=255)                 gcgcgcTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCa                                                                                                                                  >  7:295889/1‑149 (MQ=255)                   gcgcTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCAtt                                                                                                                                <  8:295889/149‑1 (MQ=255)                         ccTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGTGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           <  7:282006/148‑1 (MQ=255)                         ccTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGccc                                                                                                                            <  7:102528/147‑1 (MQ=255)                         ccTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           <  2:64312/148‑1 (MQ=255)                         ccTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           <  8:277651/148‑1 (MQ=255)                           ttGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           >  6:194837/1‑146 (MQ=255)                           ttGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           <  5:194837/146‑1 (MQ=255)                           ttGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑Cg                                                                                                                       >  1:23531/1‑149 (MQ=255)                            tGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGTAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGt                                                                                                                      >  7:94879/1‑149 (MQ=255)                             gCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑Cg                                                                                                                       <  3:161292/147‑1 (MQ=255)                              cAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTc                                                                                                                    <  3:155960/149‑1 (MQ=255)                                   gCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCgg                                                                                                               >  6:342554/1‑149 (MQ=255)                                      gATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGtt                                                                                                                     <  6:274861/140‑1 (MQ=255)                                      gATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGtt                                                                                                                     >  5:274861/1‑140 (MQ=255)                                         aCAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACtt                                                                                                         >  6:29382/1‑149 (MQ=255)                                           aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTTCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTc                                                                                                       <  1:256998/149‑1 (MQ=255)                                           aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACtt                                                                                                         >  1:274612/1‑147 (MQ=255)                                           aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACtt                                                                                                         <  2:274612/147‑1 (MQ=255)                                           aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAc                                                                                                                                        >  1:163889/1‑117 (MQ=255)                                           aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAc                                                                                                                                        <  2:163889/117‑1 (MQ=255)                                            aCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAcatg                                                                                                                                     >  7:299557/1‑116 (MQ=255)                                            aCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAcatg                                                                                                                                     <  8:299557/119‑4 (MQ=255)                                              gCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAg                                                                                                    >  3:355224/1‑149 (MQ=255)                                                cTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCACGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTTGCCATGTTTCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCt                                                                                                  >  1:68756/1‑149 (MQ=255)                                                 tGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGAATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATAGCCCTC‑CGTTCGCCGGTAACTTTCCAGCTg                                                                                                 >  3:254539/1‑149 (MQ=255)                                                  gTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTg                                                                                                                          <  6:265177/124‑1 (MQ=255)                                                  gTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTg                                                                                                                          >  5:265177/1‑124 (MQ=255)                                                     ggTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCgg                                                                                              <  5:29382/148‑1 (MQ=255)                                                      gTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCggg                                                                                             >  4:49570/1‑148 (MQ=255)                                                      gTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCggg                                                                                             >  2:89113/1‑148 (MQ=255)                                                             cAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATtggtg                                                                                     >  4:61436/1‑149 (MQ=255)                                                               ggTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTggtggt                                                                                   <  5:342554/149‑1 (MQ=255)                                                                        gAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGcc                                                                          <  3:116350/149‑1 (MQ=255)                                                                        gAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGcc                                                                          >  4:297400/1‑149 (MQ=255)                                                                         aGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCg                                                                         <  3:61436/149‑1 (MQ=255)                                                                             cATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTGCCGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGtt                                                                      <  1:89113/149‑1 (MQ=255)                                                                             cATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGtt                                                                      <  3:49570/148‑1 (MQ=255)                                                                              atctggaCTTTCAGCTTGCAGGTTTTGCAAGGGAGCGCCCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGGATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCGGCGGGTATTGGTGGGCAGCCTGCCGGTTGc                                                                    <  4:254539/143‑1 (MQ=255)                                                                              aTCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCgg                                                                                              <  2:185997/123‑1 (MQ=255)                                                                              aTCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCgg                                                                                              >  1:185997/1‑123 (MQ=255)                                                                              aTCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGc                                                                    >  2:79572/1‑149 (MQ=255)                                                                                 ggTACTTTCAGCTTGCAGGTATTGTAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCg                                                                  <  8:94879/148‑1 (MQ=255)                                                                                  gTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCTGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGAt                                                                <  2:23531/149‑1 (MQ=255)                                                                                  gTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGa                                                                 >  1:108413/1‑148 (MQ=255)                                                                                     cTTTAAGCTTGAAGGTATTGCAAGGGAGAGACCTCCATAGCAGCCGGCTTCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTTCCGGTTGCCGATaa                                                              <  7:116996/148‑1 (MQ=255)                                                                                      tttCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGcc                                                                                                                             <  4:317572/85‑1 (MQ=255)                                                                                      tttCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGcc                                                                                                                             >  3:317572/1‑85 (MQ=255)                                                                                      tttCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACg                                                            >  5:210524/1‑149 (MQ=255)                                                                                                   gTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATtg                                                                                        >  2:274717/1‑108 (MQ=255)                                                                                                   gTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATtg                                                                                        <  1:274717/108‑1 (MQ=255)                                                                                                      ttGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCa                                            <  1:265730/149‑1 (MQ=255)                                                                                                         cAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGa                                         <  1:79572/149‑1 (MQ=255)                                                                                                           aGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAgc                                       <  6:210524/149‑1 (MQ=255)                                                                                                                        cATAGCCGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGc                          <  8:252515/149‑1 (MQ=255)                                                                                                                        cATAGCAGCCGGCTGCTGGCCCTGTTTCGGGCCGCGCAACTGAGCATTGCCCTG‑CGTGCGCCGGTAACTTTCCTGCTGCGGGTTTGGGTGGTCAGCCTGCCGGTTGCCGATCACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGc                          <  1:286478/149‑1 (MQ=255)                                                                                                                                       ctgGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAgg            <  3:297400/148‑1 (MQ=255)                                                                                                                                        tgGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGttt                                >  3:100924/1‑127 (MQ=255)                                                                                                                                        tgGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGttt                                <  4:100924/127‑1 (MQ=255)                                                                                                                                          gCCATGTTGCGTGACGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCCGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGGTTCAGTGCGTTGTTCATGCAGGATGa        >  4:174195/1‑149 (MQ=255)                                                                                                                                                ttGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAGGATGACCcgcg  >  2:255204/1‑149 (MQ=255)                                                                                                                                                   cGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAgg            <  1:356353/136‑1 (MQ=255)                                                                                                                                                   cGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAgg            >  2:356353/1‑136 (MQ=255)                                                                                                                                                     |                                                                                                                                                 ACATGTCCCAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAGGATGACCCGCG  >  NC_002947/3227229‑3227520

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