Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,906,508 A→T 8.7% L112Q (CTG→CAG)  PP_4315 ← epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,906,5080AT 8.7% 271.2 / 15.3 92L112Q (CTG→CAG) PP_4315epimerase
Reads supporting (aligned to +/- strand):  ref base A (38/46);  new base T (4/4);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

TCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGTG  >  NC_002947/4906359‑4906657
                                                                                                                                                     |                                                                                                                                                     
tCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCca                                                                                                                                                       <  2:2407389/150‑1 (MQ=255)
 cAACGTTTCGCGCAGCATGGCGTCATCCATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcag                                                                                                                                                      >  2:1835439/1‑150 (MQ=255)
  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                     >  2:1981487/1‑150 (MQ=255)
   aCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagca                                                                                                                                                    >  2:993859/1‑150 (MQ=255)
   aCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagca                                                                                                                                                    <  2:2042278/150‑1 (MQ=255)
   aCGCTTCGAGCAGCATGGCGTCATCGATTAGGACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCctgct                                                                                                                                                    >  1:675164/1‑149 (MQ=255)
    cGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagcag                                                                                                                                                   >  2:1339731/1‑150 (MQ=255)
        tCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCa                                                                                                                                               <  1:2184453/150‑1 (MQ=255)
         cGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAt                                                                                                                                              <  1:663449/150‑1 (MQ=255)
          gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtg                                                                                                                                             <  1:1372434/150‑2 (MQ=255)
            gcagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTCTGCCGCCCTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgga                                                                                                                                           <  2:675164/150‑4 (MQ=255)
            ccagcaTGGCGTCATCGATCACCACCACATAATCCTCGCTAAGGTACACCTCCTTCGCCTGGCTCAAGCCCAGCGCCTTCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAg                                                                                                                                           <  1:252141/149‑1 (MQ=255)
             cagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAgc                                                                                                                                          <  1:1339731/150‑1 (MQ=255)
              agcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAgcg                                                                                                                                         <  2:2015204/150‑1 (MQ=255)
                      gTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCt                                                                                                                                 >  1:1577100/1‑150 (MQ=255)
                       tCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTg                                                                                                                                >  1:154869/1‑150 (MQ=255)
                         aTCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAc                                                                                                                              >  1:1913986/1‑150 (MQ=255)
                         aTCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAc                                                                                                                              >  1:2344761/1‑150 (MQ=255)
                          tCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCa                                                                                                                                      <  2:1186068/141‑1 (MQ=255)
                          tCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCa                                                                                                                                      >  1:1186068/1‑141 (MQ=255)
                          tCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcg                                                                                                                             <  1:1987560/150‑1 (MQ=255)
                           cGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                            <  2:2347192/150‑1 (MQ=255)
                           cGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                            <  2:154869/150‑1 (MQ=255)
                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                     >  1:671124/1‑150 (MQ=255)
                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                     >  2:2235473/1‑150 (MQ=255)
                                    gatcttAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATCGCCATTTGTCCATCAGCGCCACGGCTGACGCGCAACTCGc                                                                                                                   <  1:1398519/145‑1 (MQ=255)
                                    cGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGTCCATCAGCGCCACGGCTGACGCGCAACTCGc                                                                                                                   <  2:948507/150‑1 (MQ=255)
                                      aCATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCa                                                                                                                 <  2:2150143/150‑1 (MQ=255)
                                        aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACt                                                                                                               >  2:1098949/1‑150 (MQ=255)
                                          aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGGCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGAGCGTCTGCGCCACGGCTGACGCGCAACTCGCCagatc                                                                                                             >  2:1398519/1‑146 (MQ=255)
                                          aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                             <  2:190165/150‑1 (MQ=255)
                                            tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                           <  2:586666/150‑1 (MQ=255)
                                               tcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTa                                                                                                        >  2:473295/1‑150 (MQ=255)
                                                 gCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAtt                                                                                                      <  1:1098949/150‑1 (MQ=255)
                                                       gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCAtg                                                                                                                                             <  2:1494679/105‑2 (MQ=255)
                                                       gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCAtg                                                                                                                                             >  1:1494679/1‑104 (MQ=255)
                                                        tACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcga                                                                                               >  1:797652/1‑150 (MQ=255)
                                                         acacCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaa                                                                                              <  1:889897/150‑1 (MQ=255)
                                                             cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCAc                                                                                          >  2:1115030/1‑150 (MQ=255)
                                                              gCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACt                                                                                         >  2:1764369/1‑150 (MQ=255)
                                                                 ttCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgc                                                                                      <  2:913707/150‑1 (MQ=255)
                                                                  tCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgcg                                                                                     <  1:1061905/150‑1 (MQ=255)
                                                                    gCCTGGCTCAAGCCCAGCGCCTGCAAGGGCGCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGGATTGAAGCGAAGCACTTGCGcc                                                                                   >  1:676513/1‑150 (MQ=255)
                                                                       tGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgg                                                                                                                                            <  2:2345590/90‑3 (MQ=255)
                                                                       tGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgg                                                                                                                                            >  1:2345590/1‑88 (MQ=255)
                                                                        ggCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCt                                                                               <  1:1465609/150‑1 (MQ=255)
                                                                           tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTGGTTTAATCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                            >  1:1330180/1‑150 (MQ=255)
                                                                           tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGATGCGAAGCACTTGCGCCTTCCCAc                                                                            >  2:1946025/1‑150 (MQ=255)
                                                                           tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                            >  1:1820692/1‑150 (MQ=255)
                                                                              aGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACAGCGGGCATTGAAGCGAAGCACTTGCGCCTGCTCACCCa                                                                         >  1:740675/1‑150 (MQ=255)
                                                                                  cAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTg                                                                     >  2:843506/1‑150 (MQ=255)
                                                                                   aGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGc                                                                    <  2:355180/150‑1 (MQ=255)
                                                                                        cTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGaa                                                               <  1:993859/150‑1 (MQ=255)
                                                                                         tGCACAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCCTCAGCGCCACGGCTGACGCGCAACTCGCCCCTGCGGGTATTGGAGCGAAGCACTTGCGCCTGCTCACCCAGTGGCTCGAAc                                                              >  1:1108826/1‑150 (MQ=255)
                                                                                            aaaaGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGc                                                           >  1:291487/1‑150 (MQ=255)
                                                                                                 gCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGTAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTa                                                      >  2:1497365/1‑150 (MQ=255)
                                                                                                  cccGCTGGCATGTCGACCGCTACGGGTTTTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgcgc                                                                                    <  2:1172914/119‑1 (MQ=255)
                                                                                                  cccGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgcgc                                                                                    >  1:1172914/1‑119 (MQ=255)
                                                                                                    cGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGc                                                   <  1:679243/150‑1 (MQ=255)
                                                                                                        ggCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                               <  2:1577100/150‑1 (MQ=255)
                                                                                                         gCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAgg                                              >  1:1039982/1‑150 (MQ=255)
                                                                                                            tGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCa                                           >  2:2191792/1‑150 (MQ=255)
                                                                                                            tGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCa                                           >  2:2163487/1‑150 (MQ=255)
                                                                                                             gTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCaa                                          <  1:473295/150‑1 (MQ=255)
                                                                                                              tCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaa                                                                                              >  1:690404/1‑97 (MQ=255)
                                                                                                              tCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcgaa                                                                                              <  2:690404/97‑1 (MQ=255)
                                                                                                               cGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAgg                                        <  2:1913986/150‑1 (MQ=255)
                                                                                                               cGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAgg                                        <  1:1764369/150‑1 (MQ=255)
                                                                                                                   cGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGGTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGc                                    <  1:1835439/150‑1 (MQ=255)
                                                                                                                     cTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCAt                                  <  2:1039982/150‑1 (MQ=255)
                                                                                                                      tACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATg                                 >  2:2125376/1‑150 (MQ=255)
                                                                                                                       aCGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATgg                                <  1:2125376/150‑1 (MQ=255)
                                                                                                                           gttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGcc                                                                                   <  2:652889/95‑1 (MQ=255)
                                                                                                                           gttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGcc                                                                                   >  1:652889/1‑95 (MQ=255)
                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                     <  2:408162/150‑1 (MQ=255)
                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                     <  1:1115030/150‑1 (MQ=255)
                                                                                                                                   gCGGGGAAGTCAATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGAAACTCGCCACTGCGGGTATTGAAGCGAAGCATTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCg                    <  2:108098/150‑1 (MQ=255)
                                                                                                                                     ggggAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACTTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGAc                  >  1:1967211/1‑150 (MQ=255)
                                                                                                                                     ggggAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGAc                  <  2:1820692/150‑1 (MQ=255)
                                                                                                                                        gAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgc                                                                                      >  1:401860/1‑79 (MQ=255)
                                                                                                                                        gAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgc                                                                                      <  2:401860/79‑1 (MQ=255)
                                                                                                                                        gAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTc               >  1:2393529/1‑150 (MQ=255)
                                                                                                                                           gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTATGCCGCGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc            >  2:572462/5‑150 (MQ=255)
                                                                                                                                           gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc            >  1:261168/5‑150 (MQ=255)
                                                                                                                                           gctgaTGGCCAGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACCGCTTGGATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc            >  2:2143455/5‑150 (MQ=255)
                                                                                                                                             tgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACgg          >  2:1455115/3‑150 (MQ=255)
                                                                                                                                             tgaTGGACAGCAGACCATCACACGCACGACTGACTGGCAACTCTCCACTGCGTGTAGTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTATGCCGAGGCCAAGGTAGCATTTCCGCACAGGTCGACTTCCACgg          >  2:529905/3‑150 (MQ=255)
                                                                                                                                              cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt         >  2:1376058/1‑150 (MQ=255)
                                                                                                                                               aTGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGa                   >  2:2258406/1‑139 (MQ=255)
                                                                                                                                               aTGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGa                   <  1:2258406/139‑1 (MQ=255)
                                                                                                                                               aTGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg        <  2:797652/150‑1 (MQ=255)
                                                                                                                                                tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg       >  2:2007452/1‑150 (MQ=255)
                                                                                                                                                tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg       <  1:606727/150‑1 (MQ=255)
                                                                                                                                                 ggCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTggg      <  1:2235473/150‑1 (MQ=255)
                                                                                                                                                 ggCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTggg      <  2:740675/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg     <  2:1967211/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg     <  1:2007452/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg     <  2:2121798/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCGGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg     <  1:2020450/150‑1 (MQ=255)
                                                                                                                                                     agcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGTg  <  2:671124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGTG  >  NC_002947/4906359‑4906657

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