Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,158,869 C→G 5.4% T306S (ACT→AGT)  PP_4542 → ABC transporter ATP‑binding protein/permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,158,8690CG 5.4% 285.9 / 4.2 90T306S (ACT→AGT) PP_4542ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base C (47/38);  new base G (2/3);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CGAGAAAACCCAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATGCGCCGG  >  NC_002947/5158720‑5159012
                                                                                                                                                     |                                                                                                                                               
cGAGAAAACCCAGCTGGCCTCCCGGGTGATCACCAGCATGGCCGAGCACAACAACACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCATAGCACTGGGACGGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCAc                                                                                                                                                 >  1:1020347/1‑150 (MQ=255)
      aaCCCAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACAACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCt                                                                                                                                           >  1:88348/1‑150 (MQ=255)
       aCCCAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTg                                                                                                                                          <  1:1723246/150‑1 (MQ=255)
       aCCCAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTg                                                                                                                                          >  1:1185976/1‑150 (MQ=255)
       aCCCAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCAACTGGCCTg                                                                                                                                          <  2:857948/150‑1 (MQ=255)
        cccAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTgg                                                                                                                                         >  1:65599/1‑150 (MQ=255)
         ccAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGt                                                                                                                                        >  2:1692430/1‑150 (MQ=255)
         ccAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGt                                                                                                                                        <  2:891034/150‑1 (MQ=255)
           aGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGa                                                                                                                                      <  2:529706/150‑1 (MQ=255)
                   tCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATc                                                                                                                              >  2:1268181/1‑150 (MQ=255)
                     ccGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACAGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGt                                                                                                                            >  2:304117/1‑150 (MQ=255)
                        ggTGATCACCAGCATGGACGTCGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCaa                                                                                                                         <  2:1699785/150‑1 (MQ=255)
                             tCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGt                                                                                                                    >  1:147488/1‑150 (MQ=255)
                             tCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGt                                                                                                                    >  2:890750/1‑150 (MQ=255)
                                      tGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGa                                                                                                           <  1:189097/150‑1 (MQ=255)
                                      tGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGa                                                                                                           <  2:1580142/150‑1 (MQ=255)
                                      tGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGa                                                                                                           <  2:130818/150‑1 (MQ=255)
                                        gACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGGCCCGCCTCGTCAACATGTCGGGCTGGATc                                                                                                         >  2:1902538/1‑150 (MQ=255)
                                        gACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATc                                                                                                         <  1:1268181/150‑1 (MQ=255)
                                         aCGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCa                                                                                                        >  1:1112790/1‑150 (MQ=255)
                                         aCGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCa                                                                                                        >  2:1016851/1‑150 (MQ=255)
                                           gTGGTCATCACCACCCTCAATGGCCGGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGGCAATCGCTGATCACCGTGGGTGCGATCGCCCTGGCCACTGGCCTGGTGTTCCGCATCGTCCACAGGTCGGGCTgcgtgctg                                                                                                      >  2:529774/1‑143 (MQ=255)
                                               tcatcaCCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGAGCCGCATCGTCAACAGGTCGGTCTGGATAATgtcgg                                                                                                  >  2:892601/1‑150 (MQ=255)
                                                catcaCCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATgtgggt                                                                                                 >  2:224217/1‑150 (MQ=255)
                                                 atcaCCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCAGGCCAGTGGccaggg                                                                                                                                        >  1:833543/1‑106 (MQ=255)
                                                 atcaCCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCAGGCCAGTGGCCAggg                                                                                                                                        <  2:833543/110‑2 (MQ=255)
                                                  tcaCCACCCTCAATGGCCGGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGATATCGACTTGGATGGCTTCGGTGTGGTTCTCATCGTGAACATGTCGGGCTGGATCATgtgggtgg                                                                                               <  1:529774/150‑1 (MQ=255)
                                                     ccaccCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCa                                                                                            >  1:1461923/1‑150 (MQ=255)
                                                      caccCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCaa                                                                                           >  2:1532208/1‑150 (MQ=255)
                                                        cccTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACg                                                                                         <  1:304117/150‑1 (MQ=255)
                                                        cccTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACg                                                                                         <  2:147488/150‑1 (MQ=255)
                                                        cccTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGACTTTTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACg                                                                                         <  1:542428/150‑1 (MQ=255)
                                                          cTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGc                                                                                       >  2:1812635/1‑150 (MQ=255)
                                                          cTCAATGGCCTGCTGGTGGTCGCCACAACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGc                                                                                       <  2:88348/150‑1 (MQ=255)
                                                              aTGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATc                                                                                                         >  2:522260/1‑128 (MQ=255)
                                                              aTGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATc                                                                                                         <  1:522260/128‑1 (MQ=255)
                                                                  ccTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCga                                                                               <  1:1532208/150‑1 (MQ=255)
                                                                    tgctgGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCgaga                                                                             >  1:581478/1‑150 (MQ=255)
                                                                      ctgGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAAc                                                                           >  2:569232/1‑150 (MQ=255)
                                                                          tggtCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCg                                                                       <  1:1812635/150‑1 (MQ=255)
                                                                            gtCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGc                                                                     <  1:1016851/150‑1 (MQ=255)
                                                                             tCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACa                                                                          <  2:77148/144‑1 (MQ=255)
                                                                             tCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACa                                                                          >  1:77148/1‑144 (MQ=255)
                                                                               gCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCTTCGCCCTGGCCACTCGCTTTTTGCTCCGCCTGGGCCAGATTTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATg                                                                  <  1:1692430/150‑1 (MQ=255)
                                                                                  accaccGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTg                                                               <  1:858251/150‑1 (MQ=255)
                                                                                    caccGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCa                                                             >  1:1314867/1‑150 (MQ=255)
                                                                                     accGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAg                                                            >  2:1745965/1‑150 (MQ=255)
                                                                                            ttGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTc                                                                                                                   >  1:1271802/1‑88 (MQ=255)
                                                                                            ttGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTc                                                                                                                   <  2:1271802/88‑1 (MQ=255)
                                                                                                ccTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCa                                                 >  2:1644552/1‑150 (MQ=255)
                                                                                                  tgtggctgtggAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGa                                               <  1:488284/150‑1 (MQ=255)
                                                                                                       ctgtggAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTc                                                                                 <  2:931425/111‑1 (MQ=255)
                                                                                                       ctgtggAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTc                                                                                 >  1:931425/1‑111 (MQ=255)
                                                                                                        tgtggGGCCAATCGCTGAGCGCCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCg                                         <  2:440865/150‑1 (MQ=255)
                                                                                                               ccAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGaa                                                                            <  2:461280/108‑1 (MQ=255)
                                                                                                               ccAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGaa                                                                            >  1:461280/1‑108 (MQ=255)
                                                                                                                 aaTCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGGATGGTGAAGGACGGCCTGCAGACCATCGCCCAGCCgg                                >  1:1773618/1‑150 (MQ=255)
                                                                                                                      cTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGt                                                                >  1:1903573/1‑113 (MQ=255)
                                                                                                                      cTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGt                                                                <  2:1903573/113‑1 (MQ=255)
                                                                                                                       tGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCCGCCGgtcaccc                          >  2:787964/1‑149 (MQ=255)
                                                                                                                       tGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGgtcaccg                          >  1:1934837/1‑150 (MQ=255)
                                                                                                                          tCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCg                                                                       >  2:434495/1‑102 (MQ=255)
                                                                                                                          tCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGgtcaccgtca                       <  2:581478/150‑1 (MQ=255)
                                                                                                                          gcACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCg                                                                       <  1:434495/101‑1 (MQ=255)
                                                                                                                           cACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGgtcaccgtcac                      <  1:892601/150‑1 (MQ=255)
                                                                                                                             ccGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCAtg                                                                                                      >  1:307437/1‑68 (MQ=255)
                                                                                                                             ccGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCAtg                                                                                                      <  2:307437/68‑1 (MQ=255)
                                                                                                                             ccGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGtcaccgtcaccg                    >  2:1858597/1‑150 (MQ=255)
                                                                                                                               gTGGGTGCCATCGCCCTGGCCACTGGCCTTGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACCGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGAc                  <  2:368359/150‑1 (MQ=255)
                                                                                                                                   gTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCgag                                                                              <  2:903281/86‑1 (MQ=255)
                                                                                                                                   gTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCgag                                                                              >  1:903281/1‑86 (MQ=255)
                                                                                                                                   gTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGAc                                              <  2:970686/118‑1 (MQ=255)
                                                                                                                                   gTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGAc                                              >  1:970686/1‑118 (MQ=255)
                                                                                                                                    tGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCa                                                             <  2:799920/102‑1 (MQ=255)
                                                                                                                                    tGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCa                                                             >  1:799920/1‑102 (MQ=255)
                                                                                                                                     gCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGccc            >  1:834709/1‑150 (MQ=255)
                                                                                                                                     gCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGccc            <  1:224217/150‑1 (MQ=255)
                                                                                                                                      ccATCGCCCTGGCCACTGGCCTGGTGATCCGCGCTGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCCTCTTCGGGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCa           >  1:155600/1‑150 (MQ=255)
                                                                                                                                      ccATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCa           >  1:1634494/1‑150 (MQ=255)
                                                                                                                                       cATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGTGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCaa          >  1:1178875/1‑150 (MQ=255)
                                                                                                                                       cATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCaa          <  1:890750/150‑1 (MQ=255)
                                                                                                                                        aTCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAAt         >  1:1892233/1‑150 (MQ=255)
                                                                                                                                           gCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATgcg      >  1:384152/1‑150 (MQ=255)
                                                                                                                                            cccTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATgcgc     >  1:635511/1‑150 (MQ=255)
                                                                                                                                             cgtGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTGGGGCTAGAGCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATGCGcc    <  1:110241/148‑1 (MQ=255)
                                                                                                                                             ccTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATGCGcc    >  1:1940535/1‑150 (MQ=255)
                                                                                                                                              caggCCAGTGGCCAGGTGATCCGCATCGTCAACATGTCGGGCTGGATCAtg                                                                                                      <  2:119768/49‑1 (MQ=255)
                                                                                                                                              caggCCAGTGGCCAGGTGATCCGCATCGTCAACATGTCGGGCTGGATCAtg                                                                                                      >  1:119768/3‑51 (MQ=255)
                                                                                                                                               tGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATGCGCCgg  >  1:1603391/1‑150 (MQ=255)
                                                                                                                                                  ccACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCa                                                             >  1:790337/1‑88 (MQ=255)
                                                                                                                                                  ccACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCa                                                             <  2:790337/88‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGAGAAAACCCAGCTGGCCTCCCGGGTGATCACCAGCATGGACGTGGTCATCACCACCCTCAATGGCCTGCTGGTGGTCGCCACCACCGGCCTTGCCCTGTGGCTGTGGAGCCAATCGCTGATCACCGTGGGTGCCATCGCCCTGGCCACTGGCCTGGTGATCCGCATCGTCAACATGTCGGGCTGGATCATGTGGGTGGTCAACGGCATCTTCGAGAACATCGGCATGGTGCAGGACGGCCTGCAGACCATCGCCCAGCCGGTCACCGTCACCGACCCGCCCAATGCGCCGG  >  NC_002947/5158720‑5159012

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