Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,311,085 C→T 6.5% A266V (GCA→GTA)  syrB → syringomycin biosynthesis protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,311,0850CT 6.5% 187.7 / 2.0 64A266V (GCA→GTA) syrBsyringomycin biosynthesis protein 2
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base T (2/2);  total (35/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.52e-01

CTTCACCATGAAGAAGGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGC  >  NC_002947/4310942‑4311233
                                                                                                                                                    |                                                                                                                                                    
cTTCACCATGAAGAAGGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCa                                                                                                                                                     <  1:306584/145‑1 (MQ=255)
cTTCACCATGAAGAAGGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAg                                                                                                                                                    <  2:272760/146‑1 (MQ=255)
 ttCACCATGAAGAAGGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGacac                                                                                                                                                >  6:256655/1‑149 (MQ=255)
  tCACCATGAAGAAGGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACAcc                                                                                                                                               >  4:282683/1‑149 (MQ=255)
            gaagGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGt                                                                                                                                     >  3:289757/1‑149 (MQ=255)
                     aTTCATCATTTTCTGGTCCACCGACCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGGGCTCCGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCCGCCACACGGACCGAGTTCGGACGaa                                                                                                                            >  6:204378/1‑149 (MQ=255)
                        catcatTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGGAGACACCCTGACCGAGTTCGGCCGAACg                                                                                                                          <  1:224803/148‑1 (MQ=255)
                         atcatTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTCCTGACGAGGAGGGTCTACCCGGATGCAGACACCCTGAACGAGTTCGGACg                                                                                                                              >  1:328879/1‑148 (MQ=255)
                         atcatTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGt                                                                                                                        >  3:91412/1‑149 (MQ=255)
                         atcatTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGt                                                                                                                        <  2:116663/149‑1 (MQ=255)
                                    gTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGaa                                                                                                              <  5:74985/148‑1 (MQ=255)
                                    gTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAAc                                                                                                             <  2:48160/149‑1 (MQ=255)
                                    gTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGAGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGGAGACACCCTGCCCGAGTTCGGCCGCACGGTGAGCCTTgca                                                                                                              >  4:316121/1‑146 (MQ=255)
                                       cACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACaa                                                                                                           <  5:256655/148‑1 (MQ=255)
                                       cACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACaa                                                                                                           >  1:237882/1‑148 (MQ=255)
                                            tcctccATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGTAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACgg                                                                                                      >  8:230075/1‑148 (MQ=255)
                                            tcctccATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACgg                                                                                                      >  3:384564/1‑148 (MQ=255)
                                                 cATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGTAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGt                                                                                                >  5:90224/1‑149 (MQ=255)
                                                 cATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCAcc                                                                                                  <  2:237882/147‑1 (MQ=255)
                                                 cATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGAACGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGt                                                                                                >  5:57516/1‑149 (MQ=255)
                                                         gCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACAGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTg                                                                                        >  5:214401/1‑149 (MQ=255)
                                                             ccTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCgg                                                                                    >  4:55658/1‑149 (MQ=255)
                                                               tCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCggg                                                                                   >  6:350205/1‑148 (MQ=255)
                                                                 aTCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGtt                                                                                         >  2:307866/1‑140 (MQ=255)
                                                                 aTCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGtt                                                                                         <  1:307866/140‑1 (MQ=255)
                                                                     gggAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACg                                                                            >  4:73290/1‑149 (MQ=255)
                                                                      ggAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGc                                                                           <  5:234325/149‑1 (MQ=255)
                                                                           gacgacCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGTAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCg                                                                      <  7:230075/149‑1 (MQ=255)
                                                                               acCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCTCCGTCATCGTTGGCGGGGAAGACGCATTCGGaca                                                                  >  4:230238/1‑149 (MQ=255)
                                                                                    gccgcGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAAcc                                                             <  4:289757/149‑1 (MQ=255)
                                                                                     ccgcGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACcg                                                            >  8:162702/1‑149 (MQ=255)
                                                                                      cgcgCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACcgc                                                           <  4:91412/149‑1 (MQ=255)
                                                                                             gggCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCggg                                                                                   <  2:136451/118‑1 (MQ=255)
                                                                                              ggCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTc                                                   >  1:389516/1‑149 (MQ=255)
                                                                                               gCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCg                                                  <  7:162702/149‑1 (MQ=255)
                                                                                                cTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCTTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGc                                                 <  3:282683/149‑1 (MQ=255)
                                                                                                 tACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGc                                                 >  3:244675/1‑148 (MQ=255)
                                                                                                    gTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAAc                                             <  6:90132/149‑1 (MQ=255)
                                                                                                         gcgcTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTc                                        >  5:12409/1‑149 (MQ=255)
                                                                                                         gcgcTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCAACGTCATCGTTGGCGGGGAAGCAGCATTCGTACACAACCGCGTGCGCTCGCAAACCACTc                                        >  4:228613/1‑149 (MQ=255)
                                                                                                         gcgcTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCAACGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTc                                        >  4:223751/1‑149 (MQ=255)
                                                                                                         gcgcTACGTGCCCACCTC‑‑‑‑CGGGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTc                                        <  5:204378/149‑1 (MQ=255)
                                                                                                            cTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTg                                     >  4:12381/1‑149 (MQ=255)
                                                                                                            cTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTAACCGGATGCAGACACCCTGCCCGAGTTCGGCCGAACGGTGAGGCTTGAACAATAAGGCACAGTCATCGTTGGAGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCTAACCACTCGTg                                     >  2:65177/1‑149 (MQ=255)
                                                                                                             tACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATc                                                                                            >  4:183134/1‑93 (MQ=255)
                                                                                                                    ctaCCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGc                                   <  8:291026/141‑1 (MQ=255)
                                                                                                                    cgaCCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGc                                   >  7:291026/3‑143 (MQ=255)
                                                                                                                     cACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAg                                >  1:252721/1‑145 (MQ=255)
                                                                                                                     cACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAg                                <  2:252721/145‑1 (MQ=255)
                                                                                                                     cACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATAGTTGGCGGGGAATACGCATTCGGACACAACCACGTGCGCTCGCAAACCACTCGTGGCAAGGCGt                            <  6:214401/149‑1 (MQ=255)
                                                                                                                      aCCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGtt                           >  6:258296/1‑149 (MQ=255)
                                                                                                                        cTC‑‑‑‑CGTGAGGGTCTACCCGGATGTAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTa                         <  6:90224/149‑1 (MQ=255)
                                                                                                                        cTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGttt                          >  1:365067/1‑148 (MQ=255)
                                                                                                                        cTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGttt                          >  1:85425/1‑148 (MQ=255)
                                                                                                                          c‑‑‑‑cGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAAc                       >  6:159292/1‑149 (MQ=255)
                                                                                                                               cGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAAcc                      <  6:57516/149‑1 (MQ=255)
                                                                                                                               cGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAAc                       <  3:230238/148‑1 (MQ=255)
                                                                                                                                cgggAGGGTCTACGCGGATGCAGACACCCTGCCCGAGTTGGGCCGAACGGTGAGCCTTGAAGAATAGGGCACCGTCATCGTTGGCGGGGAAGAGGCATTCGGACTCAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTaa                        <  3:316121/144‑1 (MQ=255)
                                                                                                                                  tGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCg                     >  4:39185/1‑148 (MQ=255)
                                                                                                                                     gggTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAAtt                 <  2:389516/149‑1 (MQ=255)
                                                                                                                                       gTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTcc               <  3:223751/149‑1 (MQ=255)
                                                                                                                                       gTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTcc               <  3:228613/149‑1 (MQ=255)
                                                                                                                                            cccGGATGCAGACACCCTGACCGAGTTCGGCCGTACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTaa           <  2:328879/148‑1 (MQ=255)
                                                                                                                                            cccGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAAc          >  5:152423/1‑149 (MQ=255)
                                                                                                                                                    cAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGc  >  3:22752/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTCACCATGAAGAAGGGTGAATTCATCATTTTCTGGTCCACCCTCCTCCATGGTTCGCTTCCTCATCTGGGAGAGACGACCACGCCGCGCCTGGGCTACGTTTCGCGCTACGTGCCCACCTC‑‑‑‑CGTGAGGGTCTACCCGGATGCAGACACCCTGACCGAGTTCGGCCGAACGGTGAGCCTTGAACAATACGGCACCGTCATCGTTGGCGGGGAAGACGCATTCGGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGC  >  NC_002947/4310942‑4311233

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: