Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 670,788 C→T 5.9% Q743* (CAG→TAG) ‡ pbpC → peptidoglycan glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947670,7880CT 5.9% 449.8 / 11.9 151Q743* (CAG→TAG) ‡pbpCpeptidoglycan glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (58/84);  new base T (6/3);  total (64/87)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 157 total aligned reads displayed.

CGCAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCC  >  NC_002947/670639‑670937
                                                                                                                                                     |                                                                                                                                                     
cgcAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGcc                                                                                                                                                       >  1:2182743/1‑150 (MQ=255)
       gcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGtt                                                                                                                                                >  2:1526478/1‑150 (MQ=255)
         cGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                              <  1:2195323/150‑1 (MQ=255)
          ggCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCGATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCCTACCGGCTGGCCGtttag                                                                                                                                             <  2:2112053/150‑4 (MQ=255)
          ggCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTggg                                                                                                                                             >  2:438020/1‑150 (MQ=255)
               gTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGaaa                                                                                                                                        >  1:493220/1‑150 (MQ=255)
                   accaccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCa                                                                                                                                    <  1:380813/150‑1 (MQ=255)
                   accaccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCa                                                                                                                                    >  1:1545002/1‑150 (MQ=255)
                     caccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAgg                                                                                                                                  >  2:2217126/1‑150 (MQ=255)
                      accGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAggg                                                                                                                                 <  1:2243239/150‑1 (MQ=255)
                        cGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCATTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAggg                                                                                                                                 >  1:1811364/1‑148 (MQ=255)
                          cTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAg                                                                                                                             <  2:2321161/150‑1 (MQ=255)
                           tGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAgg                                                                                                                            <  1:952524/150‑1 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTga                                                                                                                                              <  2:1967097/128‑2 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTga                                                                                                                                              >  1:1967097/1‑127 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAg                                                                                                                        >  2:758602/1‑150 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAg                                                                                                                        >  1:1630890/1‑150 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAg                                                                                                                        >  1:1301760/1‑150 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAg                                                                                                                        >  2:126259/1‑150 (MQ=255)
                                aTCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGc                                                                                                                       <  2:1071541/150‑1 (MQ=255)
                                 tcgtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCt                                                                                                                      >  1:2078543/1‑150 (MQ=255)
                                   gtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttg                                                                                                                    <  1:1347485/150‑1 (MQ=255)
                                   gtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttg                                                                                                                    <  2:769972/150‑1 (MQ=255)
                                   gtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttg                                                                                                                    <  2:225411/150‑1 (MQ=255)
                                     cgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttgtt                                                                                                                  <  1:1353808/150‑1 (MQ=255)
                                       gCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTtgc                                                                                                                >  2:1534433/1‑150 (MQ=255)
                                        cGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTtgct                                                                                                               >  1:888221/1‑150 (MQ=255)
                                             gcCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGtt                                                                                                          <  1:2009573/150‑1 (MQ=255)
                                               ccAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTg                                                                                                        >  2:940092/1‑150 (MQ=255)
                                               ccAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTg                                                                                                        <  2:1237948/150‑1 (MQ=255)
                                                cAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGc                                                                                                       <  1:1857475/150‑1 (MQ=255)
                                                  ggCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCca                                                                                                     >  2:498605/1‑150 (MQ=255)
                                                   gCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCca                                                                                                     >  1:1891796/1‑149 (MQ=255)
                                                   gCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCca                                                                                                     <  2:1891796/149‑1 (MQ=255)
                                                   gCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTGGTTGCTGCGTTTGCcac                                                                                                    >  2:2426584/1‑150 (MQ=255)
                                                    cGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCcaca                                                                                                   <  2:2078543/150‑1 (MQ=255)
                                                     gACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAg                                                                                                  >  2:1724632/1‑150 (MQ=255)
                                                          tCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                              >  2:1354383/1‑101 (MQ=255)
                                                           cTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGTGCTGGTCGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGTCAGGACAGCTTGTTGCTGCGTTTGCTACAGGTCggg                                                                                            <  1:789551/150‑1 (MQ=255)
                                                           cTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                              >  2:1523195/1‑100 (MQ=255)
                                                           cTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                              <  1:1523195/100‑1 (MQ=255)
                                                              cgccgcCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAg                                                                                         <  1:575571/150‑1 (MQ=255)
                                                                 cgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGcc                                                                                      >  1:329622/1‑150 (MQ=255)
                                                                   cccGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGa                                                                                    >  1:1060009/1‑150 (MQ=255)
                                                                        ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                               <  1:1526478/150‑1 (MQ=255)
                                                                            caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcgc                                                                           >  2:1055496/1‑150 (MQ=255)
                                                                              gCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcgcgc                                                                         <  2:2182755/150‑1 (MQ=255)
                                                                                      ccGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTCGCCACAGGTCGGGCAGGCCGCACTGAGCGCGATggcagca                                                                 >  2:2231503/1‑145 (MQ=255)
                                                                                       cGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGaaa                                                                <  1:2217126/150‑1 (MQ=255)
                                                                                         ttGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTTTTGCTGCGtt                                                                                                          >  2:1784119/1‑106 (MQ=255)
                                                                                         ttGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTTTTGCTGCGtt                                                                                                          <  1:1784119/106‑1 (MQ=255)
                                                                                          tGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGc                                                                                           >  1:706881/1‑120 (MQ=255)
                                                                                          ggCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCg                                                             <  2:466241/149‑1 (MQ=255)
                                                                                           gCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGGCAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCgg                                                            >  1:448225/1‑150 (MQ=255)
                                                                                           gCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCgg                                                            >  1:1434390/1‑150 (MQ=255)
                                                                                             aaTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGTGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTc                                                          <  1:758602/150‑1 (MQ=255)
                                                                                              aTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCa                                                         <  1:2231503/150‑1 (MQ=255)
                                                                                                tGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGACGCGATGAGCGCGCTGGACGAAAGCGGTCAGa                                                       <  1:1242627/150‑1 (MQ=255)
                                                                                                tGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGATAGTAGTAAAGTGATGGTCAGCGGGTGGGTGAAACACCGGACTGGGACAGCTTGTTGCGGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGa                                                       <  1:534356/150‑1 (MQ=255)
                                                                                                 gCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGAc                                                      >  1:1734021/1‑150 (MQ=255)
                                                                                                  cAGGTATCTGCCTTGGGTGGTGGTGGACGGTGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACg                                                     <  1:803980/150‑1 (MQ=255)
                                                                                                   aGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCtt                                                                                                                     >  1:918784/1‑85 (MQ=255)
                                                                                                     gTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGgcg                                                  >  1:1702578/1‑150 (MQ=255)
                                                                                                       aTCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGcgcg                                                >  1:376334/1‑150 (MQ=255)
                                                                                                        tCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGGCAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCgg                                               <  2:448225/150‑1 (MQ=255)
                                                                                                        tCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCgg                                               <  1:1238668/150‑1 (MQ=255)
                                                                                                         cTGCCTTGGGGGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGACAGCGGTCAGACGGCGCggg                                              <  1:1120260/150‑1 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCgg                                               >  2:677819/1‑148 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCgg                                               <  1:677819/148‑1 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                             >  1:442805/1‑150 (MQ=255)
                                                                                                           gCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTg                                            >  1:2037509/1‑150 (MQ=255)
                                                                                                             cTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGc                                          >  2:1536259/1‑150 (MQ=255)
                                                                                                              ttGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCg                                         >  2:813202/1‑150 (MQ=255)
                                                                                                              ttGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCg                                         >  1:136885/1‑150 (MQ=255)
                                                                                                                   tggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCa                                                                                                                                    <  2:864904/54‑1 (MQ=255)
                                                                                                                   tggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                    <  1:126259/150‑1 (MQ=255)
                                                                                                                   tggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                    <  1:2353934/150‑1 (MQ=255)
                                                                                                                   tggtggtggACGGCGCTGGTGGTTCCTCAACGACCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                    <  1:1320780/150‑1 (MQ=255)
                                                                                                                    ggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  1:1534433/150‑1 (MQ=255)
                                                                                                                    ggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  1:1765923/150‑1 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                               <  2:696841/150‑1 (MQ=255)
                                                                                                                            aCGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGa                           <  2:376334/150‑1 (MQ=255)
                                                                                                                              ggCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGt                         <  2:451689/150‑1 (MQ=255)
                                                                                                                                cgcTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAg                       <  1:2384555/150‑1 (MQ=255)
                                                                                                                                cgcTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTACTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAg                       <  1:273578/150‑1 (MQ=255)
                                                                                                                                     gtggtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGtt                  <  2:1545002/150‑1 (MQ=255)
                                                                                                                                     gtggtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGCGTGGGCCGAAACCAGCGACAAGCTCGTGAATGATCAGAGTGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGtt                  >  1:534719/1‑150 (MQ=255)
                                                                                                                                       tttcacccAACGGCTGGCCGTTGAGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcg                >  2:412990/8‑150 (MQ=255)
                                                                                                                                        gtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcga               >  2:1581264/1‑150 (MQ=255)
                                                                                                                                           acccAACGGCTGGCCGTTGGGTGAAACCCAGGGGAAGGACAGCGTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga            >  1:2299429/4‑150 (MQ=255)
                                                                                                                                           acccAACGGCTGGCCGTTGAGTGAAACCCAGGGGCAGGACAGCTCGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga            >  2:245438/4‑150 (MQ=255)
                                                                                                                                            cccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGcgcg                                                <  1:2240701/111‑1 (MQ=255)
                                                                                                                                            cTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                    <  2:135301/125‑1 (MQ=255)
                                                                                                                                             tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACa          >  2:275944/1‑150 (MQ=255)
                                                                                                                                             ccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACa          >  1:631923/2‑150 (MQ=255)
                                                                                                                                               aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGCCTATTTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGc        >  1:564646/1‑150 (MQ=255)
                                                                                                                                                 cGGCCAGCCGTTGGGTGAAACCCAGGGGCAAGACAGCTTGTTGCTGCGTTTGACTCTGCGCTGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCg      <  1:625964/150‑1 (MQ=255)
                                                                                                                                                  ggCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCgg     <  2:1630890/150‑1 (MQ=255)
                                                                                                                                                    ccAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGc   <  2:1434390/150‑1 (MQ=255)
                                                                                                                                                     cAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTTCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACg                                                     <  1:1659748/99‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCC  >  NC_002947/670639‑670937

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: