Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 799,838 A→T 5.4% L213Q (CTG→CAG)  ispB ← octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947799,8380AT 5.4% 475.4 / 10.0 145L213Q (CTG→CAG) ispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (74/62);  new base T (5/3);  total (80/65)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

CACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATG  >  NC_002947/799689‑799987
                                                                                                                                                     |                                                                                                                                                     
cACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCcagca                                                                                                                                                       <  2:1913122/150‑1 (MQ=255)
cACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCcagca                                                                                                                                                       >  1:1985714/1‑150 (MQ=255)
   ccTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgt                                                                                                                                                    >  1:416546/1‑150 (MQ=255)
    cttcttGATCGCCTGGCTTACCAGCGCCGCCTGCTCGGCAGTGACTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCATGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtc                                                                                                                                                   <  2:1014516/150‑1 (MQ=255)
     ttCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcg                                                                                                                                                  >  1:76831/1‑150 (MQ=255)
      tCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgt                                                                                                                                                 >  1:602818/1‑150 (MQ=255)
       cTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCTAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgtc                                                                                                                                                <  2:449271/150‑1 (MQ=255)
       cTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGgtcgac                                                                                                                                                <  1:4154315/150‑3 (MQ=255)
       cTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgtc                                                                                                                                                <  1:1098108/150‑1 (MQ=255)
         ggATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGa                                                                                                                                              >  1:3925006/1‑150 (MQ=255)
         ggATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGa                                                                                                                                              >  2:2496367/1‑150 (MQ=255)
          gATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGAc                                                                                                                                             >  1:1419488/1‑150 (MQ=255)
           gtggCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGAcc                                                                                                                                            <  1:891648/147‑1 (MQ=255)
            tCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCa                                                                                                                                           <  1:5000970/150‑1 (MQ=255)
             cGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAg                                                                                                                                          <  2:3760378/150‑1 (MQ=255)
              gCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGc                                                                                                                                         <  2:4109065/150‑1 (MQ=255)
              gCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGc                                                                                                                                         <  2:2887704/150‑1 (MQ=255)
                cTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTg                                                                                                                                       <  1:190550/150‑1 (MQ=255)
                   gCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGaa                                                                                                                                    >  1:5082673/1‑150 (MQ=255)
                     gTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCGGGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                  >  2:5125808/1‑150 (MQ=255)
                           gcgcCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGc                                                                                                                            <  2:1702249/150‑1 (MQ=255)
                             gccgccTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGccc                                                                                                                          >  1:2115711/1‑150 (MQ=255)
                                   tGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATgg                                                                                                                    >  1:4507377/1‑150 (MQ=255)
                                      tCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCg                                                                                                                 >  1:1416385/1‑150 (MQ=255)
                                        ggCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGcc                                                                                                               >  1:2514123/1‑150 (MQ=255)
                                           aGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCaaa                                                                                                            <  2:4739604/150‑1 (MQ=255)
                                            gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                           <  1:3015409/150‑1 (MQ=255)
                                            gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                           <  2:3304752/150‑1 (MQ=255)
                                            gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                           <  1:2090110/150‑1 (MQ=255)
                                            gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                           <  2:1534484/150‑1 (MQ=255)
                                             tGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCGTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAgg                                                                                                          >  2:2904715/1‑150 (MQ=255)
                                              gCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGGCCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGt                                                                                                         >  2:650905/1‑150 (MQ=255)
                                               ccTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTa                                                                                                        >  1:2908469/1‑150 (MQ=255)
                                                         tGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCtt                                                                                              <  1:4821831/150‑1 (MQ=255)
                                                         tGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCtt                                                                                              <  1:4821834/150‑1 (MQ=255)
                                                              tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCgcgc                                                                                         <  1:2217562/150‑1 (MQ=255)
                                                              tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCgcgc                                                                                         >  2:1294649/1‑150 (MQ=255)
                                                              tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCgcgc                                                                                         >  2:3837331/1‑150 (MQ=255)
                                                              tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCgcgc                                                                                         <  1:2385068/150‑1 (MQ=255)
                                                                  tCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAg                                                                                                   <  2:2250213/136‑1 (MQ=255)
                                                                  tCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAg                                                                                                   >  1:2250213/1‑136 (MQ=255)
                                                                   cAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCtc                                                                                    >  1:1256580/1‑150 (MQ=255)
                                                                    aGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCCAAGGTACGCAGGGCTTCGCGCTGCtcg                                                                                   >  1:3278031/1‑150 (MQ=255)
                                                                    aGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCtcg                                                                                   >  2:490472/1‑150 (MQ=255)
                                                                      cGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCtcgtc                                                                                 <  1:4334802/150‑1 (MQ=255)
                                                                       ggCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCtcgtcg                                                                                <  1:3549417/150‑1 (MQ=255)
                                                                          aGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTg                                                                             <  1:2427244/150‑1 (MQ=255)
                                                                            ggTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGc                                                                           <  1:88927/150‑1 (MQ=255)
                                                                            ggTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGc                                                                           <  1:3032484/150‑1 (MQ=255)
                                                                             gTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGcc                                                                          <  1:2147765/150‑1 (MQ=255)
                                                                                ggTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTgg                                                                       <  1:557418/150‑1 (MQ=255)
                                                                                  tttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGc                                                                     >  2:3421607/1‑150 (MQ=255)
                                                                                  tttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCGTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGc                                                                     >  2:582804/1‑150 (MQ=255)
                                                                                   ttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGcc                                                                    >  2:207263/1‑150 (MQ=255)
                                                                                    tGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCa                                                                   <  2:3373509/150‑1 (MQ=255)
                                                                                    tGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCa                                                                   <  1:1854105/150‑1 (MQ=255)
                                                                                      cccTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAgc                                                                 >  1:3667988/1‑150 (MQ=255)
                                                                                        cTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAgcgc                                                               >  2:715768/1‑150 (MQ=255)
                                                                                        cTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAgcgc                                                               >  1:1679982/1‑150 (MQ=255)
                                                                                         tCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGcgcc                                                              >  2:4020914/1‑150 (MQ=255)
                                                                                         tCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGCTGGCCCAGGCCAccgac                                                              <  2:5131253/150‑3 (MQ=255)
                                                                                          cGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGcgccg                                                             <  2:2575630/150‑1 (MQ=255)
                                                                                           ggCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGcgacga                                                            <  1:3499981/150‑2 (MQ=255)
                                                                                            gCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgca                                                           <  2:5424927/150‑1 (MQ=255)
                                                                                            gCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgca                                                           <  2:5424924/150‑1 (MQ=255)
                                                                                               aGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgcagca                                                        <  2:76831/150‑1 (MQ=255)
                                                                                                 gtcgtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCAcg                                                      >  1:386786/1‑149 (MQ=255)
                                                                                                 gtcgtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACt                                                      >  2:2688496/1‑150 (MQ=255)
                                                                                                  tcgtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCGTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtg                                                     >  2:3102638/1‑150 (MQ=255)
                                                                                                   cgtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                    <  2:3676217/150‑1 (MQ=255)
                                                                                                    gtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                   <  1:582275/150‑1 (MQ=255)
                                                                                                     tcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTgg                                                  <  1:2223196/150‑1 (MQ=255)
                                                                                                       aCCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGt                                                >  2:4670152/1‑150 (MQ=255)
                                                                                                       aCCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGt                                                <  2:4378250/150‑1 (MQ=255)
                                                                                                          gACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGa                                             >  1:3210686/1‑150 (MQ=255)
                                                                                                           aCGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCACCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGaa                                            <  2:4025340/150‑1 (MQ=255)
                                                                                                                cgtGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                       <  1:2303732/148‑1 (MQ=255)
                                                                                                                  tGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGa                                     >  1:2765281/1‑150 (MQ=255)
                                                                                                                   gCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGaa                                    <  2:1397181/150‑1 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                                  >  2:4941852/1‑150 (MQ=255)
                                                                                                                      cAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAg                                 >  2:2743178/1‑150 (MQ=255)
                                                                                                                       aGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGc                                >  2:1029654/1‑150 (MQ=255)
                                                                                                                          gTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATc                             >  1:1697248/1‑150 (MQ=255)
                                                                                                                           tCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCg                            >  2:413291/1‑150 (MQ=255)
                                                                                                                                    cACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGc                   <  2:4971594/150‑1 (MQ=255)
                                                                                                                                    cACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGc                   <  2:2238228/150‑1 (MQ=255)
                                                                                                                                     aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGcc                  >  1:525056/1‑150 (MQ=255)
                                                                                                                                     aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGcc                  <  1:3734115/150‑1 (MQ=255)
                                                                                                                                     aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGcc                  >  1:4887449/1‑150 (MQ=255)
                                                                                                                                        cTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCcgcg               >  2:722205/1‑150 (MQ=255)
                                                                                                                                         ttGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGa              >  2:3639708/1‑150 (MQ=255)
                                                                                                                                          ggccttCCAGCTGGTCGACGACCTGCTGGACGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                       >  2:1563152/6‑124 (MQ=39)
                                                                                                                                           gccttCCAGCTGGTCGACGACCTGCTGGACGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAAt            >  2:3672259/5‑150 (MQ=255)
                                                                                                                                             aGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGa          >  1:317017/1‑150 (MQ=255)
                                                                                                                                               tCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGAcc        >  1:4917932/1‑150 (MQ=255)
                                                                                                                                                 cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTc      >  1:1961919/1‑150 (MQ=255)
                                                                                                                                                  agcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTcc     <  1:426402/150‑1 (MQ=255)
                                                                                                                                                   gcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgca                                                           >  1:257722/1‑95 (MQ=255)
                                                                                                                                                   gcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgca                                                           <  2:257722/95‑1 (MQ=255)
                                                                                                                                                   gcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCa    >  2:1800864/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATG  >  NC_002947/799689‑799987

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