Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,582,067 C→A 5.4% L351M (CTG→ATG)  PP_4910 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,582,0670CA 5.4% 295.2 / 4.4 92L351M (CTG→ATG) PP_4910hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/46);  new base A (3/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

TGCTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTT  >  NC_002947/5581918‑5582215
                                                                                                                                                     |                                                                                                                                                    
tgCTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCAccc                                                                                                                                                      >  2:1675823/1‑150 (MQ=255)
tgCTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCAccc                                                                                                                                                      >  2:1676200/1‑150 (MQ=255)
  cTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTg                                                                                                                                                    <  1:1667081/150‑1 (MQ=255)
   ttGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTgg                                                                                                                                                   >  1:662538/1‑150 (MQ=255)
   ttGCCGTGCAGATGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACAAGCCTGGCCGAGCTGACCGCGGGAGCGTATGGTGGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTgg                                                                                                                                                   <  2:1004858/150‑1 (MQ=255)
    tGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGt                                                                                                                                                  >  1:1544144/1‑150 (MQ=255)
     gCCGTGCAGCTGGGCTCGCCGGAACAACTGCGTAGCCTGCCATGCGCAAGCCTGGCCGAGCTTACCGATCGCGACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCCCAAGCTGAACGACGGGAGCACCCTGGTg                                                                                                                                                 <  2:1366803/150‑1 (MQ=255)
     gCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTg                                                                                                                                                 <  1:166345/150‑1 (MQ=255)
      ccGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGGCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGc                                                                                                                                                >  2:783208/1‑150 (MQ=255)
      ccGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGc                                                                                                                                                >  2:1083588/1‑150 (MQ=255)
        gTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGNCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACTGCAGCACCCTGGTGCTg                                                                                                                                              >  1:1526978/1‑150 (MQ=255)
        gTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTg                                                                                                                                              >  1:1537477/1‑150 (MQ=255)
         tGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTTACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGt                                                                                                                                             <  1:860234/150‑1 (MQ=255)
         tGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGt                                                                                                                                             <  1:645271/150‑1 (MQ=255)
           cAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTc                                                                                                                                           <  1:1393877/150‑1 (MQ=255)
               tGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTTGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCcaca                                                                                                                                       <  1:752938/150‑1 (MQ=255)
               tGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCcaca                                                                                                                                       >  2:238347/1‑150 (MQ=255)
                      gCAGGCGCTACTGCGGCGGCTGCCGCGCACACGCCTGGCCGAGCTGACCGAGCGGTAGCGTTATTGCCCGGAGCACGCGGCCTCGCGCTACGACGGCGAACGGCACACGCTGATCGTAGGCAGCACCCTTGTGGTGTTCCACAGCCGCGa                                                                                                                                <  2:693592/150‑1 (MQ=255)
                      gCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTc                                                                                                                                           >  1:1455996/1‑139 (MQ=255)
                      gCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTc                                                                                                                                           <  2:1455996/139‑1 (MQ=255)
                      gCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                >  1:1161570/1‑150 (MQ=255)
                      cgatctgCAACTGCGTCGCCTGCCACGCACTCGCCCGTCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCGTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                <  2:1095458/144‑1 (MQ=255)
                       cAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAg                                                                                                                               <  2:1544144/150‑1 (MQ=255)
                        aGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAgt                                                                                                                              <  1:1911649/150‑1 (MQ=255)
                          gAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAgtgt                                                                                                                            >  1:278327/1‑150 (MQ=255)
                             cAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGc                                                                                                                         >  2:642768/1‑150 (MQ=255)
                             cAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGc                                                                                                                         <  2:1526978/150‑1 (MQ=255)
                             cAACTGCGTCGCCTGCCACGCACACGCCTGGCACAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAgagatcg                                                                                                                         >  1:1095458/1‑144 (MQ=255)
                              aaCTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGtgtg                                                                                                                           <  2:1541677/147‑1 (MQ=255)
                              aaCTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGtgtg                                                                                                                           >  1:1541677/1‑147 (MQ=255)
                               aCTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGa                                                                                                                       >  1:1559141/1‑150 (MQ=255)
                                cTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      <  2:897086/150‑1 (MQ=255)
                                cTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      <  1:1432364/150‑1 (MQ=255)
                                cTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      <  1:476368/150‑1 (MQ=255)
                                  gCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGa                                                                                                                    <  2:1765213/150‑1 (MQ=255)
                                  gCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGa                                                                                                                    <  2:491864/150‑1 (MQ=255)
                                   cGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                <  2:502219/137‑1 (MQ=255)
                                   cGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                >  1:502219/1‑137 (MQ=255)
                                     tCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTa                                                                                                                 <  2:1783069/150‑1 (MQ=255)
                                        ccTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGt                                                                                                                                                  <  2:1235659/114‑1 (MQ=255)
                                        ccTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGt                                                                                                                                                  >  1:1235659/1‑114 (MQ=255)
                                        aggcccACGCACACACCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGGTGTTACACAGACGCGAGTGTGGCGAGGACTACCt                                                                                                              >  1:1005684/5‑150 (MQ=255)
                                            ccAGGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          >  2:1508983/1‑150 (MQ=255)
                                            ccACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          >  2:1222866/1‑150 (MQ=255)
                                            ccACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          >  2:1302952/1‑150 (MQ=255)
                                            ccACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          >  1:1047477/1‑150 (MQ=255)
                                            ccACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          >  2:552444/1‑150 (MQ=255)
                                             cACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCa                                                                                                         >  1:1257631/1‑150 (MQ=255)
                                                 cacaCGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGATAGGCACCCGCTGAACGTCGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGc                                                                                                     >  1:609407/1‑150 (MQ=255)
                                                  acacGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAg                                                                                                                                      >  1:807975/1‑116 (MQ=255)
                                                  acacGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAg                                                                                                                                      <  2:807975/116‑1 (MQ=255)
                                                          ggCCGAGCTGAACGAGCGCCACCGAGATTGCCTGGAGCACGCGGACCCGCTCTACGACGGCGAAACGCACACGCTGAACGTCGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCCTTTGAtg                                                                                            <  1:274695/150‑1 (MQ=255)
                                                           gCCGAGCTTACCGAGCGCTACCGTGATTGACTGGAGCACGCTGCCTCGCTCTACGACGGCTAAACTCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACACCCGCGAGTTTGGCGAGGACTACCTGACCACCGCCATTTGCtgt                                                                                           <  2:210754/150‑1 (MQ=255)
                                                           gCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCtgt                                                                                           <  2:444615/150‑1 (MQ=255)
                                                           gCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCtgt                                                                                           <  1:325679/150‑1 (MQ=255)
                                                             cGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTgg                                                                                         <  2:99265/150‑1 (MQ=255)
                                                                 cTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGaa                                                                                     <  1:1222866/150‑1 (MQ=255)
                                                                    aCCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCcaca                                                                                                                                       >  2:1245782/1‑97 (MQ=255)
                                                                    aCCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCcaca                                                                                                                                       <  1:1245782/97‑1 (MQ=255)
                                                                      cGAGCGCTACCGGGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGtt                                                                                <  1:65237/150‑1 (MQ=255)
                                                                           gcTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCCTTTCCTGTGGCCAACTAATgcgcg                                                                           >  2:870020/1‑150 (MQ=255)
                                                                            cTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      >  1:1709675/1‑106 (MQ=255)
                                                                            cTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      <  2:1709675/106‑1 (MQ=255)
                                                                             tACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGcc                                                                         >  2:1785621/1‑150 (MQ=255)
                                                                               ccGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCt                                                                       >  1:859336/1‑150 (MQ=255)
                                                                                 gTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTgg                                                                     <  1:238347/150‑1 (MQ=255)
                                                                                     ttGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCa                                                                 >  2:1135313/1‑150 (MQ=255)
                                                                                       gCCTGGAGCACGCGGCCTCGCTCTACGACGACGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTgcgc                                                                            >  1:8741/1‑137 (MQ=255)
                                                                                       gCCTGGAGCACGCGGCCTCGCTCTACGACGACGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTgcgc                                                                            <  2:8741/137‑1 (MQ=255)
                                                                                          tGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          >  1:1805507/1‑104 (MQ=255)
                                                                                          tGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGAcc                                                                                                          <  2:1805507/104‑1 (MQ=255)
                                                                                               cACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAAc                                                                                                       <  2:1403457/102‑1 (MQ=255)
                                                                                               cACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAAc                                                                                                       >  1:1403457/1‑102 (MQ=255)
                                                                                               cACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAg                                                       >  1:1438433/1‑150 (MQ=255)
                                                                                                     gCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAg                                                 <  1:552444/150‑1 (MQ=255)
                                                                                                     gCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAg                                                 <  1:1508983/150‑1 (MQ=255)
                                                                                                          gCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTAGCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGc                                            <  2:1093136/150‑1 (MQ=255)
                                                                                                          gCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGc                                            <  2:278327/150‑1 (MQ=255)
                                                                                                           ctctACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCt                                           <  1:1855212/150‑1 (MQ=255)
                                                                                                                cgacgGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAg                                      >  2:834606/1‑150 (MQ=255)
                                                                                                                         aaCGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTAc                             >  1:1684315/1‑150 (MQ=255)
                                                                                                                         aaCGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTAc                             >  1:1684314/1‑150 (MQ=255)
                                                                                                                          aCGCACACGCTGAACGACGGCAGCACCATGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTAcc                            >  1:356652/1‑150 (MQ=255)
                                                                                                                           cGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTAccc                           <  1:1785621/150‑1 (MQ=255)
                                                                                                                           cGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTAccc                           <  1:765784/150‑1 (MQ=255)
                                                                                                                                  gCTGAACGACGGCAGCACCCTGGTGCTGTTCCACGGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGCTCGAGGTAGCTGACAGCGGCATTACCCTGCAAcg                    >  2:125696/1‑149 (MQ=255)
                                                                                                                                          acgGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAg            <  2:1537477/150‑1 (MQ=255)
                                                                                                                                           gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCaa                      >  2:1137071/4‑139 (MQ=255)
                                                                                                                                           gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCaa                      <  1:1137071/136‑1 (MQ=255)
                                                                                                                                            aacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTg                                                                                  <  1:118617/76‑1 (MQ=255)
                                                                                                                                            aacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTg                                                                                  >  2:118617/3‑78 (MQ=255)
                                                                                                                                              cagcaCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCt        >  1:909047/1‑150 (MQ=255)
                                                                                                                                              cagcaCCCTGGTGCTGGTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGTAACTGCAGATAGGCCt        >  2:974475/1‑150 (MQ=255)
                                                                                                                                               agcaCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTg       <  1:834606/150‑1 (MQ=255)
                                                                                                                                                    ccTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGCCTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCtt  >  2:404188/1‑150 (MQ=255)
                                                                                                                                                    ccTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTAACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCtt  >  1:1739102/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTT  >  NC_002947/5581918‑5582215

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: