Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,720,658 G→C 5.4% intergenic (+34/+97) fabB → / ← PP_4176 3‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,720,6580GC 5.4% 628.8 / 18.1 183intergenic (+34/+97)fabB/PP_41763‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase
Reads supporting (aligned to +/- strand):  ref base G (81/92);  new base C (5/5);  total (86/97)
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

Only 100 of 190 total aligned reads displayed.

TGCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTTGT  >  NC_002947/4720509‑4720807
                                                                                                                                                     |                                                                                                                                                     
 gCCAATCCTGCTTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTgg                                                                                                                                                      <  1:797906/150‑1 (MQ=255)
 gCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTgg                                                                                                                                                      >  1:3532412/1‑150 (MQ=255)
 gCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTgg                                                                                                                                                      >  2:2354011/1‑150 (MQ=255)
       ccTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGAAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCggg                                                                                                                                                <  2:1944837/150‑1 (MQ=255)
         tGCTTAAAACCTAAGAAATCCCCAATATCGACACGGTCATGAGCAAAATGTTCTGCTTCGGCGGCACCAACGCCACCCTTGTGCTCAAGCGCTTGTAAGGCATGTAAGACTCTGATGATGTAAATGAATACGCCCCGACTGGTTCGGGGc                                                                                                                                              <  1:3192829/150‑1 (MQ=255)
         tGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGc                                                                                                                                              >  1:2079793/1‑150 (MQ=255)
         tGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGc                                                                                                                                              >  1:2956604/1‑150 (MQ=255)
           cttAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGt                                                                                                                                            <  1:3061193/148‑1 (MQ=255)
            gTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGAACCAGTCGGGGCGtt                                                                                                                                           <  1:1078304/150‑1 (MQ=255)
             tAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGttt                                                                                                                                          >  2:2897294/1‑150 (MQ=255)
                aaCCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttg                                                                                                                                       >  1:4404582/1‑150 (MQ=255)
                 aCCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgt                                                                                                                                      >  2:619757/1‑150 (MQ=255)
                  ccGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgtt                                                                                                                                     >  1:3137864/1‑150 (MQ=255)
                   cGAAGAACCCGCCACGATCGCCCCGGTCATGCGCAACAGCTTCGGCTTCGGCGGCCCCAACGCCACCCTGGTGCTCACGCGCTGGGAAGGCAAGTAAGACGCTGCGGCTGCAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgttt                                                                                                                                    <  2:1156999/150‑1 (MQ=255)
                   cGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgttt                                                                                                                                    <  1:1585616/150‑1 (MQ=255)
                   cGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgttt                                                                                                                                    >  2:338159/1‑150 (MQ=255)
                    gaagaaAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgtttg                                                                                                                                   >  2:3943292/1‑150 (MQ=255)
                     aagaaAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTgg                                                                                                                                  <  1:1552328/150‑1 (MQ=255)
                     aacaaaaCGCCAAGACCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTgg                                                                                                                                  <  1:1513867/147‑1 (MQ=255)
                        aaaaCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTTTTTGGGtt                                                                                                                               >  1:476626/1‑150 (MQ=255)
                         aaaCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTg                                                                                                                              <  1:4490605/150‑1 (MQ=255)
                            cGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCt                                                                                                                           <  2:2271201/150‑1 (MQ=255)
                             gCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCtg                                                                                                                          >  2:1118101/1‑150 (MQ=255)
                                aaGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGc                                                                                                                       >  1:862047/1‑150 (MQ=255)
                                    tCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGc                                                                                                                   <  2:816873/150‑1 (MQ=255)
                                    tCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGc                                                                                                                   <  2:2762668/150‑1 (MQ=255)
                                    tCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACACCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGc                                                                                                                   <  2:1527294/150‑1 (MQ=255)
                                       acacGGTCATGAGCAACGGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCctc                                                                                                                >  2:2344272/1‑150 (MQ=255)
                                         acGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCtt                                                                                                              >  2:693585/1‑150 (MQ=255)
                                            gTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           <  1:3384322/150‑1 (MQ=255)
                                            gTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           >  1:1315444/1‑150 (MQ=255)
                                            gTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           >  2:3036384/1‑150 (MQ=255)
                                              cATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCgg                                                                                                         <  1:2583726/150‑1 (MQ=255)
                                                 gAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTa                                                                                                      >  2:4053534/1‑150 (MQ=255)
                                                  aGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTaa                                                                                                     >  1:46665/1‑150 (MQ=255)
                                                    cAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   >  1:509332/1‑150 (MQ=255)
                                                     aaCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATTAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGTCCTCTTCGCGGGTAAAcc                                                                                                  <  2:2689957/150‑1 (MQ=255)
                                                      aCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGTTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                                                                 >  2:1267816/1‑150 (MQ=255)
                                                      aCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGAACCAGTCGGGGCGtt                                                                                                                                           <  2:2331382/108‑1 (MQ=39)
                                                           ttcggcttcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                            >  2:3421924/1‑150 (MQ=255)
                                                             cggcttcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                                          <  2:3101442/150‑1 (MQ=255)
                                                               gcttcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                        >  1:4077569/1‑150 (MQ=255)
                                                                  tcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTTTTTTGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                     >  1:752942/1‑150 (MQ=255)
                                                                  tcggcGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCGTCGCGGGTAAACCCGCTCCTACACAg                                                                                     >  1:2279248/1‑150 (MQ=255)
                                                                    ggcggcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGt                                                                                   <  1:4518516/150‑1 (MQ=255)
                                                                     gcggcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTc                                                                                  >  1:4031755/1‑150 (MQ=255)
                                                                          aCCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTg                                                                             <  1:3616633/150‑1 (MQ=255)
                                                                          aCCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTg                                                                             >  1:3651988/1‑150 (MQ=255)
                                                                          aCCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTg                                                                             >  1:1929204/1‑150 (MQ=255)
                                                                            cAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCa                                                                           <  1:1509715/150‑1 (MQ=255)
                                                                               cGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTg                                                                        <  2:4321305/150‑1 (MQ=255)
                                                                                  cACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCtt                                                                     >  1:428938/1‑150 (MQ=255)
                                                                                  cACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCtt                                                                     <  2:3053013/150‑1 (MQ=255)
                                                                                   aCCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTtg                                                                    <  2:1241775/150‑1 (MQ=255)
                                                                                    cccTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTtgt                                                                   <  1:2412332/150‑1 (MQ=255)
                                                                                    cccTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGATTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTtgt                                                                   <  2:3585825/150‑1 (MQ=255)
                                                                                     ccTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCTACACAGGTCATTTGCAGTGCTtgtg                                                                  >  1:184426/1‑150 (MQ=255)
                                                                                       tGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCGTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGc                                                                >  1:1781464/1‑150 (MQ=255)
                                                                                         gTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCa                                                              >  1:3422322/1‑150 (MQ=255)
                                                                                          tGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTTGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGGCAt                                                             >  1:4160748/1‑150 (MQ=255)
                                                                                            cTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTg                                                           <  2:2810527/150‑1 (MQ=255)
                                                                                             tCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGACAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGc                                                          <  2:165889/150‑1 (MQ=255)
                                                                                             tCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGc                                                          >  1:1476771/1‑150 (MQ=255)
                                                                                             tCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGc                                                          <  2:3250446/150‑1 (MQ=255)
                                                                                             tCAAGAGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGACTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGATTGTGGCCATTGc                                                          >  2:215309/1‑150 (MQ=255)
                                                                                              cAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCa                                                         <  1:3307268/150‑1 (MQ=255)
                                                                                               aaGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAg                                                        >  1:2199514/1‑150 (MQ=255)
                                                                                                aGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAgg                                                       >  2:1303866/1‑150 (MQ=255)
                                                                                                 gcgcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAATCGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGc                                                      >  2:2190364/1‑150 (MQ=255)
                                                                                                       ggAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATgg                                                <  1:3712810/150‑1 (MQ=255)
                                                                                                       ggAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATgg                                                <  2:1372995/150‑1 (MQ=255)
                                                                                                        gAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGgc                                               <  2:1774270/150‑1 (MQ=255)
                                                                                                          aGGCAAGTAAGACGCTGCTGCTGTAAATGACAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGgcgc                                             <  2:2607171/150‑1 (MQ=255)
                                                                                                              aaGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCt                                         >  2:364190/1‑150 (MQ=255)
                                                                                                                 tAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTgg                                      <  2:4182376/150‑1 (MQ=255)
                                                                                                                 tAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTgg                                      >  2:3529134/1‑150 (MQ=255)
                                                                                                                 tAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTgg                                      <  2:3303716/150‑1 (MQ=255)
                                                                                                                  aaGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTggg                                     <  1:4332485/150‑1 (MQ=255)
                                                                                                                  aaGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTGTTTGTTTGGGTTGCCGGTGCTGGCCTCGTCGCGGGTAAACCCGCTCCGACACAGGTCATTTGCAGTGCGTGTGGCCCTTGCAGGCTACTGGCGCGCCGTggg                                     >  1:810482/1‑150 (MQ=38)
                                                                                                                     aCGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCaa                                  <  1:2887479/150‑1 (MQ=255)
                                                                                                                      cGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAac                                 <  1:895443/150‑1 (MQ=255)
                                                                                                                            tgctgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCAATGCAGGCTAATGGCGCGCCATGGGCAACACCGcc                           >  1:1995116/1‑150 (MQ=255)
                                                                                                                            tgctgTAAATGAAAACGCCCCGACTGGGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGcc                           >  2:3927621/1‑150 (MQ=255)
                                                                                                                             gctgTAACTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCa                          <  1:2354011/150‑1 (MQ=255)
                                                                                                                               tgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCtgt                                                                                                                         >  2:2715996/1‑53 (MQ=255)
                                                                                                                               tgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCtgt                                                                                                                         <  1:2715996/53‑1 (MQ=255)
                                                                                                                               tgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGGAAACCCGCTCCTACACAGGTCATTTTCAGTGCGTGTGGCCCTTGCAGGCTAAATGCGCGCCTTGGGCAACACCCCCCGa                        >  1:757853/1‑150 (MQ=38)
                                                                                                                                  aaaTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGaaaa                     <  2:1511511/150‑1 (MQ=255)
                                                                                                                                   aaTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTTCAGGCTAATTGCGCGCCTTGGGCAACACCGCCAGAAAAt                    >  1:3282977/1‑150 (MQ=255)
                                                                                                                                   aaTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAAt                    >  2:474232/1‑150 (MQ=255)
                                                                                                                                   aaTGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGGACACCCGCTCCTACACACGTCATTTGCAGTGCTTGTGGCCATTTCAGGCTAATGGCGCGCCTTGGGCCACACCGCCAGAAAAt                    >  1:2699652/1‑150 (MQ=39)
                                                                                                                                       aaaaCGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATtctc                <  2:4404582/150‑1 (MQ=255)
                                                                                                                                        aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTcc               >  2:1133035/1‑150 (MQ=38)
                                                                                                                                        aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTcc               >  2:3169076/1‑150 (MQ=38)
                                                                                                                                        aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCGCTTCGCGGGGAAACCCGCTCCTACACCGGGCCTTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCGTGGGGAACACCGCCAGCAAAATCTcc               >  1:56043/1‑150 (MQ=255)
                                                                                                                                          aCGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCg             <  2:2266545/150‑1 (MQ=255)
                                                                                                                                            gCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGaaa                      <  1:2060096/139‑1 (MQ=255)
                                                                                                                                            gCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGaaa                      >  2:2060096/1‑139 (MQ=255)
                                                                                                                                              cccGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAAGTCTCCCGGGcc         >  2:1469647/1‑150 (MQ=255)
                                                                                                                                                     ggTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTtgt  <  1:693585/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTTGT  >  NC_002947/4720509‑4720807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: