Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,726,082 C→G 6.5% G245G (GGC→GGG)  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,0820CG 6.5% 196.5 / 3.3 61G245G (GGC→GGG) PP_3291metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base G (2/2);  total (29/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 = 8.30e-01

TGGTCAGGACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGC  >  NC_002947/3725934‑3726228
                                                                                                                                                     |                                                                                                                                                  
tGGTCAGGACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGc                                                                                                                                                    >  2:195387/1‑149 (MQ=255)
     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                                                                               >  7:262111/1‑148 (MQ=255)
                                                                      cG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATc                                                                               >  8:160092/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         >  2:311488/5‑149 (MQ=255)
                                                                                                                                            agcagTTGGGCGCGGCCAACGGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGa         >  1:286539/5‑149 (MQ=255)
                                                                                                                                              ccGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGAcc       >  6:35720/1‑149 (MQ=255)
                                                                                                                                                    gCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACACCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGc  <  8:333488/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGTCAGGACGGCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGC  >  NC_002947/3725934‑3726228

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: ‑