Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,726,087 C→G 6.7% A247G (GCC→GGC)  PP_3291 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,726,0870CG 6.7% 190.4 / 3.5 60A247G (GCC→GGC) PP_3291metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

AGGACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAGT  >  NC_002947/3725939‑3726233
                                                                                                                                                     |                                                                                                                                                  
aGGACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCgcgc                                                                                                                                                    >  7:43381/1‑149 (MQ=255)
  gACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGccc                                                                                                                                                  <  6:250196/149‑1 (MQ=255)
    cggcggACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCaa                                                                                                                                                <  5:119583/149‑1 (MQ=255)
      gcggACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACt                                                                                                                                              <  8:230883/149‑1 (MQ=255)
        ggACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGGCGCGGCCAACGGc                                                                                                                                            <  6:263647/149‑1 (MQ=255)
        ggACGTGCTGTTGATCGAATCTACCTACGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGGCGCGGCCAACGGc                                                                                                                                            <  8:339003/149‑1 (MQ=255)
           cGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCa                                                                                                                                         <  5:380299/149‑1 (MQ=255)
               cTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGt                                                                                                                                      <  6:2100/148‑1 (MQ=255)
                    taTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCTACCG‑ATCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTAc                                                                                                                                 <  7:169400/147‑1 (MQ=255)
                    gATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTAc                                                                                                                                 <  7:7499/148‑1 (MQ=255)
                     aTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTAcc                                                                                                                                <  8:101343/148‑1 (MQ=255)
                     aTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTAcc                                                                                                                                <  5:324543/148‑1 (MQ=255)
                     aTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCt                                                                                                                               >  8:217976/1‑149 (MQ=255)
                     aTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCt                                                                                                                               <  8:43381/149‑1 (MQ=255)
                     aTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACATTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCt                                                                                                                               <  5:330114/149‑1 (MQ=255)
                      tCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTg                                                                                                                              >  3:151353/1‑149 (MQ=255)
                      tCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTg                                                                                                                              >  3:203407/1‑149 (MQ=255)
                               ccTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGc                                                                                                                       <  1:195387/147‑1 (MQ=255)
                                cTATGGCGATCGCCAGCATCCCGCCGAATCGACCGAAGCGGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGGGTCTTGGCATTACCCGCGTCCCCTCGTTCGCCGTTGGCCGCGCCGACCTGGGCATGTATTACATGTACCCGCt                                                                                                                      >  1:252622/1‑148 (MQ=255)
                                 tATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCaa                                                                                                                   >  4:148491/1‑149 (MQ=255)
                                     gCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAgcgc                                                                                                               >  7:340661/1‑149 (MQ=255)
                                               gCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTg                                                                                                      >  8:140126/1‑148 (MQ=255)
                                                   cccGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCgg                                                                                                         >  4:174966/1‑141 (MQ=255)
                                                   cccGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCgg                                                                                                         <  3:174966/141‑1 (MQ=255)
                                                      gccgAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGa                                                                                              <  4:203407/149‑1 (MQ=255)
                                                       ccgAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGat                                                                                             <  7:140126/149‑1 (MQ=255)
                                                       ccgAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGat                                                                                             <  5:95047/149‑1 (MQ=255)
                                                         gAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCACGCGCGACCGGCTGATCCCTGatat                                                                                           >  4:117781/1‑149 (MQ=255)
                                                            tCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATccc                                                                                        >  1:43450/1‑149 (MQ=255)
                                                              gACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGa                                                                                      <  2:194381/149‑1 (MQ=255)
                                                               aCCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGa                                                                                      >  5:32508/1‑148 (MQ=255)
                                                                ccG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATc                                                                                    <  7:217976/149‑1 (MQ=255)
                                                                ccG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATc                                                                                    >  1:352529/1‑149 (MQ=255)
                                                                 cG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATc                                                                                    >  8:160092/1‑148 (MQ=255)
                                                                 cG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATc                                                                                    >  7:262111/1‑148 (MQ=255)
                                                                 cG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCt                                                                                   <  5:265857/149‑1 (MQ=255)
                                                                      cAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCt                                                                               <  3:270295/149‑1 (MQ=255)
                                                                        gTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTg                                                                              >  7:333488/1‑148 (MQ=255)
                                                                         tACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGaa                                                                            <  6:32508/149‑1 (MQ=255)
                                                                         tACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGa                                                                             >  8:327091/1‑148 (MQ=255)
                                                                           ccTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGaa                                                                            >  5:256624/1‑147 (MQ=255)
                                                                           ccTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACa                                                                          >  5:258840/1‑149 (MQ=255)
                                                                           ccTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCACTGATATCCCGATCTACCTGAACa                                                                          <  8:262111/149‑1 (MQ=255)
                                                                               gCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTcc                                                                      >  3:243743/1‑149 (MQ=255)
                                                                               gCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTc                                                                       <  6:258840/148‑1 (MQ=255)
                                                                                 cAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGa                                                                    <  3:148491/149‑1 (MQ=255)
                                                                                  aaGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGAt                                                                   <  2:352529/149‑1 (MQ=255)
                                                                                        atcaACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCAc                                                             >  4:165452/1‑149 (MQ=255)
                                                                                           aaCGAAACGCTCCTCCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTCATATCCCGATCTACCTCAACAGTCCGATGGCCACCGa                                                          >  2:22929/1‑149 (MQ=255)
                                                                                                  cGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCAcc                                                   <  4:151353/149‑1 (MQ=255)
                                                                                                        tGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGATGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATGCCTCATATCCCGATCTCCCTGAACTGTCCGAAGGCCACCGACGCAAGCGCCCTg                                             >  4:234308/1‑149 (MQ=255)
                                                                                                              aTGGCATTACCCTCGTAACCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACcag                                       >  2:97131/1‑149 (MQ=255)
                                                                                                                 gCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACcagcag                                    <  3:165452/149‑1 (MQ=255)
                                                                                                                  cATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACcagcag                                    >  2:359837/1‑148 (MQ=255)
                                                                                                                   aTTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTa                                  <  7:327091/149‑1 (MQ=255)
                                                                                                                                    gtTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCgg                  <  7:160092/148‑1 (MQ=255)
                                                                                                                                     ttCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCt                >  7:315849/1‑149 (MQ=255)
                                                                                                                                       agcagTTGGGCGCGGCCAACGGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGa              >  1:286539/5‑149 (MQ=255)
                                                                                                                                       agcagTTGGGCGCGGCCAACGGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGa              >  2:311488/5‑149 (MQ=255)
                                                                                                                                         ccGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGAcc            >  6:35720/1‑149 (MQ=255)
                                                                                                                                               gCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACACCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGc       <  8:333488/148‑1 (MQ=255)
                                                                                                                                                 cgcgCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGa    <  3:117781/149‑1 (MQ=255)
                                                                                                                                                  gcgcCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAg   >  6:62528/1‑149 (MQ=255)
                                                                                                                                                   cgcCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAGt  <  2:43450/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGGACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAGT  >  NC_002947/3725939‑3726233

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑