Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,726,093 T→G 6.5% L249R (CTG→CGG)  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,0930TG 6.5% 178.4 / 3.1 62L249R (CTG→CGG) PP_3291metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base G (2/2);  total (30/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 = 6.06e-01

GCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAGT  >  NC_002947/3725945‑3726233
                                                                                                                                                     |                                                                                                                                            
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:43381/149‑1 (MQ=255)
               aTCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCt                                                                                                                               >  8:217976/1‑149 (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                                                                                                         <  3:174966/141‑1 (MQ=255)
                                             cccGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCgg                                                                                                         >  4:174966/1‑141 (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                                                                                    >  1:352529/1‑149 (MQ=255)
                                                          ccG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATc                                                                                    <  7:217976/149‑1 (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              >  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)
                                                                                                                                           cgcgCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGa    <  3:117781/149‑1 (MQ=255)
                                                                                                                                            gcgcCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAg   >  6:62528/1‑149 (MQ=255)
                                                                                                                                             cgcCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAGt  <  2:43450/149‑1 (MQ=255)
                                                                                                                                                ccAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGa                         <  5:304143/123‑1 (MQ=255)
                                                                                                                                                ccAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGa                         >  6:304143/1‑123 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCGGACGTGCTGTTGATCGAATCTACCTATGGCGATCGCCAGCATCCCGCCGAATCGACCG‑AGCAGTACCTTGCCAAGGTCATCAACGAAACGCTCCTGCGTCATGGCATTACCCTCGTACCCTCGTTCGCCGTTGGCCGCGCCCAACTGCTCATGTATTACCTGTACACGCTCAAGCGCGACCGGCTGATCCCTGATATCCCGATCTACCTGAACAGTCCGATGGCCACCGACGCCACCGCCCTGTACCAGCAGTACCGCAGTGAGCATCGGCTGACCGCAGCCGAGT  >  NC_002947/3725945‑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: ‑