Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,982,614 G→A 9.7% R189H (CGC→CAC)  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,6140GA 9.7% 172.2 / 4.9 61R189H (CGC→CAC) PP_1775metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (3/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TACCCTGGTCTTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGTTCCT  >  NC_002947/1982466‑1982762
                                                                                                                                                    |                                                                                                                                                     
tACCCTGGTCTTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTgcg                                                                                                                                                       >  8:170888/1‑149 (MQ=255)
   ccTGGTATTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCAccg                                                                                                                                                    >  6:329676/1‑148 (MQ=255)
        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                                                                                                                                >  2:340825/1‑149 (MQ=255)
                       tCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                >  3:332308/1‑149 (MQ=255)
                       tCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTg                                                                                                                                >  5:147935/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                                                                                                                                       <  2:76925/122‑6 (MQ=255)
                                           aTCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCACCGGGCGCAGcagcgg                                                                                                                                       >  1:76925/1‑117 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             <  4:369508/142‑1 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             >  3:279273/1‑142 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             <  4:279273/142‑1 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGtgctggtgctgg                                                                                                             >  3:369508/1‑142 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGaaa                                                                                                          <  7:170888/145‑1 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGaaa                                                                                                          <  6:349501/145‑1 (MQ=255)
                                                ctgGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGaaa                                                                                                          <  6:147935/145‑1 (MQ=255)
                                                           ttATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCgg                                                                                                                                      <  1:101792/107‑1 (MQ=255)
                                                           ttATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCgg                                                                                                                                      >  2:101792/1‑107 (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)
                                                                                                                                                    |                                                                                                                                                     
TACCCTGGTCTTGCCTGTCGTATTCGCCTGCAGCGTGCAGGGCATCTGCTGGGTTCCGCTTATGTCGAGTGCGATATGCAGCATGAGCAAACCAACAGCCGGTACGTGTTTTCTGGCGATCTGGGCGCTTGTGGCAACCCGCTGCTGCGCCCGGTGCAGCCCCCGGAGCGTGCTGACGTGCTGGTGCTGGAAAGCACCTACGGCGACCGCCTGCACCCGCAGGCAGACGACCGCCGTCAGCGGCTGGAAGACGCCATCGACCGAGCCCTCGCCGACAAGGGCACCCTCCTGGTTCCT  >  NC_002947/1982466‑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: