Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 670,785 C→T 11.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,7850CT11.7% 221.8 / 15.3 86Q743* (CAG→TAG) ‡pbpCpeptidoglycan glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base T (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 = 8.38e-01

CAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGC  >  NC_002947/670638‑670933
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                              >  4:320574/1‑94 (MQ=255)
                                                              ccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCTGGCCGTTg                                                                                                                                              <  3:320574/94‑1 (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                                                                                                                <  7:188951/112‑1 (MQ=255)
                                                                          caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCttgttg                                                                                                                >  8:188951/1‑112 (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                                  <  3:377938/149‑1 (MQ=255)
                                                                                                                   gtggtggACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                  <  7:150211/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                      >  7:16277/3‑138 (MQ=255)
                                                                                                                                          cccAACGGCTGGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAg                      <  8:16277/136‑1 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGGGTTCGAGGTCAGGGAGTAGGTGTTcgccga           >  4:78077/1‑148 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga           >  3:165428/1‑148 (MQ=255)
                                                                                                                                           tCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgacga           >  5:297701/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)
                                                                                                                                                   |                                                                                                                                                    
CAGGCGCCGGCCAGTGCACCACCGCTGTCGATCGTCGGCGTGCGCCCAGGCGACAATCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGC  >  NC_002947/670638‑670933

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: