Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,125,368	C→A	6.1%	A54D (GCC→GAC)	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,368	0	C	A	6.1%	211.5 / 3.2	65	A54D (GCC→GAC)	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base C (34/27);  new base A (2/2);  total (36/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 = 8.49e-01

ACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTC  >  NC_002947/5125220‑5125513                                                                                                                                                       |                                                                                                                                                   accgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGc                                                                                                                                                    <  2:156294/149‑1 (MQ=255)          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                                                                                                                        <  3:151670/130‑1 (MQ=255)                                                cGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagcagc                                                                                                                        >  4:151670/1‑130 (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                                                                                          >  7:84258/1‑149 (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                                                                                          <  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                                                                                                         <  2:165561/128‑1 (MQ=255)                                                                 tGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccgc                                                                                                         >  1:165561/1‑128 (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  >  1:118538/1‑149 (MQ=255)                                                                                                                                                    tGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTc  >  3:123125/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                   ACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC‑TTCGGTC‑GCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTC  >  NC_002947/5125220‑5125513

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