Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,791 C→A 5.8% G386V (GGC→GTC)  PP_0222 ← DszA family monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,7910CA 5.8% 459.5 / 12.9 138G386V (GGC→GTC) PP_0222DszA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (79/51);  new base A (4/4);  total (83/55)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 138 total aligned reads displayed.

CGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGGAT  >  NC_002947/274645‑274939
                                                                                                                                                  |                                                                                                                                                    
cGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCgg                                                                                                                                                   >  1:1701575/1‑150 (MQ=255)
    cTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTc                                                                                                                                               >  1:2757534/1‑150 (MQ=255)
     tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                              >  1:4246272/1‑150 (MQ=255)
     tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                              <  1:3921561/150‑1 (MQ=255)
     tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                              <  2:4814971/150‑1 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCAt                                                                                                                                             >  2:1696166/1‑150 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCAt                                                                                                                                             >  2:1696189/1‑150 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaa                                                                                                                                             <  2:2055969/150‑2 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaa                                                                                                                                             <  2:4699607/150‑2 (MQ=255)
       gAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATc                                                                                                                                            >  1:326136/1‑150 (MQ=255)
       gAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaac                                                                                                                                            <  2:3817725/150‑3 (MQ=255)
        aaCAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGATGAAACCGGCCTGGACGGTTTCaacc                                                                                                                                           <  1:3438924/150‑4 (MQ=255)
         aCAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGa                                                                                                                                          >  2:1740908/1‑150 (MQ=255)
              ttttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAc                                                                                                                                     <  2:1785953/150‑1 (MQ=255)
               tttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGTGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcg                                                                                                                                    >  2:672128/1‑150 (MQ=255)
                ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                                   <  1:370691/150‑1 (MQ=255)
                  cgcgCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGa                                                                                                                                 <  1:1634712/150‑1 (MQ=255)
                    cgcAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATc                                                                                                                               >  1:865467/1‑150 (MQ=255)
                      cAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCaa                                                                                                                             <  1:3910173/150‑1 (MQ=255)
                      cAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCaa                                                                                                                             <  2:2259523/150‑1 (MQ=255)
                             acTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCACCTGCTCg                                                                                                                      >  2:2210058/2‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  1:4000537/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  1:4266262/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  1:4632538/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  2:1670011/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  2:4477636/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  1:2363513/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                    >  1:2341595/1‑150 (MQ=255)
                                  ggTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCAc                                                                                                                 <  2:1295264/150‑1 (MQ=255)
                                  accgaaGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCAc                                                                                                                 >  1:1775673/6‑150 (MQ=255)
                                      tAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGt                                                                                                             >  1:254742/1‑150 (MQ=255)
                                        ggCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTc                                                                                                           >  2:2425392/1‑150 (MQ=255)
                                            gTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                                       <  1:4260332/150‑1 (MQ=255)
                                            gTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                                       <  2:3184518/150‑1 (MQ=255)
                                             tCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcg                                                                                                      >  1:3376194/1‑150 (MQ=255)
                                               ttGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcg                                                                                                    >  1:5003642/1‑150 (MQ=255)
                                               ttGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcg                                                                                                    >  2:4431024/1‑150 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                  >  1:4916523/1‑150 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                  <  2:1612845/150‑1 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                  >  1:2072727/1‑150 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                  >  1:2072643/1‑150 (MQ=255)
                                                               gctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGt                                                                                    <  2:655674/150‑1 (MQ=255)
                                                                ctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTa                                                                                   <  1:1377299/150‑1 (MQ=255)
                                                                ctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTa                                                                                   >  1:302287/1‑150 (MQ=255)
                                                                ctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTa                                                                                   >  2:5102519/1‑150 (MQ=255)
                                                                 tgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAg                                                                                  >  1:3251173/1‑150 (MQ=255)
                                                                 tgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAg                                                                                  >  1:304552/1‑150 (MQ=255)
                                                                           ggATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCa                                                                        >  1:3137626/1‑150 (MQ=255)
                                                                              tGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAgg                                                                     >  1:2418662/1‑150 (MQ=255)
                                                                                 ccAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGt                                                                  >  1:41503/1‑150 (MQ=255)
                                                                                 ccAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGt                                                                  >  1:2932258/1‑150 (MQ=255)
                                                                                  cAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTg                                                                 <  1:2872644/150‑1 (MQ=255)
                                                                                  cAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTg                                                                 <  2:1339800/150‑1 (MQ=255)
                                                                                    gCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCt                                                               >  1:4032758/1‑150 (MQ=255)
                                                                                     cTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTg                                                              >  2:3722182/1‑150 (MQ=255)
                                                                                     cTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTg                                                              >  2:158688/1‑150 (MQ=255)
                                                                                      tCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGt                                                             >  2:3665201/1‑150 (MQ=255)
                                                                                       cGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGtt                                                            <  2:932012/150‑1 (MQ=255)
                                                                                       cGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGtt                                                            <  1:2379316/150‑1 (MQ=255)
                                                                                        gATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTg                                                           <  2:2933929/150‑1 (MQ=255)
                                                                                        gATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTg                                                           >  2:2628428/1‑150 (MQ=255)
                                                                                         aTGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGc                                                          >  1:1776229/1‑150 (MQ=255)
                                                                                             aGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACa                                                      >  1:2490220/1‑150 (MQ=255)
                                                                                                ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                                   <  2:1363798/150‑1 (MQ=255)
                                                                                                ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                                   <  2:1035117/150‑1 (MQ=255)
                                                                                                ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                                   <  1:578369/150‑1 (MQ=255)
                                                                                                ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                                   <  2:1364118/150‑1 (MQ=255)
                                                                                                 cTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgcg                                                  >  2:3557331/1‑150 (MQ=255)
                                                                                                  tCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgcgc                                                 >  2:4484217/1‑150 (MQ=255)
                                                                                                  tCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgcgc                                                 >  1:2106061/1‑150 (MQ=255)
                                                                                                    gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACa                                                      <  2:3679648/143‑1 (MQ=255)
                                                                                                    gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                               >  2:3461181/1‑150 (MQ=255)
                                                                                                     aaGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcgg                                              >  1:2342517/1‑150 (MQ=255)
                                                                                                         tctcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCggg                                          >  2:3912873/1‑150 (MQ=255)
                                                                                                           tcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTc                                        >  2:3557569/1‑150 (MQ=255)
                                                                                                                   gACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGtt                                <  2:2473533/150‑1 (MQ=255)
                                                                                                                   gACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGtt                                <  1:3572272/150‑1 (MQ=255)
                                                                                                                    aCGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTg                               >  1:1601845/1‑150 (MQ=255)
                                                                                                                      ggTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTc                             >  2:2987056/1‑150 (MQ=255)
                                                                                                                        tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                           >  2:2592067/1‑150 (MQ=255)
                                                                                                                        tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                           <  2:90198/150‑1 (MQ=255)
                                                                                                                        tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                           <  2:2757534/150‑1 (MQ=255)
                                                                                                                         gCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTg                          <  2:338383/150‑1 (MQ=255)
                                                                                                                         gCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTg                          <  2:4221398/150‑1 (MQ=255)
                                                                                                                           gggTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCa                        >  2:3180416/1‑150 (MQ=255)
                                                                                                                              tCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACt                     >  1:3577629/1‑150 (MQ=255)
                                                                                                                              gatcgatGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACt                     >  2:2779477/7‑150 (MQ=255)
                                                                                                                                aGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTgg                   >  1:111365/1‑150 (MQ=255)
                                                                                                                                aGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTgg                   <  2:4093824/150‑1 (MQ=255)
                                                                                                                                   ttGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc                >  2:420329/1‑150 (MQ=255)
                                                                                                                                   atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc                >  1:3244420/2‑150 (MQ=255)
                                                                                                                                   atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc                >  2:1084784/2‑150 (MQ=255)
                                                                                                                                    tGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGcgcg               >  2:2438431/1‑150 (MQ=255)
                                                                                                                                       aaCCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTgg            >  1:3372113/1‑150 (MQ=255)
                                                                                                                                         ccGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCt          <  2:3572594/150‑1 (MQ=255)
                                                                                                                                          cGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTg         >  1:416641/1‑150 (MQ=255)
                                                                                                                                          cGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTg         >  1:416516/1‑150 (MQ=255)
                                                                                                                                                ggCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGGa   >  1:4655655/1‑150 (MQ=255)
                                                                                                                                                ggCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGGa   <  1:796342/150‑1 (MQ=255)
                                                                                                                                                 gCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGGAt  <  1:1504590/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGGAT  >  NC_002947/274645‑274939

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: