Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,631 A→C 6.7% E259A (GAA→GCA)  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,6310AC 6.7% 252.7 / 6.5 89E259A (GAA→GCA) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base C (3/3);  total (46/43)
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.74e-01

CCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGG  >  NC_002947/1526484‑1526777
                                                                                                                                                    |                                                                                                                                                  
cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCCTCCCGTGCCGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACTGCGCGCGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                                  <  2:691467/150‑1 (MQ=255)
cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTTCGATGCCCATGGCATCGCAg                                                                                                                                                  <  2:552226/150‑1 (MQ=255)
cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAg                                                                                                                                                  <  1:369544/150‑1 (MQ=255)
 ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGCTGCCATTGGCATCGAAgg                                                                                                                                                 >  2:1407569/1‑150 (MQ=255)
        aCGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGGCAAGGAACAGGAAGTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCAACGAAGGCTggggg                                                                                                                                          >  2:1032949/1‑150 (MQ=255)
          gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGAGCGCCCCCAGCCTTCGATGCCAATGGCATCGCAGGcgcgcgcg                                                                                                                                         <  2:1263395/150‑8 (MQ=37)
          gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                        >  1:1886341/1‑150 (MQ=255)
          gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                        >  1:1430240/1‑150 (MQ=255)
             gcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTTGATAGCGCCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATTCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                                     <  2:818471/150‑1 (MQ=255)
             gcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                                     <  2:920104/150‑1 (MQ=255)
               aCGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTgg                                                                                                                                   >  1:735287/1‑150 (MQ=255)
                   gtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGt                                                                                                                               <  1:1808195/150‑1 (MQ=255)
                       tACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATg                                                                                                                           >  2:87038/1‑150 (MQ=255)
                         cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGGATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                         <  1:614171/150‑1 (MQ=255)
                                   gCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAg                                                                                                               >  2:642622/1‑150 (MQ=255)
                                    ccAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGCCCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                                              <  1:1180609/150‑1 (MQ=255)
                                    ccAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                                              >  2:1307259/1‑150 (MQ=255)
                                      aaGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGc                                                                                                            <  1:993194/150‑1 (MQ=255)
                                      aaGTACATCGTCAACGACCGCGATTACCGGAGCCCGTGCGGGCTGGATAGCGCCAACGACCAGGAACTGATCGACCTGTCCGCCCGCGCCTGCGATTCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGc                                                                                                            <  1:1039026/150‑1 (MQ=255)
                                         tACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCgg                                                                                                         <  2:525995/150‑1 (MQ=255)
                                            aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                             >  1:1832494/1‑111 (MQ=255)
                                            aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                             <  2:1832494/111‑1 (MQ=255)
                                            aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTc                                                                                                      >  1:1164441/1‑150 (MQ=255)
                                             tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                                     <  1:87038/150‑1 (MQ=255)
                                                tCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                         >  1:1831843/1‑127 (MQ=255)
                                                tCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                         <  2:1831843/127‑1 (MQ=255)
                                                tCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgc                                                                                                  >  2:276340/1‑150 (MQ=255)
                                                 cAACGATACCCAGTACCGAATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgct                                                                                                 <  1:1411207/150‑1 (MQ=255)
                                                      aTACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGaa                                                                                                                                                   >  2:516928/1‑95 (MQ=255)
                                                      aTACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGaa                                                                                                                                                   <  1:516928/95‑1 (MQ=255)
                                                        aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                          <  1:1805966/150‑1 (MQ=255)
                                                        aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                          >  1:461076/1‑150 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                        >  1:169793/1‑150 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                        >  2:149757/1‑150 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                        >  2:843460/1‑150 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGAGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                        >  1:99856/1‑150 (MQ=255)
                                                           cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAATGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCACTGGCATCGAAGGCTGGGGGCGGTTGCACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTc                                                                                       <  1:670050/150‑1 (MQ=255)
                                                                ccGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAacac                                                                                  >  1:101191/1‑150 (MQ=255)
                                                                  gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                             <  1:14207/121‑1 (MQ=255)
                                                                  gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                             >  2:14207/1‑121 (MQ=255)
                                                                   cATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACATGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGc                                                                               >  2:1881750/1‑150 (MQ=255)
                                                                         gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAgg                                                                         <  2:1676849/150‑1 (MQ=255)
                                                                          tGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGc                                                                        <  1:1577217/150‑1 (MQ=255)
                                                                            cGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                                      <  2:1430240/150‑1 (MQ=255)
                                                                             gggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATg                                                                     >  1:929004/1‑150 (MQ=255)
                                                                              ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                                    <  1:1582164/150‑1 (MQ=255)
                                                                                   gATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGACGGCTGGGCGCTGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGAt                                                               <  1:1175635/150‑1 (MQ=255)
                                                                                    aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATc                                                              <  1:149757/150‑1 (MQ=255)
                                                                                           cAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGCACATAgacc                                                       <  2:1461378/150‑4 (MQ=255)
                                                                                            aaGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTg                                                      >  1:269374/1‑150 (MQ=255)
                                                                                             aggaacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTgg                                                     >  1:373407/1‑150 (MQ=255)
                                                                                             aggaacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTgg                                                     >  2:209138/1‑150 (MQ=255)
                                                                                             aggaacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTgg                                                     <  2:1015870/150‑1 (MQ=255)
                                                                                              ggaacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGt                                                    <  2:1165875/150‑1 (MQ=255)
                                                                                                 acaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGcc                                                 >  1:52221/1‑150 (MQ=255)
                                                                                                   aggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                                               <  1:1248775/150‑1 (MQ=255)
                                                                                                   aggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                                               <  2:2845/150‑1 (MQ=255)
                                                                                                    ggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCAt                                              >  1:332654/1‑150 (MQ=255)
                                                                                                       acTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggc                                           >  2:1522142/1‑150 (MQ=255)
                                                                                                          gATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggcggc                                        <  2:1886341/150‑1 (MQ=255)
                                                                                                                ccTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGc                                  <  1:1522142/150‑1 (MQ=255)
                                                                                                                   gACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGt                                                                                                        <  1:769525/77‑1 (MQ=255)
                                                                                                                   gACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGt                                                                                                        >  2:769525/1‑77 (MQ=255)
                                                                                                                   gACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGc                               <  2:1613850/150‑1 (MQ=255)
                                                                                                                   gACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGc                               >  2:1872508/1‑150 (MQ=255)
                                                                                                                    aCCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCt                              <  1:549539/150‑1 (MQ=255)
                                                                                                                      cgcgcgcgCCTGCGATGCCCTTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTgg                            >  2:595201/1‑150 (MQ=255)
                                                                                                                       gcgcgcgcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGc                           >  2:542194/1‑150 (MQ=255)
                                                                                                                          cgcgcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                        <  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                 >  1:1408112/1‑143 (MQ=255)
                                                                                                                                         cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTcc                 <  2:1408112/143‑1 (MQ=255)
                                                                                                                                                aTCGAAGGCTGGGGGCGGTTGGCCGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg   >  2:1549424/1‑150 (MQ=255)
                                                                                                                                                aTCGAAGGCTGGGGGCGGTTGGCCGTGATCGACACTCTGTGGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg   <  2:152346/150‑1 (MQ=255)
                                                                                                                                                aTCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg   >  1:598056/1‑150 (MQ=255)
                                                                                                                                                aTCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg   >  2:746322/1‑150 (MQ=255)
                                                                                                                                                aTCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg   <  1:755711/150‑1 (MQ=255)
                                                                                                                                                 tCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTgg  >  2:333286/1‑150 (MQ=255)
                                                                                                                                                   gAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGa                                                                >  2:1746519/1‑86 (MQ=255)
                                                                                                                                                   gAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGa                                                                <  1:1746519/86‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGG  >  NC_002947/1526484‑1526777

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: