Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,312,198	G→A	6.1%	P820L (CCA→CTA)	rapA ←	RNA polymerase‑associated protein RapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,312,198	0	G	A	6.1%	193.6 / 2.2	65	P820L (CCA→CTA)	rapA	RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTGGCCAGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGCCTGACAGCAC  >  NC_002947/1312052‑1312340                                                                                                                                                     |                                                                                                                                               gCTTGGCCAGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTggc                                                                                                                                              <  5:109318/149‑1 (MQ=255) gCTTGGCCAGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTggc                                                                                                                                              >  7:12947/1‑149 (MQ=255)     ggCCAGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTggcggc                                                                                                                                           <  7:125343/148‑1 (MQ=255)      gCCAGCGCATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTGGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCa                                                                                                                                          <  7:288529/148‑1 (MQ=255)       ccAGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAgg                                                                                                                                        <  2:260462/149‑1 (MQ=255)         aGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTa                                                                                                                                      <  4:101131/149‑1 (MQ=255)            acaTCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCg                                                                                                                                   <  6:298645/149‑1 (MQ=255)             caTCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCgg                                                                                                                                  >  1:56098/1‑149 (MQ=255)              aTCGCGCTGGGCCTGGACGAACTTGTTGGAGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGc                                                                                                                                 >  5:22617/1‑149 (MQ=255)                 gcgcTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGccc                                                                                                                               <  8:123116/148‑1 (MQ=255)                         ccTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAgg                                                                                                                      <  8:135596/149‑1 (MQ=255)                            ggACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTa                                                                                                                   >  6:3790/1‑149 (MQ=255)                                          gcgctggcgcGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCCTTGGCGTCGAGCAGGCA‑GCGCCGTGCCGCTGGCGGCAGGTAGCGGCCCTGTTGCAGGCTACGGGGTGCccccac                                                                                                     >  6:191185/1‑149 (MQ=255)                                            gctggcgcGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCt                                                                                                   >  8:105282/1‑149 (MQ=255)                                               ggcgcGCGGTACACTCTCGAGCTGGGCGTTGCGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCATTTGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCg                                                                                                 <  1:92857/149‑1 (MQ=255)                                                gcgcgcGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGCAAACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGATGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGc                                                                                                >  1:67679/1‑149 (MQ=255)                                                gcgcgcGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGc                                                                                                <  2:67679/148‑1 (MQ=255)                                                          cacTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAAcag                                                                                      <  5:3790/148‑1 (MQ=255)                                                          cacTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAAcag                                                                                      <  8:12947/148‑1 (MQ=255)                                                           acTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAAcagca                                                                                    <  7:105282/149‑1 (MQ=255)                                                             tctcGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAAca                                                                                       >  4:282224/1‑144 (MQ=255)                                                             tctcGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAAca                                                                                       <  3:282224/144‑1 (MQ=255)                                                             tctcGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGt                                                                                  >  7:66542/1‑149 (MQ=255)                                                               tcGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGtt                                                                                 >  2:207310/1‑148 (MQ=255)                                                                 gAGCTGGTCGTTGAGGGGTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCa                                                                               <  5:11979/148‑1 (MQ=255)                                                                   gCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGc                                                                            <  8:226045/149‑1 (MQ=255)                                                                   gCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGc                                                                            <  5:12933/149‑1 (MQ=255)                                                                       ggcGTTGAGGGTTTCGAAAGC‑AACGAGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTAGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGTGGTGCCACCACCTCGATGACGAACAGCAGTTCAAGCAGGa                                                                        >  7:196169/3‑149 (MQ=255)                                                                       gTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGa                                                                        <  5:317230/149‑1 (MQ=255)                                                                       gTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTAGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGa                                                                        >  7:191168/1‑149 (MQ=255)                                                                         cGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACa                                                                      <  3:288477/149‑1 (MQ=255)                                                                               gggTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCt                                                                >  6:123757/1‑149 (MQ=255)                                                                                      gAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGttt                                                           >  2:142818/1‑147 (MQ=255)                                                                                      gAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGttt                                                           >  1:84669/1‑147 (MQ=255)                                                                                       aaaGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGc                                                                            >  7:40718/1‑129 (MQ=255)                                                                                       aaaGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGc                                                                            <  8:40718/129‑1 (MQ=255)                                                                                       aaaGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAg                                                        >  1:171088/1‑149 (MQ=255)                                                                                        aaGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAgc                                                       <  2:84669/149‑1 (MQ=255)                                                                                          gC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGa                                                         >  8:113745/1‑145 (MQ=255)                                                                                          gC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGa                                                         <  7:113745/145‑1 (MQ=255)                                                                                           c‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCt                                                    >  5:177668/1‑149 (MQ=255)                                                                                           c‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCt                                                    >  8:306366/1‑149 (MQ=255)                                                                                             aaCGCGTGACGCCAGGTCGTTGCCATTGGCGTCGCGCAGGCA‑GCGCAGTACCGCTGGCGGCAGGGAGCGGCCCAGTTGCAGGCTACGGGGGGCCACAACCTCGCTGAAGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCtt                                                   >  3:219970/1‑149 (MQ=255)                                                                                              aCGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCa                                                                           <  6:278858/124‑1 (MQ=255)                                                                                              aCGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCttt                                                  <  2:56098/149‑1 (MQ=255)                                                                                                gcgTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGt                                                >  4:107013/1‑149 (MQ=255)                                                                                                     aCGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCtt                                            <  8:66542/148‑1 (MQ=255)                                                                                                        ccAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATc                                        <  1:207310/149‑1 (MQ=255)                                                                                                               gTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCa                                 >  5:91300/1‑149 (MQ=255)                                                                                                               gTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGAGCTTTGTTCTTGATCAACGCCa                                 >  3:57782/1‑149 (MQ=255)                                                                                                                 tGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCAcg                               >  1:21655/1‑149 (MQ=255)                                                                                                                    cATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCAcggcg                            >  1:65547/1‑149 (MQ=255)                                                                                                                      gtGGCGTCGAGCAGGCA‑GCGGAGTGGCGCTGGCGGCTGGTAGCGGCCGAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGgt                          <  2:236578/148‑1 (MQ=255)                                                                                                                         gCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGtt                       >  8:22448/1‑149 (MQ=255)                                                                                                                           gTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGc                     <  7:306366/149‑1 (MQ=255)                                                                                                                           gTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGc                     <  6:22617/149‑1 (MQ=255)                                                                                                                             cGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGcc                    >  8:212213/1‑148 (MQ=255)                                                                                                                               aGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCAt                 <  4:57782/149‑1 (MQ=255)                                                                                                                               aGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCAt                 <  4:32997/149‑1 (MQ=255)                                                                                                                                gcaggca‑gCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATc                >  1:68866/1‑149 (MQ=255)                                                                                                                                    gca‑gCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGc            <  3:107013/149‑1 (MQ=255)                                                                                                                                    gca‑gCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGc            <  2:65547/149‑1 (MQ=255)                                                                                                                                     ca‑gCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGcc           <  4:58664/149‑1 (MQ=255)                                                                                                                                          gcAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCt                                                    >  3:215261/1‑104 (MQ=255)                                                                                                                                          gcAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCt                                                    <  4:215261/104‑1 (MQ=255)                                                                                                                                                ccGCTAGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCAGCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGCCTGACAGCAc  <  8:191168/148‑1 (MQ=255)                                                                                                                                                ccGCTAGCGGCAGGTAGCGGACCAGTTGCAGGTTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCTGTGTTGCCCATCGAGCCTGACAGCAc  <  8:196169/148‑1 (MQ=255)                                                                                                                                                  gCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGCCTGACAGCAc  >  6:124050/1‑146 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCTTGGCCAGCACATCGCGCTGGGCCTGGACGAACTTGTTGGCGCTGGCGCGCGGTACACTCTCGAGCTGGTCGTTGAGGGTTTCGAAAGC‑AACGCGTGACGCCAGGTCGTTGCCATTGGCGTCGAGCAGGCA‑GCGCAGTGCCGCTGGCGGCAGGTAGCGGCCCAGTTGCAGGCTACGGGGTGCCACCACCTCGCTGACGAACAGCAGTTCAAGCAGGACAGTGCCTGGTTTGAGCGCTTTGTTCTTGATCAACGCCACGGCGGTGTTGCCCATCGAGCCTGACAGCAC  >  NC_002947/1312052‑1312340

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