Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,497,841 T→G 8.0% M141R (ATG→AGG)  PP_1310 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,497,8410TG 8.0% 230.5 / 12.5 88M141R (ATG→AGG) PP_1310hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/41);  new base G (4/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

TGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCG  >  NC_002947/1497692‑1497990
                                                                                                                                                     |                                                                                                                                                     
tGCAGCATGGCTGCGTGCCAGGCACCCCTGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAt                                                                                                                                                       <  1:1588082/150‑1 (MQ=255)
  cagcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCAGGATCGCCCCGGCGGTCATgc                                                                                                                                                     <  2:302984/150‑1 (MQ=255)
  cagcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGAGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATgc                                                                                                                                                     <  1:1674239/150‑1 (MQ=255)
    gcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACATGTTCATGATCGCCCCGGCGGTCATgcgc                                                                                                                                                   <  1:1842122/150‑1 (MQ=255)
     caTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTTGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGgcgca                                                                                                                                                  >  1:1796845/1‑149 (MQ=255)
     caTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATgc                                                                                                                                                     >  2:1622199/1‑147 (MQ=255)
     caTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATgc                                                                                                                                                     <  1:1622199/147‑1 (MQ=255)
      aTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCt                                                                                                                                                 <  2:1578379/150‑1 (MQ=255)
           tGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTAGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAccg                                                                                                                                            >  2:861873/1‑150 (MQ=255)
             cGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAccgcc                                                                                                                                          <  1:1092943/150‑1 (MQ=255)
               tGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTTGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgccgg                                                                                                                                        <  2:1796845/150‑2 (MQ=255)
                 ccAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTCACCGCCNAAt                                                                                                                                      >  1:1366392/1‑150 (MQ=255)
                  cAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAtt                                                                                                                                     >  2:1261940/1‑150 (MQ=255)
                  cAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAtt                                                                                                                                     <  2:1631066/150‑1 (MQ=255)
                  cAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGGCGGGTGTCTGCGCCGTGTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAtt                                                                                                                                     >  1:1377223/1‑150 (MQ=255)
                    ggCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTcc                                                                                                                                   <  1:1459331/150‑1 (MQ=255)
                    ggCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGACTGACCGCCGAATTcc                                                                                                                                   <  1:130251/150‑1 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgc                                                                                                                                           >  1:1648971/1‑137 (MQ=255)
                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTGCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgc                                                                                                                                           <  2:1648971/137‑1 (MQ=255)
                          ccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGc                                                                                                                             >  1:744018/1‑150 (MQ=255)
                          ccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGc                                                                                                                             >  2:1586810/1‑150 (MQ=255)
                           cGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGaa                                                                                                                                       >  2:675420/1‑139 (MQ=255)
                           cGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGaa                                                                                                                                       <  1:675420/139‑1 (MQ=255)
                           cGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGcc                                                                                                                            >  2:56326/1‑150 (MQ=255)
                           cGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCCGCGGTCATGAGCCTGACCGCCGAATTCCCGATGcc                                                                                                                            <  1:912254/150‑1 (MQ=255)
                                  ccGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTCTTCAAGATCGCTCGGGCGTTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAgcgc                                                                                                                     <  1:711253/150‑1 (MQ=255)
                                      tGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGGCCGCCGAATTCCCGATGCCGCAGCGCTGcc                                                                                                                 >  1:866662/1‑150 (MQ=255)
                                        gTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCg                                                                                                               >  1:501599/1‑150 (MQ=255)
                                           cTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTc                                                                                                            >  1:1631839/1‑150 (MQ=255)
                                                gggCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgt                                                                                                       >  1:1596128/1‑150 (MQ=255)
                                                 ggCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtg                                                                                                      <  2:1502317/150‑1 (MQ=255)
                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTCCCTTCTCAcctggg                                                                                                     >  2:49319/1‑146 (MQ=255)
                                                   cTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtggt                                                                                                    <  2:1041475/150‑1 (MQ=255)
                                                    ttttCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtggtg                                                                                                   >  1:1898459/1‑150 (MQ=255)
                                                            cgGTGGGTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGCTGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGa                                                                                           <  2:455967/150‑1 (MQ=255)
                                                             ggtggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGAc                                                                                          >  1:542976/1‑150 (MQ=255)
                                                              gtggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                  >  2:1449530/1‑89 (MQ=37)
                                                              gtggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTGCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                  <  1:1449530/93‑5 (MQ=37)
                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGAt                                                                                    <  2:1611608/150‑1 (MQ=255)
                                                                      gTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGTAGCACCTGTTCATGAACGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                                                 <  1:566055/150‑1 (MQ=255)
                                                                       ttGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAGTTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAg                                                                                >  1:1156617/1‑150 (MQ=255)
                                                                         gCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtc                                                                              <  1:769699/150‑1 (MQ=255)
                                                                           caccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAagtagt                                                                            >  1:634079/1‑150 (MQ=255)
                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                                           >  2:268630/1‑150 (MQ=255)
                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                                           >  1:26660/1‑150 (MQ=255)
                                                                              caGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCATCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTcgc                                                                         <  1:49319/150‑1 (MQ=255)
                                                                              caGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTAGTcgc                                                                         <  2:634079/150‑1 (MQ=255)
                                                                               aGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTcgcg                                                                        >  2:565824/1‑150 (MQ=255)
                                                                                    ggCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAGTTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCa                                                                   <  2:1156617/150‑1 (MQ=255)
                                                                                     gCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGcc                                                                                                                                                  >  2:1681075/1‑70 (MQ=255)
                                                                                     gCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGcc                                                                                                                                                  <  1:1681075/70‑1 (MQ=255)
                                                                                             ccTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTa                                                          <  1:186938/150‑1 (MQ=255)
                                                                                             ccTTGAAACTGCAGGCATCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTa                                                          >  2:973587/1‑150 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                         <  1:56326/150‑1 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGAGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCTTTGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                         <  2:219019/150‑1 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGgcgca                                                                                                                                                  <  2:230253/61‑2 (MQ=255)
                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGgcgca                                                                                                                                                  >  1:230253/1‑60 (MQ=255)
                                                                                              cTTGAAACTGAAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                         >  2:465869/1‑150 (MQ=255)
                                                                                               ttGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTTCCGATGCCGCAGCGCTGCGCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCCGCTGGTTTACg                                                        >  2:331059/1‑150 (MQ=255)
                                                                                                  aaaCTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                     <  1:1094136/150‑1 (MQ=255)
                                                                                                      tGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTc                                                 >  1:1107131/1‑150 (MQ=255)
                                                                                                      ggCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGCCCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCGGCTGGACGCCGATGAAGTCGTCGCGCCGCAGCTCGTTTACTAATGGTc                                                 <  2:268041/149‑1 (MQ=255)
                                                                                                         aGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGa                                              <  2:1631839/150‑1 (MQ=255)
                                                                                                           gCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACa                                            <  2:542976/150‑1 (MQ=255)
                                                                                                             gTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGc                                          >  2:1351153/1‑150 (MQ=255)
                                                                                                               cGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCt                                        <  1:1586810/150‑1 (MQ=255)
                                                                                                                  gTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGTGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTc                                     <  1:973587/150‑1 (MQ=255)
                                                                                                                   ttGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTcg                                    >  1:1753269/1‑150 (MQ=255)
                                                                                                                      cagcaCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGt                                 >  2:1310736/1‑150 (MQ=255)
                                                                                                                      cagcaCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGt                                 >  2:31445/1‑150 (MQ=255)
                                                                                                                         caCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTccc                              >  1:1722383/1‑150 (MQ=255)
                                                                                                                          aCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTcccc                             >  2:1439096/1‑150 (MQ=255)
                                                                                                                              gTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGAGTCCCCATGa                         <  2:26660/150‑1 (MQ=255)
                                                                                                                               ttCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                     >  1:72771/1‑121 (MQ=255)
                                                                                                                               ttCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAAt                                                     <  2:72771/121‑1 (MQ=255)
                                                                                                                                tCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgc                       <  2:432728/150‑1 (MQ=255)
                                                                                                                                  aTGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgctg                     >  2:604355/1‑150 (MQ=255)
                                                                                                                                   tGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAgctgc                    <  1:31445/150‑1 (MQ=255)
                                                                                                                                      tCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGTTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGa                 >  1:1475093/1‑150 (MQ=255)
                                                                                                                                        gCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtc                                                                              <  1:991346/87‑1 (MQ=255)
                                                                                                                                        gCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtc                                                                              >  2:991346/1‑87 (MQ=255)
                                                                                                                                         ccccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCg                                                                                      <  1:1923021/78‑1 (MQ=255)
                                                                                                                                         ccccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCg                                                                                      >  2:1923021/1‑78 (MQ=255)
                                                                                                                                         ccccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGt                                                             <  1:606232/103‑1 (MQ=255)
                                                                                                                                         ccccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGt                                                             >  2:606232/1‑103 (MQ=255)
                                                                                                                                           ccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTg            >  2:378590/1‑150 (MQ=255)
                                                                                                                                             ggcggTCAGGCGCATGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGCGGTAGTGGGGGGCCGACCGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGc          >  2:1483095/1‑150 (MQ=255)
                                                                                                                                              gcggTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         <  2:501599/150‑1 (MQ=255)
                                                                                                                                              gcggTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         <  2:202512/150‑1 (MQ=255)
                                                                                                                                                     tGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCg  <  1:1852773/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCG  >  NC_002947/1497692‑1497990

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: