Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 670,786 A→G 11.6% Q743R (CAG→CGG) ‡ pbpC → peptidoglycan glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947670,7860AG11.6% 217.2 / 17.3 86Q743R (CAG→CGG) ‡pbpCpeptidoglycan glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/39);  new base G (5/5);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCC  >  NC_002947/670638‑670934
                                                                                                                                                    |                                                                                                                                                    
cctgcggcgGGCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCa                                                                                                                                                      <  2:24794/146‑1 (MQ=255)
   gcgcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGcc                                                                                                                                                   >  5:382956/1‑149 (MQ=255)
    cgcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCg                                                                                                                                                  <  6:70882/149‑1 (MQ=255)
     gcCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGt                                                                                                                                                 >  6:185401/1‑149 (MQ=255)
       cGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTg                                                                                                                                               <  2:264062/149‑1 (MQ=255)
        ggCCAGTGCTCCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTTAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGGGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTgg                                                                                                                                              <  7:136363/149‑1 (MQ=255)
        ggCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTga                                                                                                                                              <  5:185401/149‑2 (MQ=255)
        ggCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTga                                                                                                                                              <  2:10695/149‑2 (MQ=255)
         gCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTggg                                                                                                                                             <  7:29359/149‑1 (MQ=255)
          ccAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTgagg                                                                                                                                            <  7:101930/149‑4 (MQ=255)
          ccAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGt                                                                                                                                            <  8:289668/149‑1 (MQ=255)
          ccAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGt                                                                                                                                            <  2:281617/149‑1 (MQ=255)
                   caccGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGAGCCGCCCGGCCACCCGCACTGAGCCGTTGCAATTCCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGCTTCCTCAACCGCCAGCCCTTCCGTGAAACCCAg                                                                                                                                   >  2:224127/1‑149 (MQ=255)
                      cGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAgggg                                                                                                                                >  1:261535/1‑149 (MQ=255)
                         tGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAg                                                                                                                             >  2:139533/1‑149 (MQ=255)
                             gATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGa                                                                                                                           >  4:237476/1‑147 (MQ=255)
                               tcgtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGc                                                                                                                       <  1:392945/149‑1 (MQ=255)
                                 gtcgGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGATGGGTGAAACCCAGGGGCAGGACAGCtt                                                                                                                     <  1:224127/149‑1 (MQ=255)
                                     gCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttgttg                                                                                                                 <  1:372917/149‑1 (MQ=255)
                                            cccAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGtt                                                                                                          >  8:2215/1‑149 (MQ=255)
                                             ccAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGAGttt                                                                                                         >  2:183090/1‑149 (MQ=255)
                                               aGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGc                                                                                                       <  7:313051/149‑1 (MQ=255)
                                                 gCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCca                                                                                                     >  2:269613/1‑149 (MQ=255)
                                                     cAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAgg                                                                                                 >  2:279808/1‑149 (MQ=255)
                                                      aaTCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGt                                                                                                <  2:325117/149‑1 (MQ=255)
                                                      aaTCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGt                                                                                                <  3:237476/149‑1 (MQ=255)
                                                        tCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCg                                                                                              >  4:327432/1‑149 (MQ=255)
                                                          tGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCggg                                                                                            >  3:308040/1‑149 (MQ=255)
                                                            cgccgcCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCggg                                                                                            <  3:146798/147‑1 (MQ=255)
                                                             gccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGCCGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAg                                                                                         >  5:78705/1‑149 (MQ=255)
                                                              ccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTg                                                                                                                                               <  3:320574/94‑1 (MQ=255)
                                                              ccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTg                                                                                                                                               >  4:320574/1‑94 (MQ=255)
                                                                gccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGcc                                                                                      >  3:132488/1‑149 (MQ=255)
                                                                   cGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGaa                                                                                   <  4:308040/149‑1 (MQ=255)
                                                                   cGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGACGTTGGGTGAAACCCAGTGGCAGGACAGCTTGTTGCTGCGTTTGACACAGGTCGGGCAGGCCGaa                                                                                   >  1:311264/1‑149 (MQ=255)
                                                                       caccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCa                                                                               <  7:129074/149‑1 (MQ=255)
                                                                         ccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgc                                                                             <  1:269613/149‑1 (MQ=255)
                                                                          caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttgttg                                                                                                                 >  8:188951/1‑112 (MQ=255)
                                                                          caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttgttg                                                                                                                 <  7:188951/112‑1 (MQ=255)
                                                                          caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcg                                                                            >  8:60734/1‑149 (MQ=255)
                                                                           aGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGt                                                                                                                                            <  7:260025/84‑1 (MQ=255)
                                                                           aGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGt                                                                                                                                            >  8:260025/1‑84 (MQ=255)
                                                                              aCTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCt                                                                        <  6:382956/149‑1 (MQ=255)
                                                                                tGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTg                                                                       >  6:217851/1‑148 (MQ=255)
                                                                                  aGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCa                                                                                                                              <  2:250739/91‑1 (MQ=255)
                                                                                   gCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGAc                                                                    >  3:390502/1‑148 (MQ=255)
                                                                                      gTTGCATTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTTGTTCCTCAACTGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACTGCTTGTTGCTGCGTTTGCCACAGTTCGGTCAGGCCGAACTCAGCGCGCTGGACGaaa                                                                <  1:54934/149‑1 (MQ=255)
                                                                                           aaTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGt                                                           >  2:260855/1‑149 (MQ=255)
                                                                                            aTTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTc                                                          <  7:45575/149‑1 (MQ=255)
                                                                                              tGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGt                                                           >  4:297424/1‑146 (MQ=255)
                                                                                               gCAGGTATCTGCCTTGGGTCGTGGTGCACGGCTCTGGTGCTTACTCAACTGCAAGCCGTTGCGAGAAACACAGGGGAAGGACAGCTAGTTGCTGCGTATGCCACAGGTCGCACAGGCAGAACTCAGCGCGCTGGACGAAAGCGGTCAg                                                        >  4:331439/1‑148 (MQ=255)
                                                                                                 aGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACg                                                     >  6:386818/1‑149 (MQ=255)
                                                                                                   gTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGgc                                                   >  3:84139/1‑149 (MQ=255)
                                                                                                      tCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGgcgc                                                 <  6:201060/148‑1 (MQ=255)
                                                                                                         gCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                             >  1:17201/1‑149 (MQ=255)
                                                                                                          ccTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCCCAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTg                                            >  3:336733/1‑149 (MQ=255)
                                                                                                          ccTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCAGAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                             <  3:331439/148‑1 (MQ=255)
                                                                                                          ccTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCgg                                               <  7:312776/146‑1 (MQ=255)
                                                                                                           cTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTgg                                           >  8:150211/1‑149 (MQ=255)
                                                                                                              gggTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGTAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGt                                        >  4:310226/1‑149 (MQ=255)
                                                                                                                   gtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  3:327432/149‑1 (MQ=255)
                                                                                                                   gtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  7:150211/149‑1 (MQ=255)
                                                                                                                   gtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  3:377938/149‑1 (MQ=255)
                                                                                                                   gtgggtgACGGAGCTGGTGGTTACTAAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  2:311264/143‑1 (MQ=255)
                                                                                                                   ggggtggCCGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                   <  7:60734/147‑1 (MQ=255)
                                                                                                                      gtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTg                                                                       >  2:251188/1‑110 (MQ=255)
                                                                                                                       tggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGTGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCa                               <  1:139533/149‑1 (MQ=255)
                                                                                                                         gACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGt                             <  4:84139/149‑1 (MQ=255)
                                                                                                                            ggCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAg                          >  2:373639/1‑149 (MQ=255)
                                                                                                                             gcgcTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGTAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGt                         <  3:310226/149‑1 (MQ=255)
                                                                                                                                   gtggtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGgt                     <  1:183090/147‑1 (MQ=255)
                                                                                                                                    tggtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGtt                  >  8:106859/1‑149 (MQ=255)
                                                                                                                                      gtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTc                 <  2:261535/148‑1 (MQ=255)
                                                                                                                                       ttCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcga               >  2:68365/1‑149 (MQ=255)
                                                                                                                                          cccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAg                       <  8:16277/136‑1 (MQ=255)
                                                                                                                                          cccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAg                       >  7:16277/3‑138 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGGGTTCGAGGTCAGGGAGTAGGTGTTcgccga            >  4:78077/1‑148 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga            >  5:297701/1‑148 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga            >  3:165428/1‑148 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga            >  6:387498/1‑148 (MQ=255)
                                                                                                                                           ccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga            >  5:74827/2‑148 (MQ=255)
                                                                                                                                           ccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGAGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCgacgac           >  8:202925/2‑149 (MQ=255)
                                                                                                                                             aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACTGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGGTCCAGGTCAGTGAGTAGGTGTTCGACGACa          <  4:44539/148‑1 (MQ=255)
                                                                                                                                             aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAg         >  7:132209/1‑149 (MQ=255)
                                                                                                                                               cGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGcc       <  8:132209/149‑1 (MQ=255)
                                                                                                                                                 gCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCgg     <  3:297424/149‑1 (MQ=255)
                                                                                                                                                   cAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGc   <  7:31846/149‑1 (MQ=255)
                                                                                                                                                    aGCCGTTGGGTGAAAGCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGcc  <  7:106859/149‑1 (MQ=255)
                                                                                                                                                    aGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGCGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGcc  <  1:279808/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCC  >  NC_002947/670638‑670934

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