Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,503,212	G→T	5.9%	G281C (GGT→TGT)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,503,212	0	G	T	5.9%	210.5 / 3.3	68	G281C (GGT→TGT)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base T (2/2);  total (34/34)
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.35e-01

TACCTGGGCCGCCTGGACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGATTTTTGC  >  NC_002947/4503068‑4503357                                                                                                                                                 |                                                                                                                                                   tACCTGGGCCGCCTGGACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGa                                                                                                                                                <  3:105977/149‑1 (MQ=255)    cTGGGCCGCCTGGACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGt                                                                                                                                             >  6:384278/1‑149 (MQ=255)     tGGGCCGCCTGGACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGt                                                                                                                                             >  4:176419/1‑148 (MQ=255)              tGGACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGccac                                                                                                                                    >  6:136377/1‑148 (MQ=255)                gACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGccacca                                                                                                                                  <  3:129272/148‑1 (MQ=255)                     tACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAg                                                                                                                             <  7:147324/148‑1 (MQ=255)                      aCCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGgt                                                                                                                           >  1:129297/1‑149 (MQ=255)                        cGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGTCATCGGTGACGTTACCGCCACCAGCAAGgtgt                                                                                                                         >  1:160831/1‑149 (MQ=255)                           cATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTcgc                                                                                                                      <  8:326381/149‑1 (MQ=255)                           cATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTcgc                                                                                                                      <  7:176147/149‑1 (MQ=255)                                aCGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAg                                                                                                                 >  2:104210/1‑149 (MQ=255)                                     gCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGt                                                                                                            <  3:176419/149‑1 (MQ=255)                                     gCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGt                                                                                                            >  3:128577/1‑149 (MQ=255)                                         cgcTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGTCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGaa                                                                                                        <  2:160831/149‑1 (MQ=255)                                         cgcTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGaa                                                                                                        <  1:163641/149‑1 (MQ=255)                                              ttCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAggg                                                                                                     <  1:362636/147‑1 (MQ=255)                                                 cctgTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGa                                                                                                 <  4:249494/145‑1 (MQ=255)                                                 cAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGcc                                                                                                                   >  6:263106/1‑130 (MQ=255)                                                 cAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGcc                                                                                                                   <  5:263106/130‑1 (MQ=255)                                                  aGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGAtg                                                                                               >  7:14844/1‑149 (MQ=255)                                                     tGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTTGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCTAAGGGCCGATGTgg                                                                                            <  4:304197/149‑1 (MQ=255)                                                       ccGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGAt                                                                                          <  5:136377/149‑1 (MQ=255)                                                       ccGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGACGGTCGAAGGGCCGCTGTGGAt                                                                                          >  3:315237/1‑149 (MQ=255)                                                         gACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATg                                                                                         >  1:275239/1‑148 (MQ=255)                                                              cGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCTGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCg                                                                                   >  6:133967/1‑148 (MQ=255)                                                              cGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCTGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCg                                                                                   >  2:62398/1‑148 (MQ=255)                                                              cGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑Gc                                                                                    >  8:342167/1‑147 (MQ=255)                                                                cTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGAt                                                                                <  1:96395/149‑1 (MQ=255)                                                                   cACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCTGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCaa                                                                             <  1:62398/149‑1 (MQ=255)                                                                   cACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCTGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCaa                                                                             <  5:133967/149‑1 (MQ=255)                                                                     cTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAgg                                                                           >  7:117618/1‑149 (MQ=255)                                                                        tGGCCTGGGCTTCTCGTGGGCGCTGGTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAggg                                                                          <  7:256590/147‑1 (MQ=255)                                                                         ggCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAggggg                                                                        <  2:275239/148‑1 (MQ=255)                                                                           ccTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCa                                                                              <  8:76881/140‑1 (MQ=255)                                                                           ccTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCa                                                                              >  7:76881/1‑140 (MQ=255)                                                                                gCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTg                                                                <  2:15754/149‑1 (MQ=255)                                                                                gCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTg                                                                <  2:129297/149‑1 (MQ=255)                                                                                  ttCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGc                                                                      <  6:332442/141‑1 (MQ=255)                                                                                  ttCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGc                                                                      >  5:332442/1‑141 (MQ=255)                                                                                   tctcGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCACGGGGGGCGTGCGGGt                                                              >  8:255508/1‑148 (MQ=255)                                                                                    ctcGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTc                                                             >  4:70504/1‑148 (MQ=255)                                                                                     tcGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGAt                                                                                                <  4:238567/113‑1 (MQ=255)                                                                                       gTGGGAGCTTTTCATTCAGATGCTTTTGATCTACCTGTTCACTTCTCTGGAAGCCATAGGTGACGTTACCGCCACCAGCAAGGTGACGCGCCAGCCGGTCGAAGGGCCGATGAGGATGCA‑GCAGATCAAGCTGGGCGTGCTGGTAAACg                                                         >  4:105194/1‑149 (MQ=255)                                                                                        tccgatcTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAAc                                                          <  3:109627/141‑1 (MQ=255)                                                                                               tGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCaa                                                  >  8:128321/1‑148 (MQ=255)                                                                                                 ttCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGGGTGATGGTCatg                                                          >  4:109627/1‑136 (MQ=255)                                                                                                 ttCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAAc                                                 <  8:117618/147‑1 (MQ=255)                                                                                                 ttCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTc                                               <  3:105194/149‑1 (MQ=255)                                                                                                    aTTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGAt                                                                                >  8:333171/1‑113 (MQ=255)                                                                                                    aTTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGAt                                                                                <  7:333171/113‑1 (MQ=255)                                                                                                    aTTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCgct                                            >  8:298387/1‑149 (MQ=255)                                                                                                    aTTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCgct                                            >  5:138348/1‑149 (MQ=255)                                                                                                    aTTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCgct                                            >  2:292512/1‑149 (MQ=255)                                                                                                       ccGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCgctgc                                          <  7:128321/148‑1 (MQ=255)                                                                                                               ggTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTc                                 >  4:344906/1‑149 (MQ=255)                                                                                                                        ccTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCATCCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGtt                             <  1:104210/145‑1 (MQ=255)                                                                                                                          tGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGa                                                                                           <  3:97985/81‑1 (MQ=255)                                                                                                                          tGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGa                                                                                           >  4:97985/1‑81 (MQ=255)                                                                                                                          tGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGCTCGAAGGGCCGATGTGCATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGTTCTGTTCAACAcc                      >  5:209714/1‑149 (MQ=255)                                                                                                                                 ttCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTccc                 >  7:226519/1‑147 (MQ=255)                                                                                                                                   cGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTa               <  5:165391/147‑1 (MQ=255)                                                                                                                                    gCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCt            >  1:997/1‑149 (MQ=255)                                                                                                                                            ccATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGGCGGTCGTAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCTACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGAttttt    <  5:384278/149‑1 (MQ=255)                                                                                                                                            ccATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGAttttt    <  3:70504/149‑1 (MQ=255)                                                                                                                                             cATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGATTTTTg   >  3:172134/1‑149 (MQ=255)                                                                                                                                              aTCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGTGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGATTTTTGc  <  7:298387/149‑1 (MQ=255)                                                                                                                                              aTCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGATTTTTGc  >  3:365585/1‑149 (MQ=255)                                                                                                                                              aTCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGATTTTTGc  >  4:172354/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TACCTGGGCCGCCTGGACTTTACCGGCATGCACGAGGCCGCGCTGTTCCAGGTGCCGACACCGCTGCACTTTGGCCTGGGCTTCTCGTGGGCGCTGTTCATTCCGATGCTGGTGATCTACCTGGTCACTTCGCTGGAAGCCATCGGTGACGTTACCGCCACCAGCAAGGTGTCGCGCCAGCCGGTCGAAGGGCCGATGTGGATGCA‑GCGGATCAAGGGGGGCGTGCTGGTCAACGGCGCCAACTCGCTGCTGGCCGGTCTGTTCAACACCTTCCCTAGCTCGATTTTTGC  >  NC_002947/4503068‑4503357

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