Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,947 T→C 7.8% H63H (CAT→CAC PP_1239 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,416,9470TC 7.8% 182.1 / 4.8 64H63H (CAT→CACPP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (2/3);  total (30/34)
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.37e-01

CACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAGCG  >  NC_002947/1416800‑1417094
                                                                                                                                                   |                                                                                                                                                   
cACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGTGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATTCCGACCACGTGCATgg                                                                                                                                                   <  2:530677/150‑1 (MQ=255)
 aCGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGCCCCAGTCGAACATGCCGACCACGTGCATggg                                                                                                                                                  <  1:1107387/150‑1 (MQ=255)
   gCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGGCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGt                                                                                                                                                <  1:385272/150‑1 (MQ=255)
   gCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGt                                                                                                                                                <  2:1003728/150‑1 (MQ=255)
     tGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCGCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCg                                                                                                                                              <  1:629351/150‑1 (MQ=255)
      gATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCgg                                                                                                                                             >  2:48503/1‑150 (MQ=255)
       aTCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGc                                                                                                                                            <  2:608316/150‑1 (MQ=255)
       aTCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGCACACCCACGAACATGCCGACCACGTGCATGGGGTCGGc                                                                                                                                            >  1:784584/1‑150 (MQ=255)
        tCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGca                                                                                                                                           <  2:145696/150‑2 (MQ=255)
          gCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTg                                                                                                                                         >  1:412475/1‑150 (MQ=255)
            cAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGGTGGTCACCCACGAACATGCAGACCACGAGCATGGGTGCGGCTTGCt                                                                                                                                       >  2:1083401/1‑150 (MQ=255)
                 gtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGcat                                                                                                                                          <  1:1417271/142‑3 (MQ=255)
                 gtgACACGTTAATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGcat                                                                                                                                          >  2:1417271/1‑140 (MQ=255)
                    acacGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgc                                                                                                                               <  1:48503/150‑1 (MQ=255)
                    acacGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgc                                                                                                                               <  2:617750/150‑1 (MQ=255)
                    acacGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgc                                                                                                                               >  2:539692/1‑150 (MQ=255)
                        gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTa                                                                                                                           >  1:351351/1‑150 (MQ=255)
                        gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTa                                                                                                                           >  2:869317/1‑150 (MQ=255)
                         ttCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTAc                                                                                                                          <  2:1164606/150‑1 (MQ=255)
                          tCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACa                                                                                                                         <  1:563786/150‑1 (MQ=255)
                            aTCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAAt                                                                                                                       <  1:757148/150‑1 (MQ=255)
                               cTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTa                                                                                                                    <  2:271578/150‑1 (MQ=255)
                                        tGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGTTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTAc                                                                                                           <  2:412475/150‑1 (MQ=255)
                                           ggCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGTAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTc                                                                                                        >  1:1380015/1‑150 (MQ=255)
                                                   ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAgg                                                                                                >  2:507217/1‑150 (MQ=255)
                                                    cTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAggg                                                                                               >  2:1129189/1‑150 (MQ=255)
                                                     tGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGa                                                                                              >  2:1143888/1‑150 (MQ=255)
                                                      gCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGAc                                                                                             <  2:1029577/150‑1 (MQ=255)
                                                              ccGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGcgcg                                                                                     >  2:840247/1‑150 (MQ=255)
                                                              ccGAGCAGCGCCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACgg                                                                                                                                <  1:756385/107‑1 (MQ=255)
                                                              ccGAGCAGCGCCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACgg                                                                                                                                >  2:756385/1‑107 (MQ=255)
                                                                    ggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGCCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATgc                                                                               >  2:754288/1‑150 (MQ=255)
                                                                                       ggTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGc                                                            <  2:686119/150‑1 (MQ=255)
                                                                                       ggTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGc                                                            <  1:1083401/150‑1 (MQ=255)
                                                                                        gTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGcc                                                           >  2:1516/1‑150 (MQ=255)
                                                                                         tttCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCg                                                          <  1:840247/150‑1 (MQ=255)
                                                                                          ttCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCTGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCgg                                                         <  2:784584/150‑1 (MQ=255)
                                                                                             tGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGa                                                                                              >  2:579894/1‑110 (MQ=255)
                                                                                             tGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGa                                                                                              <  1:579894/110‑1 (MQ=255)
                                                                                               cagctcagctCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTCCAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTtc                                                    <  1:810479/150‑1 (MQ=255)
                                                                                                  ctcagctcagTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCCACCGCAGCCAAGTGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTcg                                                 <  2:349228/150‑1 (MQ=255)
                                                                                                   tcagctcagTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTcgc                                                <  2:351351/150‑1 (MQ=255)
                                                                                                        tcagTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGtgtg                                           >  2:650343/1‑150 (MQ=255)
                                                                                                             gcagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGtt                                                       <  2:1153966/133‑1 (MQ=255)
                                                                                                             gcagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGtt                                                       >  1:1153966/1‑133 (MQ=255)
                                                                                                              cagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCaaa                                     <  1:1516/150‑1 (MQ=255)
                                                                                                                 tgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGcc                                  <  1:539692/150‑1 (MQ=255)
                                                                                                                      gTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCat                             >  2:1398301/1‑150 (MQ=255)
                                                                                                                           ccACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCaaa                                     <  2:164230/137‑1 (MQ=255)
                                                                                                                           ccACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCaaa                                     >  1:164230/1‑137 (MQ=255)
                                                                                                                           ccACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCgg                        >  1:1080697/1‑150 (MQ=255)
                                                                                                                            cACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACgg                                                                                                                                <  2:162873/45‑1 (MQ=255)
                                                                                                                            cACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACgg                                                                                                                                >  1:162873/1‑45 (MQ=255)
                                                                                                                            cACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGc                       >  2:366204/1‑150 (MQ=255)
                                                                                                                               gAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCt                                 >  1:1209741/1‑137 (MQ=255)
                                                                                                                               gAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCt                                 <  2:1209741/137‑1 (MQ=255)
                                                                                                                                 aCATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAGCCTGCATATCGGCAGCCTg                 >  1:1082882/1‑150 (MQ=255)
                                                                                                                                   aagCCGACCCCATGCACGTGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTgg                >  1:553310/3‑150 (MQ=255)
                                                                                                                                   aTGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTgg                >  2:386361/1‑150 (MQ=255)
                                                                                                                                     gCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGaa              <  1:507217/150‑1 (MQ=255)
                                                                                                                                     cccGACCCCATGCACGTGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGaa              <  2:553310/149‑1 (MQ=255)
                                                                                                                                         ccaaCGTGCATGGGGTCGGTTTGCTGTCCCGGCCCTACAATTTACCGGTTTACCTCAGCCAAGGGACTTTGCGGGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCCCGCGTGTGGCCAAGCCTGCATATCGGCAGCCTGGAAGTAAc         <  2:1082882/147‑1 (MQ=255)
                                                                                                                                              gTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCa     >  1:630206/1‑150 (MQ=255)
                                                                                                                                                 cATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAgcg  >  2:1248089/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAGCG  >  NC_002947/1416800‑1417094

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