Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,622 T→A 6.9% I256N (ATT→AAT)  ddlB → D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,526,6220TA 6.9% 267.1 / 7.8 88I256N (ATT→AAT) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base T (41/40);  new base A (3/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 = 9.93e-01

GCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACT  >  NC_002947/1526473‑1526769
                                                                                                                                                      |                                                                                                                                                   
gccgccGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCAt                                                                                                                                                     >  2:1808195/1‑150 (MQ=255)
 ccgccgATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCAtt                                                                                                                                                    >  1:920104/1‑150 (MQ=255)
  cgccgATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTg                                                                                                                                                   <  2:1888864/150‑1 (MQ=255)
  cgccgATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTg                                                                                                                                                   <  1:1224160/150‑1 (MQ=255)
     cgATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCa                                                                                                                                                <  2:1176651/150‑1 (MQ=255)
      gATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCAt                                                                                                                                               >  2:1248775/1‑150 (MQ=255)
       aTCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATc                                                                                                                                              >  1:1015870/1‑150 (MQ=255)
           cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCCTCCCGTGCCGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACTGCGCGC‑GCCTTCGATGCCAATGGCATCGCAg                                                                                                                                          <  2:691467/150‑1 (MQ=255)
           cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTTCGATGCCCATGGCATCGCAg                                                                                                                                          <  2:552226/150‑1 (MQ=255)
           cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAg                                                                                                                                          <  1:369544/150‑1 (MQ=255)
            ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGCTGCCATTGGCATCGAAgg                                                                                                                                         >  2:1407569/1‑150 (MQ=255)
                   aCGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGGCAAGGAACAGGAAGTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCAACGAAGGCTggggg                                                                                                                                  >  2:1032949/1‑150 (MQ=255)
                     gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGAGCGCCCCCAGCCTTCGATGCCAATGGCATCGCAGGcgcgcgcg                                                                                                                                 <  2:1263395/150‑8 (MQ=37)
                     gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                >  1:1430240/1‑150 (MQ=255)
                     gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                >  1:1886341/1‑150 (MQ=255)
                        gcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTTGATAGCGCCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATTCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                             <  2:818471/150‑1 (MQ=255)
                        gcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                             <  2:920104/150‑1 (MQ=255)
                          aCGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTgg                                                                                                                           >  1:735287/1‑150 (MQ=255)
                              gtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGt                                                                                                                       <  1:1808195/150‑1 (MQ=255)
                                  tACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATg                                                                                                                   >  2:87038/1‑150 (MQ=255)
                                    cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGGATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                 <  1:614171/150‑1 (MQ=255)
                                              gCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAg                                                                                                       >  2:642622/1‑150 (MQ=255)
                                               ccAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGCCCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                                      <  1:1180609/150‑1 (MQ=255)
                                               ccAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                                      >  2:1307259/1‑150 (MQ=255)
                                                 aaGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGc                                                                                                    <  1:993194/150‑1 (MQ=255)
                                                 aaGTACATCGTCAACGACCGCGATTACCGGAGCCCGTGCGGGCTGGATAGCGCCAACGACCAGGAACTGATCGACCTGTCCGCCCGC‑GCCTGCGATTCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGc                                                                                                    <  1:1039026/150‑1 (MQ=255)
                                                    tACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCgg                                                                                                 <  2:525995/150‑1 (MQ=255)
                                                       aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                     <  2:1832494/111‑1 (MQ=255)
                                                       aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                     >  1:1832494/1‑111 (MQ=255)
                                                       aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTc                                                                                              >  1:1164441/1‑150 (MQ=255)
                                                        tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                             <  1:87038/150‑1 (MQ=255)
                                                           tCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                 <  2:1831843/127‑1 (MQ=255)
                                                           tCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                 >  1:1831843/1‑127 (MQ=255)
                                                           tCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgc                                                                                          >  2:276340/1‑150 (MQ=255)
                                                            cAACGATACCCAGTACCGAATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgct                                                                                         <  1:1411207/150‑1 (MQ=255)
                                                                 aTACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGaa                                                                                                                                           >  2:516928/1‑95 (MQ=255)
                                                                 aTACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGaa                                                                                                                                           <  1:516928/95‑1 (MQ=255)
                                                                   aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                  >  1:461076/1‑150 (MQ=255)
                                                                   aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                  <  1:1805966/150‑1 (MQ=255)
                                                                     ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                >  1:169793/1‑150 (MQ=255)
                                                                     ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                >  2:149757/1‑150 (MQ=255)
                                                                     ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                >  2:843460/1‑150 (MQ=255)
                                                                     ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGAGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                >  1:99856/1‑150 (MQ=255)
                                                                      cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAATGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCACTGGCATCGAAGGCTGGGGGCGGTTGCACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTc                                                                               <  1:670050/150‑1 (MQ=255)
                                                                           ccGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAacac                                                                          >  1:101191/1‑150 (MQ=255)
                                                                             gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                     <  1:14207/121‑1 (MQ=255)
                                                                             gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                     >  2:14207/1‑121 (MQ=255)
                                                                              cATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACATGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGc                                                                       >  2:1881750/1‑150 (MQ=255)
                                                                                    gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAgg                                                                 <  2:1676849/150‑1 (MQ=255)
                                                                                     tGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGc                                                                <  1:1577217/150‑1 (MQ=255)
                                                                                       cGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                              <  2:1430240/150‑1 (MQ=255)
                                                                                        gggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATg                                                             >  1:929004/1‑150 (MQ=255)
                                                                                         ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                            <  1:1582164/150‑1 (MQ=255)
                                                                                              gATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGACGGCTGGGCGCTGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGAt                                                       <  1:1175635/150‑1 (MQ=255)
                                                                                               aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATc                                                      <  1:149757/150‑1 (MQ=255)
                                                                                                      cAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGCACATAgacc                                               <  2:1461378/150‑4 (MQ=255)
                                                                                                       aaGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTg                                              >  1:269374/1‑150 (MQ=255)
                                                                                                        aggaacaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTgg                                             <  2:1015870/150‑1 (MQ=255)
                                                                                                        aggaacaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTgg                                             >  1:373407/1‑150 (MQ=255)
                                                                                                        aggaacaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTgg                                             >  2:209138/1‑150 (MQ=255)
                                                                                                         ggaacaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGt                                            <  2:1165875/150‑1 (MQ=255)
                                                                                                            acaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGcc                                         >  1:52221/1‑150 (MQ=255)
                                                                                                              aggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                                       <  2:2845/150‑1 (MQ=255)
                                                                                                              aggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                                       <  1:1248775/150‑1 (MQ=255)
                                                                                                               ggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCAt                                      >  1:332654/1‑150 (MQ=255)
                                                                                                                  acTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggc                                   >  2:1522142/1‑150 (MQ=255)
                                                                                                                     gATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggcggc                                <  2:1886341/150‑1 (MQ=255)
                                                                                                                           ccTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGc                          <  1:1522142/150‑1 (MQ=255)
                                                                                                                              gACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGt                                                                                                >  2:769525/1‑77 (MQ=255)
                                                                                                                              gACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGt                                                                                                <  1:769525/77‑1 (MQ=255)
                                                                                                                              gACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGc                       <  2:1613850/150‑1 (MQ=255)
                                                                                                                              gACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGc                       >  2:1872508/1‑150 (MQ=255)
                                                                                                                               aCCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCt                      <  1:549539/150‑1 (MQ=255)
                                                                                                                                 cgcgcgc‑gCCTGCGATGCCCTTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTgg                    >  2:595201/1‑150 (MQ=255)
                                                                                                                                  gcgcgc‑gcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGc                   >  2:542194/1‑150 (MQ=255)
                                                                                                                                     cgc‑gcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                <  1:245496/150‑1 (MQ=255)
                                                                                                                                          ccTTCGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGa                                                                                   >  1:27009/1‑79 (MQ=38)
                                                                                                                                          ccTTCGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACt            >  1:1257419/1‑150 (MQ=255)
                                                                                                                                           ccgCGGTGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACtt           >  2:1919013/3‑150 (MQ=255)
                                                                                                                                           cTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACtt           >  2:304674/1‑150 (MQ=255)
                                                                                                                                            ttcGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTc          >  1:1073039/3‑150 (MQ=255)
                                                                                                                                            tGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTc          >  1:384522/1‑150 (MQ=255)
                                                                                                                                            tGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTc          <  2:52221/150‑1 (MQ=255)
                                                                                                                                              cGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGc                                          <  1:43624/116‑1 (MQ=255)
                                                                                                                                              cGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGc                                          >  2:43624/1‑116 (MQ=255)
                                                                                                                                              cGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCa        <  2:332654/150‑1 (MQ=255)
                                                                                                                                               gATGCCATTGGCATCGAAGGCTGGGGGCGGTTGTACGTGATGCGGGACGAGCGGGGCGGGTTTTGGCTTCTGTGAGTCAACACGGCACCAGGCATGCCAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAg       >  2:850929/1‑150 (MQ=255)
                                                                                                                                                    cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCGCCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACt  >  1:301597/1‑150 (MQ=255)
                                                                                                                                                    cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTcc         <  2:1408112/143‑1 (MQ=255)
                                                                                                                                                    cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTcc         >  1:1408112/1‑143 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACT  >  NC_002947/1526473‑1526769

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