Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,227,379	C→G	8.6%	A73G (GCC→GGC)	PP_2828 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,227,379	0	C	G	8.6%	206.8 / 9.2	70	A73G (GCC→GGC)	PP_2828	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base G (2/4);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CCAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAGGATGACCCGCG  >  NC_002947/3227236‑3227520                                                                                                                                                |                                                                                                                                               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                                                                                                                           <  8:277651/148‑1 (MQ=255)                  ccTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           <  2:64312/148‑1 (MQ=255)                    ttGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           <  5:194837/146‑1 (MQ=255)                    ttGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCt                                                                                                                           >  6:194837/1‑146 (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                                                                                                         <  2:274612/147‑1 (MQ=255)                                    aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACtt                                                                                                         >  1:274612/1‑147 (MQ=255)                                    aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAc                                                                                                                                        >  1:163889/1‑117 (MQ=255)                                    aaCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAc                                                                                                                                        <  2:163889/117‑1 (MQ=255)                                     aCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAcatg                                                                                                                                     <  8:299557/119‑4 (MQ=255)                                     aCGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGCGGCACGCAAcatg                                                                                                                                     >  7:299557/1‑116 (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                                                                          >  4:297400/1‑149 (MQ=255)                                                                 gAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGcc                                                                          <  3:116350/149‑1 (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                                <  4:100924/127‑1 (MQ=255)                                                                                                                                 tgGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGttt                                >  3:100924/1‑127 (MQ=255)                                                                                                                                   gCCATGTTGCGTGACGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCCGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGGTTCAGTGCGTTGTTCATGCAGGATGa        >  4:174195/1‑149 (MQ=255)                                                                                                                                         ttGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAGGATGACCcgcg  >  2:255204/1‑149 (MQ=255)                                                                                                                                            cGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAgg            >  2:356353/1‑136 (MQ=255)                                                                                                                                            cGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAgg            <  1:356353/136‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               CCAGGCGCAGCGCGCTTCCTTGCAGTTGCTGATACAACGCCTGTTGGTCGCTGCAGGTGGGCCGGAGTACATCGGTACTTTCAGCTTGCAGGTATTGCAAGGGAGCGACCGCCATAGCAGCCGGCTGCTGGCCATGTTGCGTGCCGCGCAACTGAGCATTGCCCTG‑CGTTCGCCGGTAACTTTCCAGCTGCGGGTATTGGTGGTCAGCCTGCCGGTTGCCGATAACGCCACGCTGGAAGCTCACGAGCGCAGTTTCAGTGCGTTGTTCATGCAGGATGACCCGCG  >  NC_002947/3227236‑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: ‑