Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,125,370	G→C	7.1%	G55R (GGT→CGT) ‡	PP_4512 →	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,125,370	0	G	C	7.1%	212.2 / 4.3	66	G55R (GGT→CGT) ‡	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base C (3/2);  total (37/29)
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

GGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGG  >  NC_002947/5125229‑5125517                                                                                                                                                |                                                                                                                                                     ggCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCgg                                                                                                                                                <  2:213588/148‑1 (MQ=255) ggCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCg                                                                                                                                                 <  1:118786/147‑1 (MQ=255)      aCTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGGCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCGGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGc                                                                                                                                          <  5:240/149‑1 (MQ=255)       cTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGc                                                                                                                                          >  5:100529/1‑148 (MQ=255)                ggCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGa                                                                                                                                 <  7:14879/148‑1 (MQ=255)                  cgcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagca                                                                                                                              <  7:65328/149‑1 (MQ=255)                   gcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagcag                                                                                                                             >  5:224024/1‑149 (MQ=255)                    cgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagcagc                                                                                                                            >  5:210735/1‑149 (MQ=255)                       tGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAg                                                                                                                         >  5:268937/1‑149 (MQ=255)                                gCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTCGGAGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgcgcg                                                                                                                 >  1:101771/1‑149 (MQ=255)                                gCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgcgcgc                                                                                                                >  3:63613/1‑149 (MQ=255)                                 ccTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgcc                                                                                                               >  3:129655/1‑149 (MQ=255)                                   ttACCGGCGTCTTCGGCAAACTGGACACCAGCGCCACCCAGGCGATCATAGTC‑TTCGGTA‑GCGCCAACATCGCCGTACTCCCGCTGGCCCTGTTCGCCAGCCTGGCCCGTCAGGGCTGGCAACGCCTGAGCAGCCAGGAAGCGcgccgc                                                                                                             >  4:109693/1‑149 (MQ=255)                                     aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccgc                                                                                                             <  8:17674/147‑1 (MQ=255)                                     aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCct                                                                                                           >  7:52296/1‑149 (MQ=255)                                     aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGc                                                                                                                                          <  8:289768/118‑1 (MQ=255)                                     aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGc                                                                                                                                          >  7:289768/1‑118 (MQ=255)                                       cGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagcagc                                                                                                                            >  4:151670/1‑130 (MQ=255)                                       cGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagcagc                                                                                                                            <  3:151670/130‑1 (MQ=255)                                         gCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCaa                                                                                                                                        <  4:198667/116‑1 (MQ=255)                                         gCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCct‑gctg                                                                                                      <  6:224024/149‑1 (MQ=255)                                           gTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑gctgct                                                                                                    <  4:129655/149‑1 (MQ=255)                                               tCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggt                                                                                                >  1:322840/1‑149 (MQ=255)                                                cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAG‑CAGGCGGCGCGCATCCTGGCTGCTca                                                                                                  <  6:100529/146‑2 (MQ=255)                                                 ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtgg                                                                                              >  6:249899/1‑149 (MQ=255)                                                 ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtgg                                                                                              >  6:257107/1‑149 (MQ=255)                                                 ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtgg                                                                                              >  7:84258/1‑149 (MQ=255)                                                 ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtgg                                                                                              <  3:101428/149‑1 (MQ=255)                                                     aaCTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGtgtg                                                                                          >  7:82935/1‑149 (MQ=255)                                                     aaCTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGtgtg                                                                                          >  3:12785/1‑149 (MQ=255)                                                      aCTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGtgc                                                                                         <  5:257107/149‑1 (MQ=255)                                                        tGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccgc                                                                                                             >  1:165561/1‑128 (MQ=255)                                                        tGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccgc                                                                                                             <  2:165561/128‑1 (MQ=255)                                                        tGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGtgctg                                                                                       >  6:53680/1‑149 (MQ=255)                                                          gCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTgctgct                                                                                     <  8:52296/149‑1 (MQ=255)                                                           ccaccaGCGCCAGCCCGGCGATCATCGTCTTTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTgctgct                                                                                     >  1:36598/1‑149 (MQ=255)                                                                aGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggccgg                                                                               <  6:210735/149‑1 (MQ=255)                                                                    ccAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCAc                                                                           >  5:258716/1‑149 (MQ=255)                                                                     cAGCCCGGCGATCATCGGC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACt                                                                          <  4:12785/149‑1 (MQ=255)                                                                      aGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTg                                                                         <  3:109693/149‑1 (MQ=255)                                                                       gCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTgg                                                                        <  6:277314/149‑1 (MQ=255)                                                                         ccGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGt                                                                      <  8:84258/149‑1 (MQ=255)                                                                                tcatcGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttc                                                               <  4:63613/149‑1 (MQ=255)                                                                                  atcGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGtgctg                                                                                       >  1:109792/1‑123 (MQ=255)                                                                                  atcGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGtgctg                                                                                       <  2:109792/123‑1 (MQ=255)                                                                                  atcGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttctt                                                             >  5:15134/1‑149 (MQ=255)                                                                                   tcGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttctt                                                             <  1:18783/148‑1 (MQ=255)                                                                                   tcGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttctt                                                             >  2:18783/1‑148 (MQ=255)                                                                                     gTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCAt                                                          >  2:64689/1‑149 (MQ=255)                                                                                      tC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATc                                                         >  1:338702/1‑149 (MQ=255)                                                                                                 gcgcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCaaa                                                <  5:249899/149‑1 (MQ=255)                                                                                                 gcgcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCaaa                                                <  5:174180/149‑1 (MQ=255)                                                                                                       ccTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGggc                                          <  2:338702/149‑1 (MQ=255)                                                                                                         ttCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGggcgg                                        >  7:267971/1‑149 (MQ=255)                                                                                                            gCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGt                                     <  2:101771/149‑1 (MQ=255)                                                                                                              cGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAg                                   >  3:261418/1‑149 (MQ=255)                                                                                                              cGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAg                                   >  1:295023/1‑149 (MQ=255)                                                                                                                   tCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATc                              >  8:80803/1‑149 (MQ=255)                                                                                                                         tGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCAcc                        >  5:280686/1‑149 (MQ=255)                                                                                                                            cccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTc                     >  1:250190/1‑149 (MQ=255)                                                                                                                                 gttgCCAGCCCGGACCGGCCAGGCTGGCGACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCtt                >  7:4100/4‑149 (MQ=255)                                                                                                                                      cAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCa           <  5:53680/149‑1 (MQ=255)                                                                                                                                         ccTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCt        >  7:279311/1‑149 (MQ=255)                                                                                                                                           tGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTc      >  3:123125/1‑149 (MQ=255)                                                                                                                                           tGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTc      >  1:118538/1‑149 (MQ=255)                                                                                                                                               cGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCgg  >  7:266547/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     GGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGG  >  NC_002947/5125229‑5125517

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