Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 670,788 C→T 5.7% 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.7% 717.6 / 21.0 228Q743* (CAG→TAG) ‡pbpCpeptidoglycan glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (101/114);  new base T (5/8);  total (106/122)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 234 total aligned reads displayed.

CGCAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGC  >  NC_002947/670639‑670936
                                                                                                                                                     |                                                                                                                                                    
cgcAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGcc                                                                                                                                                      >  2:3270804/1‑150 (MQ=255)
cgcAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGcc                                                                                                                                                      >  2:20863/1‑150 (MQ=255)
   aGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGc                                                                                                                                                   <  2:4469214/150‑1 (MQ=255)
   aGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGc                                                                                                                                                   >  1:3939968/1‑150 (MQ=255)
      cgcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGt                                                                                                                                                >  1:5131151/1‑150 (MQ=255)
       gcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGtt                                                                                                                                               <  1:5318686/150‑1 (MQ=255)
       gcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGtt                                                                                                                                               <  1:4162894/150‑1 (MQ=255)
        ccGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTg                                                                                                                                              <  1:154649/150‑1 (MQ=255)
         cGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                             <  2:647691/150‑1 (MQ=255)
           gCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGt                                                                                                                                           <  2:2029252/150‑1 (MQ=255)
           gCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGt                                                                                                                                           <  1:4671447/150‑1 (MQ=255)
             cAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGa                                                                                                                                         <  1:2987724/150‑1 (MQ=255)
              aGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGaa                                                                                                                                        <  1:4892254/150‑1 (MQ=255)
              aGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGaa                                                                                                                                        >  2:3836507/1‑150 (MQ=255)
              aGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGaa                                                                                                                                        >  1:3198029/1‑150 (MQ=255)
               gTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTgaggaacc                                                                                                                                       <  2:1465857/150‑8 (MQ=255)
                   accaccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCa                                                                                                                                   <  1:346752/150‑1 (MQ=255)
                     caccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGACACCCAgg                                                                                                                                 >  1:4047471/1‑150 (MQ=255)
                     caccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAgg                                                                                                                                 >  2:5173629/1‑150 (MQ=255)
                          cTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAg                                                                                                                            <  1:5173629/150‑1 (MQ=255)
                               gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAg                                                                                                                       >  2:3788256/1‑150 (MQ=255)
                                  cgtcgGCGTCCGCCCCGGCGCCCATCTGCGCCGCCCGGCCCCCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCCCCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCtt                                                                                                                    <  2:4337698/150‑1 (MQ=255)
                                   gtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttg                                                                                                                   <  1:1555014/150‑1 (MQ=255)
                                    tcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAg                                                                                                                            >  1:1259625/1‑140 (MQ=255)
                                    tcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttgt                                                                                                                  >  1:1035490/1‑150 (MQ=255)
                                       gCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTtgc                                                                                                               >  2:2556680/1‑150 (MQ=255)
                                             gcCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGtt                                                                                                         >  1:1405839/1‑150 (MQ=255)
                                              cccAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAggg                                                                                                                                >  1:3607595/1‑126 (MQ=255)
                                              cccAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAggg                                                                                                                                <  2:3607630/126‑1 (MQ=255)
                                                cAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGc                                                                                                      <  1:2227899/150‑1 (MQ=255)
                                                 aGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGcc                                                                                                     >  2:3046781/1‑150 (MQ=255)
                                                   gCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCcac                                                                                                   >  2:180021/1‑150 (MQ=255)
                                                     gACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTTCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAg                                                                                                 >  1:3355657/1‑150 (MQ=255)
                                                     gACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAg                                                                                                 >  2:2469599/1‑150 (MQ=255)
                                                          tCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCgg                                                                                            >  1:4260222/1‑150 (MQ=255)
                                                          tCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCgg                                                                                            <  1:3422264/150‑1 (MQ=255)
                                                             gcgcCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTTCTGCGTTTGCCACAGGTCGGGCa                                                                                         <  2:167542/150‑1 (MQ=255)
                                                               gccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAgg                                                                                       >  1:4981598/1‑150 (MQ=255)
                                                               gccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAgg                                                                                       >  2:4250054/1‑150 (MQ=255)
                                                                ccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGc                                                                                      <  1:1176082/150‑1 (MQ=255)
                                                                 cgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                             >  2:2687357/1‑94 (MQ=255)
                                                                 cgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                             <  1:2687357/94‑1 (MQ=255)
                                                                      ggCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACt                                                                                <  1:1927123/150‑1 (MQ=255)
                                                                        ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                              <  1:4804407/150‑1 (MQ=255)
                                                                        ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                              <  2:1182899/150‑1 (MQ=255)
                                                                        ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                              <  2:3183720/150‑1 (MQ=255)
                                                                        ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                              <  2:105097/150‑1 (MQ=255)
                                                                        ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                              <  2:4654983/150‑1 (MQ=255)
                                                                           ccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcg                                                                           >  1:591200/1‑150 (MQ=255)
                                                                            caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCtt                                                                                                                    <  1:748227/108‑1 (MQ=255)
                                                                            caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCtt                                                                                                                    >  2:748227/1‑108 (MQ=255)
                                                                            caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcgc                                                                          >  2:1509199/1‑150 (MQ=255)
                                                                             aGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcgcg                                                                         <  1:3046781/150‑1 (MQ=255)
                                                                                       cGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGaaa                                                               <  2:3925004/150‑1 (MQ=255)
                                                                                       cGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGaaa                                                               <  1:1803697/150‑1 (MQ=255)
                                                                                        gTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAg                                                              >  1:2072613/1‑150 (MQ=255)
                                                                                         ttGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGc                                                             <  1:3788256/150‑1 (MQ=255)
                                                                                           gCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCgg                                                           <  1:4087996/150‑1 (MQ=255)
                                                                                           gCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCgg                                                           <  2:3557191/150‑1 (MQ=255)
                                                                                            cAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGt                                                          <  1:464824/150‑1 (MQ=255)
                                                                                             gccctcAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTc                                                         <  1:2498888/145‑1 (MQ=255)
                                                                                             aaTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTc                                                         >  2:103725/1‑150 (MQ=255)
                                                                                              aTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCa                                                        <  1:801618/150‑1 (MQ=255)
                                                                                               ttGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAg                                                       <  1:2903106/150‑1 (MQ=255)
                                                                                                        tCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCgg                                              <  1:2639073/150‑1 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                            >  2:2175566/1‑150 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                            >  1:3021048/1‑150 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                            >  1:1883323/1‑150 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                            >  2:2943656/1‑150 (MQ=255)
                                                                                                          tGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                            >  2:3726102/1‑150 (MQ=255)
                                                                                                           gCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTg                                           >  1:4495036/1‑150 (MQ=255)
                                                                                                           gCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTg                                           >  2:3534548/1‑150 (MQ=255)
                                                                                                             cTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGc                                         <  2:1921699/150‑1 (MQ=255)
                                                                                                             cTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCGGGACGAAAGCGGTCAGACGGCGCGGGTGGc                                         <  2:5131151/150‑1 (MQ=255)
                                                                                                                  gtggtggtgGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTcc                                    <  2:53895/150‑1 (MQ=255)
                                                                                                                   tggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                   <  2:4823156/150‑1 (MQ=255)
                                                                                                                   tggtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                   <  2:4981598/150‑1 (MQ=255)
                                                                                                                     gtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAgg                                 >  2:2414735/1‑150 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                  >  2:4792166/1‑146 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                              >  2:872899/1‑150 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                              >  1:3737827/1‑150 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                              <  2:1264502/150‑1 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                              >  2:2746010/1‑150 (MQ=255)
                                                                                                                        gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                              >  1:667461/1‑150 (MQ=255)
                                                                                                                           gACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTg                           >  2:5347693/1‑150 (MQ=255)
                                                                                                                               gcgcTTGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAgga                       <  1:473851/150‑3 (MQ=255)
                                                                                                                               gcgcTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTa                       <  1:3201173/150‑1 (MQ=255)
                                                                                                                                     gtggtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGtt                 >  1:4300776/1‑150 (MQ=255)
                                                                                                                                       tttcacccAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcg               >  1:742921/8‑150 (MQ=255)
                                                                                                                                        ttcacccAACGGCTGGCCGTTGAGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcga              >  2:5091556/7‑150 (MQ=255)
                                                                                                                                        gtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcga              >  1:4427016/1‑150 (MQ=255)
                                                                                                                                         ttCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgac             <  2:3198048/150‑1 (MQ=255)
                                                                                                                                             tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTATGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACa         >  2:2709377/1‑150 (MQ=255)
                                                                                                                                             ccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACa         >  2:4863621/2‑150 (MQ=255)
                                                                                                                                              cAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAg        >  2:5215132/1‑150 (MQ=255)
                                                                                                                                               aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGt                        <  1:4525477/133‑1 (MQ=255)
                                                                                                                                               aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGt                        >  2:4525477/1‑133 (MQ=255)
                                                                                                                                               aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGc       >  1:2791504/1‑150 (MQ=255)
                                                                                                                                                    ccAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGc  >  2:619061/1‑150 (MQ=255)
                                                                                                                                                     cAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTc                                                         >  2:3010078/1‑94 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGC  >  NC_002947/670639‑670936

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