Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,428,030 C→T 7.4% P660L (CCG→CTG)  PP_3042 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,428,0300CT 7.4% 187.9 / 2.9 66P660L (CCG→CTG) PP_3042terminase large subunit
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCTGTTTGATGAGCCGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGA  >  NC_002947/3427885‑3428178
                                                                                                                                                 |                                                                                                                                                    
aGTCTGTTTGATGAGCCGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCgg                                                                                                                                                    >  3:235770/1‑148 (MQ=255)
aGTCTGTTTGATGAGCCGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGc                                                                                                                                                   <  6:304626/149‑1 (MQ=255)
 gTCTGTTTGATGAGCCGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCa                                                                                                                                                  >  6:216053/1‑149 (MQ=255)
         gatgaGCCGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCAc                                                                                                                                           <  2:96153/148‑1 (MQ=255)
               ccGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGc                                                                                                                                     >  1:236946/1‑148 (MQ=255)
                    tcCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAcgcc                                                                                                                               >  1:190874/1‑149 (MQ=255)
                      cTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAcgccgc                                                                                                                             >  7:145448/1‑149 (MQ=255)
                      cTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAcgccg                                                                                                                              >  5:22918/1‑148 (MQ=255)
                      cTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAcgccg                                                                                                                              >  4:89705/1‑148 (MQ=255)
                       tGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAcgccg                                                                                                                              >  3:5084/1‑147 (MQ=255)
                       tGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAcgccg                                                                                                                              >  2:316674/1‑147 (MQ=255)
                          gAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCt                                                                                                                         >  6:244679/1‑149 (MQ=255)
                           aaGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTc                                                                                                                        >  4:50678/1‑149 (MQ=255)
                            aGCAACCGGTCTTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTcc                                                                                                                       <  6:41694/149‑1 (MQ=255)
                            aGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTc                                                                                                                        >  3:121956/1‑148 (MQ=255)
                                  cGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCagcagc                                                                                                                 <  1:316674/149‑1 (MQ=255)
                                  cGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCagcagc                                                                                                                 <  2:190874/149‑1 (MQ=255)
                                     tCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCgg                                                                                                               >  1:309037/1‑148 (MQ=255)
                                       gTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGcgt                                                                                                            <  4:243682/149‑3 (MQ=255)
                                       gTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTa                                                                                                            >  5:59846/1‑149 (MQ=255)
                                        ttGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCgg                                                                                                               <  6:49176/145‑1 (MQ=255)
                                        ttGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCgg                                                                                                               >  5:49176/1‑145 (MQ=255)
                                          gAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCTGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCt                                                                                                         >  3:298039/1‑149 (MQ=255)
                                          gAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCTGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCt                                                                                                         >  3:303149/1‑149 (MQ=255)
                                           aGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCTGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCt                                                                                                         <  4:298039/148‑1 (MQ=255)
                                           aGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCTGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCt                                                                                                         <  4:303149/148‑1 (MQ=255)
                                             cACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCTGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                      <  1:267143/149‑1 (MQ=255)
                                             cACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                      <  3:89705/149‑1 (MQ=255)
                                               ccAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAaga                                                                                                    <  5:216053/149‑1 (MQ=255)
                                               ccAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAaga                                                                                                    <  6:22918/149‑1 (MQ=255)
                                                 aaCACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAc                                                                                                                                  >  7:125762/1‑117 (MQ=255)
                                                 aaCACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAAc                                                                                                                                  <  8:125762/117‑1 (MQ=255)
                                                      gATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                              >  3:115323/1‑148 (MQ=255)
                                                       aTGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGa                                                                                            <  5:244679/149‑1 (MQ=255)
                                                          acgacgATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCTGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACAt                                                                                         <  8:128691/149‑1 (MQ=255)
                                                            gacgaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATgg                                                                                       >  7:55695/1‑149 (MQ=255)
                                                            gacgaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATgg                                                                                       <  2:236946/149‑1 (MQ=255)
                                                               gaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCCCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                     >  7:36660/1‑148 (MQ=255)
                                                                    cccGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATAcc                                                                               <  4:121956/149‑1 (MQ=255)
                                                                       ggACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAg                                                                            <  8:145448/149‑1 (MQ=255)
                                                                        gacgacGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAgg                                                                           <  8:55695/149‑1 (MQ=255)
                                                                        gacgacGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACAGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAgg                                                                           <  3:50678/149‑1 (MQ=255)
                                                                          cgacgaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGca                                                                         >  8:44327/1‑149 (MQ=255)
                                                                          cgacgaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGca                                                                         <  6:258226/149‑1 (MQ=255)
                                                                           gacgaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGca                                                                         >  8:253036/1‑148 (MQ=255)
                                                                                   aGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTcgcg                                                                <  6:259396/149‑1 (MQ=255)
                                                                                          ccGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCGGAag                                                                                                     <  6:67074/105‑1 (MQ=255)
                                                                                          ccGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCGGAag                                                                                                     >  5:67074/1‑105 (MQ=255)
                                                                                           cGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGa                                                        <  8:315653/149‑1 (MQ=255)
                                                                                             ccAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGc                                                      <  6:59846/149‑1 (MQ=255)
                                                                                                    caccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGc                                                      >  4:98529/1‑142 (MQ=255)
                                                                                                    caccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGc                                                      <  3:98529/142‑1 (MQ=255)
                                                                                                     accGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGAtt                                              <  2:309037/149‑1 (MQ=255)
                                                                                                      ccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACCCCTCGCGGCCGTCGAGCGTGCggtcg                                             >  5:54864/1‑145 (MQ=255)
                                                                                                       cGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTgc                                            <  7:253036/149‑1 (MQ=255)
                                                                                                                   ctccacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACg                                <  8:36660/149‑1 (MQ=255)
                                                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGTCATGGCATATACCCAGGCACACCTCGCgg                                                               <  1:211247/115‑1 (MQ=255)
                                                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                     >  4:227045/1‑93 (MQ=255)
                                                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                     <  3:227045/93‑1 (MQ=255)
                                                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCgg                                                               >  2:211247/1‑115 (MQ=255)
                                                                                                                       acgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgtt                            >  3:347365/1‑149 (MQ=255)
                                                                                                                             cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCg                       <  4:5084/148‑1 (MQ=255)
                                                                                                                              actccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                      >  8:3540/1‑148 (MQ=255)
                                                                                                                                  cacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCCCGTGGCGAACGGGTTGTTCGATACAg                 >  1:65155/1‑149 (MQ=255)
                                                                                                                                  cacCGCCGCCTCGACCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGAAGCGGCTATCTGATGAGACGCTGACATGGCATTTACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAg                 >  2:242741/1‑149 (MQ=255)
                                                                                                                                     cgccgcCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGa              >  3:62772/1‑149 (MQ=255)
                                                                                                                                                 cGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGa  <  4:115323/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGTCTGTTTGATGAGCCGTCTCCTGCGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGA  >  NC_002947/3427885‑3428178

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: ‑