Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,982,617 C→G 8.2% P190R (CCG→CGG)  PP_1775 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,982,6170CG 8.2% 185.4 / 6.1 61P190R (CCG→CGG) PP_1775metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/3);  total (28/33)
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.10e-01

TCTTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGTTCCT  >  NC_002947/1982474‑1982762
                                                                                                                                               |                                                                                                                                                  
tCTTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGc                                                                                                                                               <  4:65625/149‑1 (MQ=255)
  ttGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAg                                                                                                                                             >  1:188660/1‑149 (MQ=255)
    gCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGcc                                                                                                                                           >  4:33677/1‑149 (MQ=255)
      cTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGcccc                                                                                                                                         <  4:251862/149‑1 (MQ=255)
       tGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGccccc                                                                                                                                        >  5:349501/1‑149 (MQ=255)
         tCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGAACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCACCGGGCGCAGcagcgg                                                                                                                                       <  5:329676/148‑6 (MQ=255)
            tATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGc                                                                                                                                   <  5:371225/149‑1 (MQ=255)
            tATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGGTCCGCTTTTGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGc                                                                                                                                   <  7:57708/149‑1 (MQ=255)
               tCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                >  5:147935/1‑149 (MQ=255)
               tCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                >  2:340825/1‑149 (MQ=255)
               tCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                >  3:332308/1‑149 (MQ=255)
                 gCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCACCGGGCGCAGcag                                                                                                                                          >  7:296343/1‑135 (MQ=255)
                 gCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCACCGGGCGCAGcag                                                                                                                                          <  8:296343/137‑3 (MQ=255)
                   cTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGccc                                                                                                                                                    >  4:206590/1‑125 (MQ=255)
                            gCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                >  7:38126/1‑136 (MQ=255)
                            gCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                <  8:38126/136‑1 (MQ=255)
                              aGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACgtgctggtgc                                                                                                                <  1:340825/149‑1 (MQ=255)
                                ggCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACgtgctggtgctg                                                                                                              <  2:188660/149‑1 (MQ=255)
                                   aTCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCACCGGGCGCAGcagcgg                                                                                                                                       >  1:76925/1‑117 (MQ=255)
                                   aTCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCACCGGGCGCAGcagcgg                                                                                                                                       <  2:76925/122‑6 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             <  4:279273/142‑1 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             <  4:369508/142‑1 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             >  3:369508/1‑142 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             >  3:279273/1‑142 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGaaa                                                                                                          <  6:349501/145‑1 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGaaa                                                                                                          <  6:147935/145‑1 (MQ=255)
                                        ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGaaa                                                                                                          <  7:170888/145‑1 (MQ=255)
                                                   ttATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCgg                                                                                                                                      >  2:101792/1‑107 (MQ=255)
                                                   ttATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCgg                                                                                                                                      <  1:101792/107‑1 (MQ=255)
                                                      tGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGcc                                                                                        <  3:33677/149‑1 (MQ=255)
                                                        tCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTg                                                                                      >  4:313195/1‑149 (MQ=255)
                                                           aGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCAc                                                                                   <  3:313195/149‑1 (MQ=255)
                                                             tGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCAccc                                                                                 >  7:122632/1‑149 (MQ=255)
                                                                gATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGCACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCgc                                                                               <  1:78651/149‑1 (MQ=255)
                                                                 atatGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCgcag                                                                             >  2:259252/1‑149 (MQ=255)
                                                                  tatGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCgcagg                                                                            >  8:267762/1‑149 (MQ=255)
                                                                    tGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCTATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACgtg                                                                                                                       <  2:241146/104‑1 (MQ=255)
                                                                    tGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACgtg                                                                                                                       >  1:241146/1‑104 (MQ=255)
                                                                      cagcaTGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGTAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAga                                                                        >  2:133483/1‑149 (MQ=255)
                                                                      cagcaTGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAgac                                                                       >  1:112333/1‑149 (MQ=255)
                                                                       agcaTGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGGGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAgac                                                                       <  5:36468/149‑1 (MQ=255)
                                                                             aGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCt                                                                                       >  1:208283/1‑127 (MQ=255)
                                                                             aGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCt                                                                                       <  2:208283/127‑1 (MQ=255)
                                                                             aGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGAccgc                                                                 >  2:332967/1‑149 (MQ=255)
                                                                                aaaCCAACAGCCGGTACGTGCTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGt                                                              >  8:360974/1‑149 (MQ=255)
                                                                                        aGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTg                                                     <  2:112333/149‑1 (MQ=255)
                                                                                         gCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACgtg                                                                                                                       <  4:84189/83‑1 (MQ=255)
                                                                                         gCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACgtg                                                                                                                       >  3:84189/1‑83 (MQ=255)
                                                                                           cGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGa                                                   <  4:309588/149‑1 (MQ=255)
                                                                                             gTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAg                                                 <  1:259252/149‑1 (MQ=255)
                                                                                             gTACGTGTTTTCGGGCGATCTGGGCGCTTGTGGCAACCGGCTGGTGGGCCCGTTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTAGGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAg                                                 <  8:122632/149‑1 (MQ=255)
                                                                                                          ggCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCg                                    >  1:334091/1‑149 (MQ=255)
                                                                                                           gCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGAcc                                                                                           >  2:153759/1‑93 (MQ=255)
                                                                                                           gCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGAc                                      <  7:321306/146‑1 (MQ=255)
                                                                                                                           tgGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAg                    <  1:148581/148‑1 (MQ=255)
                                                                                                                                    gctgctGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACcctcc          >  5:324710/1‑149 (MQ=255)
                                                                                                                                    gctgctGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCGCGCCGACAAGGGCACcctcc          >  4:226086/1‑149 (MQ=255)
                                                                                                                                      tgctgcGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCATcctcct         >  6:68976/1‑148 (MQ=255)
                                                                                                                                        ctgcGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCATcctcct         <  5:68976/146‑1 (MQ=255)
                                                                                                                                        ctgcGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGt      <  2:334091/149‑1 (MQ=255)
                                                                                                                                           cgcCCGGTGTAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGTTcc   <  1:133483/149‑1 (MQ=255)
                                                                                                                                            gcCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGTTCCt  <  7:360974/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCTTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGTTCCT  >  NC_002947/1982474‑1982762

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: