Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 60,907 A→G 18.8% D26G (GAC→GGC)  PP_0052 → beta‑lactamase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294760,9070AG18.8% 143.9 / 26.2 69D26G (GAC→GGC) PP_0052beta‑lactamase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (6/7);  total (32/37)
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.94e-01

CATCTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCGGTGGGTGCTGGAAAC  >  NC_002947/60764‑61048
                                                                                                                                               |                                                                                                                                             
catCTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTg                                                                                                                                          <  5:47510/149‑1 (MQ=255)
 atCTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTg                                                                                                                                          >  6:180461/1‑148 (MQ=255)
   cTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCg                                                                                                                                            >  5:243971/1‑144 (MQ=255)
   cTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCg                                                                                                                                            <  6:243971/144‑1 (MQ=255)
   cTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGaaa                                                                                                                                       >  2:349741/1‑149 (MQ=255)
     tGATAAAGCTCCTGTCACTGCCTACGCCTGTAGAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCCCGTGGCCGCCTTCTACGCCGCGGCGACGTCGACGCTCAGCTACCTGGGCATGGACCGTGCAAcc                                                                                                                                     >  2:291736/1‑149 (MQ=255)
      gATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAAccc                                                                                                                                    >  6:94822/1‑149 (MQ=255)
           aGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGATGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACaa                                                                                                                               >  4:335968/1‑149 (MQ=255)
           aGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACaa                                                                                                                               >  4:124854/1‑149 (MQ=255)
           aGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACaa                                                                                                                               >  4:159693/1‑149 (MQ=255)
             cTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATg                                                                                                                             >  7:216847/1‑149 (MQ=255)
                cTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATgcg                                                                                                                           >  7:69586/1‑148 (MQ=255)
                   tCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCAt                                                                                                                        <  3:279527/148‑1 (MQ=255)
                         ccTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATc                                                                                                                   <  8:216847/147‑1 (MQ=255)
                               ccTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGt                                                                                                           <  3:124854/149‑1 (MQ=255)
                                cTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTg                                                                                                          >  6:147524/1‑149 (MQ=255)
                                   tACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCt                                                                                                        <  7:276619/148‑1 (MQ=255)
                                     cAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGa                                                                                                     <  4:180975/149‑1 (MQ=255)
                                      aaGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAAccc                                                                                                                                    >  2:259070/1‑117 (MQ=255)
                                      aaGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAAccc                                                                                                                                    <  1:259070/117‑1 (MQ=255)
                                      aaGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGAt                                                                                                    <  3:298022/149‑1 (MQ=255)
                                       aGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGATGCCTTCTACGACGAGGCGACCTCGACCATCAGGTACCTGGTCATGGACCGTGAAACGCGACAATTCGCATTGATCGACAGCGTGCTGGAtt                                                                                                   <  3:335968/149‑1 (MQ=255)
                                        ggcggggcggATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTAACTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTa                                                                                                  <  5:222029/149‑1 (MQ=255)
                                             ggcggATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGAcc                                                                                             <  7:273159/149‑1 (MQ=255)
                                                 gATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCaa                                                                                         <  6:40660/149‑1 (MQ=255)
                                                   tAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGt                                                                                       >  2:62126/1‑149 (MQ=255)
                                                            gaACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCa                                                                              >  3:337430/1‑149 (MQ=255)
                                                                gTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTg                                                                          >  4:224942/1‑149 (MQ=255)
                                                                     gatCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCAACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAgcg                                                                     >  8:290225/1‑149 (MQ=255)
                                                                         atcGTCAACACCCTTCACGGGGACGCCTGCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGCACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGTATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTc                                                                 <  5:8767/149‑1 (MQ=255)
                                                                         atcGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTc                                                                 <  5:94822/149‑1 (MQ=255)
                                                                          tcGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCg                                                                >  5:296558/1‑149 (MQ=255)
                                                                          tcGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCg                                                                <  5:180461/149‑1 (MQ=255)
                                                                            ggCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCgg                                                               <  4:337430/148‑1 (MQ=255)
                                                                               aacaacCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGa                                                           >  5:203890/1‑149 (MQ=255)
                                                                                acaacCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGAt                                                          <  6:287777/149‑1 (MQ=255)
                                                                                    ccTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGcc                                                      <  8:69586/149‑1 (MQ=255)
                                                                                     cTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCt                                                     >  3:338049/1‑149 (MQ=255)
                                                                                        cACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTg                                                                                                          <  6:37896/93‑1 (MQ=255)
                                                                                        cACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTg                                                                                                          >  5:37896/1‑93 (MQ=255)
                                                                                           gTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGa                                               <  3:224942/149‑1 (MQ=255)
                                                                                                 gCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCgtgt                                         <  4:338049/149‑1 (MQ=255)
                                                                                                   cTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGtgtg                                        >  1:191848/1‑148 (MQ=255)
                                                                                                      cTACGTCGAGGCGTCCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACACTGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAAc                                     <  4:243171/148‑1 (MQ=255)
                                                                                                       tACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGa                                   >  4:270512/1‑149 (MQ=255)
                                                                                                         cgacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGc                                 <  6:296558/149‑1 (MQ=255)
                                                                                                         cgacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCTCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAg                                  <  1:291736/148‑1 (MQ=255)
                                                                                                         cgacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGa                                   <  1:349741/147‑1 (MQ=255)
                                                                                                          gacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGa                                       >  5:70897/1‑142 (MQ=255)
                                                                                                          gacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGa                                       <  6:70897/142‑1 (MQ=255)
                                                                                                          gacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCt                                >  5:89110/1‑149 (MQ=255)
                                                                                                          gacgaGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCt                                >  6:44714/1‑149 (MQ=255)
                                                                                                              aGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAAt                            <  6:89110/149‑1 (MQ=255)
                                                                                                                 cGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCgtgt                                         >  6:353906/1‑133 (MQ=255)
                                                                                                                 cGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCgtgt                                         <  5:353906/133‑1 (MQ=255)
                                                                                                                     cTCGACCATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAgtg                     >  4:123195/1‑149 (MQ=255)
                                                                                                                        gACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGc                  >  3:27812/1‑149 (MQ=255)
                                                                                                                         aCCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCg                 <  5:44714/149‑1 (MQ=255)
                                                                                                                          ccATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCgg                >  4:96825/1‑149 (MQ=255)
                                                                                                                           cATCAGCTACCTGGTCATGGGCCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCggt               <  1:62126/149‑1 (MQ=255)
                                                                                                                             tCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCggtgg             <  2:191848/149‑1 (MQ=255)
                                                                                                                               aGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCggtggg            <  5:147524/148‑1 (MQ=255)
                                                                                                                                gCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCggtgggtg          >  3:62732/1‑149 (MQ=255)
                                                                                                                                    ccTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCGGTGGGTGCt        <  3:270512/147‑1 (MQ=255)
                                                                                                                                     cTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCGGTGGGTGCTgg      >  6:276890/1‑148 (MQ=255)
                                                                                                                                        gTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCGGTGGGTGCTGGaaa   <  3:159693/148‑1 (MQ=255)
                                                                                                                                        gTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCGGTGGGTGCTGGAAAc  <  7:290225/149‑1 (MQ=255)
                                                                                                                                         tCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGc                          >  6:92094/1‑124 (MQ=255)
                                                                                                                                         tCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGc                          <  5:92094/124‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CATCTTGATAAAGCTCCTGTCACTGCCTACGCCTGTACAAGGCGGGGCGGATAAATTGGAGAACGTGATGATCATCGGCAACAACCTTCACGTGGACGCCTTCTACGACGAGGCGACCTCGACCATCAGCTACCTGGTCATGGACCGTGAAACCCGACAATGCGCATTGATCGACAGCGTGCTGGATTATGACCCCAAGTCCGGGCGCACCTGCAGCGCTTCGGCTGATCGCCTGATCGAGCGTGTGAACGAGCTGAATGCCAGTGTGCGGTGGGTGCTGGAAAC  >  NC_002947/60764‑61048

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: