Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,788 A→T 6.1% L387Q (CTG→CAG)  PP_0222 ← DszA family monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,7880AT 6.1% 442.6 / 12.7 132L387Q (CTG→CAG) PP_0222DszA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (76/48);  new base T (4/4);  total (80/52)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 133 total aligned reads displayed.

CGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGG  >  NC_002947/274645‑274937
                                                                                                                                               |                                                                                                                                                     
cGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCgg                                                                                                                                                 >  1:1701575/1‑150 (MQ=255)
   ccTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGttt                                                                                                                                              <  2:4511901/150‑1 (MQ=255)
    cTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTc                                                                                                                                             >  1:2757823/1‑150 (MQ=255)
    cTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTc                                                                                                                                             >  1:2757534/1‑150 (MQ=255)
     tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                            >  1:4246272/1‑150 (MQ=255)
     tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                            <  1:789568/150‑1 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCAt                                                                                                                                           >  2:1696189/1‑150 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCAt                                                                                                                                           >  2:1696166/1‑150 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaa                                                                                                                                           <  2:4699607/150‑2 (MQ=255)
      ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaa                                                                                                                                           <  2:2055969/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)
              ttttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAc                                                                                                                                   <  2:1785953/150‑1 (MQ=255)
                ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                                 <  1:370691/150‑1 (MQ=255)
                ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                                 <  2:252710/150‑1 (MQ=255)
                ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                                 <  1:3322928/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)
                             ccTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCg                                                                                                                    <  2:2076985/150‑1 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                  >  1:4000537/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                  >  1:4632538/1‑150 (MQ=255)
                               ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                                  >  1:4266262/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)
                                         gCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCa                                                                                                        >  2:4320253/1‑150 (MQ=255)
                                            gTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                                     <  1:4260332/150‑1 (MQ=255)
                                            gTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                                     >  2:796342/1‑150 (MQ=255)
                                             tCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcg                                                                                                    >  2:1628811/1‑150 (MQ=255)
                                               ttGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcg                                                                                                  >  2:4431024/1‑150 (MQ=255)
                                                tGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAACTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcgc                                                                                                 <  2:1699055/150‑1 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                >  1:2072727/1‑150 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                <  2:1612845/150‑1 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                >  1:2072643/1‑150 (MQ=255)
                                                 gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                                >  1:4916523/1‑150 (MQ=255)
                                                      ggCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCa                                                                                           >  1:3416150/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)
                                                                           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                                                               <  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)
                                                                                              gTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAg                                                   >  1:1571870/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)
                                                                                                    gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACa                                                    >  1:3679648/1‑143 (MQ=255)
                                                                                                    gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                             >  1:3141096/1‑150 (MQ=255)
                                                                                                    gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                             >  2:3461181/1‑150 (MQ=255)
                                                                                                    gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                             >  2:2946655/1‑150 (MQ=255)
                                                                                                     aaGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcgg                                            >  1:2342517/1‑150 (MQ=255)
                                                                                                        ctctCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcggcgg                                         <  1:1740908/150‑1 (MQ=255)
                                                                                                         tctcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCggg                                        >  1:4808296/1‑150 (MQ=255)
                                                                                                         tctcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCggg                                        >  2:3912873/1‑150 (MQ=255)
                                                                                                           tcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTc                                      >  1:3953016/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:2757534/150‑1 (MQ=255)
                                                                                                                        tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                         <  2:2757823/150‑1 (MQ=255)
                                                                                                                        tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                         <  2:90198/150‑1 (MQ=255)
                                                                                                                         gCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTg                        <  2:4221398/150‑1 (MQ=255)
                                                                                                                         gCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTg                        <  2:338383/150‑1 (MQ=255)
                                                                                                                           gggTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCa                      >  2:3180416/1‑150 (MQ=255)
                                                                                                                              tCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACt                   >  1:3029441/1‑150 (MQ=255)
                                                                                                                              gatcgatGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACt                   >  2:2779477/7‑150 (MQ=255)
                                                                                                                                aGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTgg                 <  2:4093824/150‑1 (MQ=255)
                                                                                                                                aGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTgg                 >  1:111365/1‑150 (MQ=255)
                                                                                                                                   atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc              >  1:3244420/2‑150 (MQ=255)
                                                                                                                                   atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc              >  2:1084784/2‑150 (MQ=255)
                                                                                                                                      aaaCCGGCCTGGACGGTTTCAACGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTg           >  2:3246097/1‑150 (MQ=255)
                                                                                                                                       aaCCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTgg          >  1:3372113/1‑150 (MQ=255)
                                                                                                                                         ccGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCt        <  2:3572594/150‑1 (MQ=255)
                                                                                                                                         ccGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCt        <  1:2438431/150‑1 (MQ=255)
                                                                                                                                          cGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTg       >  1:416516/1‑150 (MQ=255)
                                                                                                                                          cGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTg       >  1:416641/1‑150 (MQ=255)
                                                                                                                                           gTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGa      <  2:4916523/150‑1 (MQ=255)
                                                                                                                                               aGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTgg  >  2:2372434/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGACTGG  >  NC_002947/274645‑274937

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: