Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,882,143 C→G 6.3% T80S (ACC→AGC) ‡ PP_3425 → RND family transporter MFP subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,882,1430CG 6.3% 860.9 / 35.3 255T80S (ACC→AGC) ‡PP_3425RND family transporter MFP subunit
Reads supporting (aligned to +/- strand):  ref base C (114/125);  new base G (9/7);  total (123/132)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 260 total aligned reads displayed.

CAGGCTCCCGCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCC  >  NC_002947/3881995‑3882292
                                                                                                                                                    |                                                                                                                                                     
       ccGCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGgcg                                                                                                                                               <  1:2633257/150‑1 (MQ=255)
        cGCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGgcgc                                                                                                                                              <  1:1693730/150‑1 (MQ=255)
           tGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCa                                                                                                                                           <  2:4424793/150‑1 (MQ=255)
           tGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCa                                                                                                                                           <  1:1242884/150‑1 (MQ=255)
           tGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCa                                                                                                                                           >  1:1999907/1‑150 (MQ=255)
            gcgcCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAg                                                                                                                                          <  2:2925755/150‑1 (MQ=255)
            gcgcCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGGCGCCTTCGGTGAAGGCGacctg                                                                                                                                          <  2:5605036/150‑3 (MQ=255)
                cGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCa                                                                                                                                      <  2:4244491/150‑1 (MQ=255)
                 gAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCaa                                                                                                                                     <  1:1222643/150‑1 (MQ=255)
                 gAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCaa                                                                                                                                     <  2:1395055/150‑1 (MQ=255)
                 gAAGGTAAGTGGGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCaa                                                                                                                                     <  2:4321856/150‑1 (MQ=255)
                      tAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAg                                                                                                                                <  2:1677728/150‑1 (MQ=255)
                           gtgGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGAc                                                                                                                           >  2:1789456/1‑150 (MQ=255)
                              gCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACctg                                                                                                                        >  2:4091270/1‑150 (MQ=255)
                                 aaGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACctgctg                                                                                                                     >  1:3783485/1‑150 (MQ=255)
                                   ggTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGtg                                                                                                                   >  1:4542266/1‑149 (MQ=255)
                                     tGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTcc                                                                                                                 <  2:2897609/150‑1 (MQ=255)
                                     tGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTcc                                                                                                                 <  1:1042487/150‑1 (MQ=255)
                                      gATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCa                                                                                                                <  2:2611736/150‑1 (MQ=255)
                                       aTCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGTCGCCCAGGTCAAGAAATGCGACCTGCTGTTCAAg                                                                                                               >  1:4199703/1‑150 (MQ=255)
                                       aTCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAg                                                                                                               >  2:3803188/1‑150 (MQ=255)
                                       aTCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGACCAGGTCAAGAAAGGCGACCTGCTGTTCCAg                                                                                                               >  1:1204441/1‑150 (MQ=255)
                                        tCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGa                                                                                                              <  1:3215875/150‑1 (MQ=255)
                                        tCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGa                                                                                                              <  2:1011897/150‑1 (MQ=255)
                                          gAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATc                                                                                                            >  1:3385412/1‑150 (MQ=255)
                                             caacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGAc                                                                                                         <  1:2724770/150‑1 (MQ=255)
                                              aacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCa                                                                                                                                           >  1:4657303/1‑115 (MQ=255)
                                              aacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGAcc                                                                                                        >  1:3333095/1‑150 (MQ=255)
                                              aacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGAcc                                                                                                        >  1:746322/1‑150 (MQ=255)
                                                 cGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCcgc                                                                                                     >  1:3544721/1‑150 (MQ=255)
                                                    tCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGcc                                                                                                  <  1:2445661/150‑1 (MQ=255)
                                                         gAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTc                                                                                             <  2:4036213/150‑1 (MQ=255)
                                                            tGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAg                                                                                          <  1:2051048/150‑1 (MQ=255)
                                                             gggACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGCCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAgg                                                                                         <  2:2828412/150‑1 (MQ=255)
                                                              ggACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGACGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGc                                                                                        >  2:522777/1‑150 (MQ=255)
                                                                aCGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTg                                                                                      <  2:3261358/150‑1 (MQ=255)
                                                                  gAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGaa                                                                                    <  2:3208771/150‑1 (MQ=255)
                                                                  gAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGaa                                                                                    <  2:4660126/150‑1 (MQ=255)
                                                                   aGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAg                                                                                   <  1:4916688/150‑1 (MQ=255)
                                                                     ttCACCGGCCGGCTCGTAGGCCCGGTAACCGTCGTTGTGCGCCCGGGGGTCTCCGGCCAGATCGACCTGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTc                                                                                 <  2:640453/150‑1 (MQ=255)
                                                                     ttCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTc                                                                                 >  1:4333860/1‑150 (MQ=255)
                                                                     ttCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTc                                                                                 >  2:819994/1‑150 (MQ=255)
                                                                      tCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTcc                                                                                <  1:2347665/150‑1 (MQ=255)
                                                                      tCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTcc                                                                                <  1:1602799/150‑1 (MQ=255)
                                                                      tCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTcc                                                                                >  1:1708830/1‑150 (MQ=255)
                                                                        aCCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTccgc                                                                              >  1:5107117/1‑150 (MQ=255)
                                                                             ccgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCt                                                                         <  2:3485089/150‑1 (MQ=255)
                                                                             ccgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCt                                                                         >  2:1987861/1‑150 (MQ=255)
                                                                              cgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTc                                                                        >  1:1900192/1‑150 (MQ=255)
                                                                               gccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCg                                                                       >  1:5478873/1‑150 (MQ=255)
                                                                                  tCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCTGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAg                                                                    >  2:2802589/1‑150 (MQ=255)
                                                                                  tCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAg                                                                    <  1:2311550/150‑1 (MQ=255)
                                                                                  tCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAg                                                                    <  1:980629/150‑1 (MQ=255)
                                                                                   cGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCGAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGc                                                                   <  2:868693/150‑1 (MQ=255)
                                                                                    gAAGCCAAGGAAACCGTCGAAGTGCGCCCGCGGGTGTCCGGCCAGATCGACCAGGTCGGCTTCACCGAAGGCGCGCAGGTCAAGAAAGTCGACCTGCTGTTCAAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGcc                                                                  >  2:4719191/1‑150 (MQ=255)
                                                                                        ccccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGTCCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGTCGTCCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCaag                                                              >  2:5072393/1‑149 (MQ=255)
                                                                                             gAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCTTCGGTGAAGGCGACCTGGt                                                                                                                                        <  2:4995369/71‑1 (MQ=18)
                                                                                             gAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCTTCGGTGAAGGCGACCTGGt                                                                                                                                        >  1:4995369/1‑71 (MQ=18)
                                                                                               aaCCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACa                                                       <  1:5575164/150‑1 (MQ=255)
                                                                                                aCCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACag                                                      >  2:5592126/1‑150 (MQ=255)
                                                                                                aCCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACag                                                      >  2:3556502/1‑150 (MQ=255)
                                                                                                aCCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACag                                                      <  1:5338199/150‑1 (MQ=255)
                                                                                                 ccGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACagg                                                     <  2:5541308/150‑1 (MQ=255)
                                                                                                 ccGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACagg                                                     <  1:1402236/150‑1 (MQ=255)
                                                                                                   gTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggcc                                                   >  2:740755/1‑150 (MQ=255)
                                                                                                    tCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggcca                                                  >  1:2917318/1‑150 (MQ=255)
                                                                                                      gAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaag                                                <  2:551811/150‑1 (MQ=255)
                                                                                                       aaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaagg                                               >  2:2399602/1‑150 (MQ=255)
                                                                                                       aaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaagg                                               >  2:2180666/1‑150 (MQ=255)
                                                                                                       aaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaagg                                               >  1:4392450/1‑150 (MQ=255)
                                                                                                       aaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaagg                                               >  1:4248267/1‑150 (MQ=255)
                                                                                                         gTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaaggcc                                             >  1:87240/1‑150 (MQ=255)
                                                                                                           gcgcCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCAc                                           <  1:1538082/150‑1 (MQ=255)
                                                                                                            cgcCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCAcc                                          >  1:25326/1‑150 (MQ=255)
                                                                                                              cccGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGc                                        >  1:4823401/1‑150 (MQ=255)
                                                                                                                   gggTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATcc                                   <  2:4690102/150‑1 (MQ=255)
                                                                                                                     gTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGc                                 <  2:1204875/150‑1 (MQ=255)
                                                                                                                           ggCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGcc                           >  2:4515654/1‑150 (MQ=255)
                                                                                                                                gATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGa                      >  1:2816075/1‑150 (MQ=255)
                                                                                                                                gATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGa                      >  2:2520149/1‑150 (MQ=255)
                                                                                                                                 aTCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCGACAGGCCAAGGCCACCGCCCTCCGCAGCGCCAACGaa                     >  2:100955/1‑150 (MQ=255)
                                                                                                                                   cGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGc                   >  1:3996731/1‑150 (MQ=255)
                                                                                                                                   cGACCAGGGCGCCTTCACCGAAGCCGCCCAGGTCAAGACACGCGACCTGCTGGTCCAGACCGACCCGCGCCCGTTCCAGGCTGACGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGc                   <  2:4542266/150‑1 (MQ=255)
                                                                                                                                    gACCTGGGCGCCTTCGGTGAAGGCGACCTGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGcc                  >  1:4784305/1‑150 (MQ=255)
                                                                                                                                    gACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGcc                  >  1:3492148/1‑150 (MQ=255)
                                                                                                                                     aCCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGccc                 >  1:1764875/1‑150 (MQ=255)
                                                                                                                                     aCCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCACCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGccc                 <  2:750449/150‑1 (MQ=255)
                                                                                                                                       cAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgc               <  2:961342/150‑1 (MQ=255)
                                                                                                                                        aGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgcc              <  1:934707/150‑1 (MQ=255)
                                                                                                                                        aGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgcc              >  1:4696498/1‑150 (MQ=255)
                                                                                                                                        aGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgcc              >  1:3281798/1‑150 (MQ=255)
                                                                                                                                         ggTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgccg             >  1:880136/1‑150 (MQ=255)
                                                                                                                                         ggTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgccg             >  1:3662735/1‑150 (MQ=255)
                                                                                                                                          ggcGCCTTCGGTGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGt            >  1:2914/3‑150 (MQ=255)
                                                                                                                                          gTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGt            <  2:4887763/150‑1 (MQ=255)
                                                                                                                                             gCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGc         <  2:1488550/150‑1 (MQ=255)
                                                                                                                                               cTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGa       <  2:3335872/150‑1 (MQ=255)
                                                                                                                                               cTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGa       >  2:2206286/1‑150 (MQ=255)
                                                                                                                                               cTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCAACGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGa       >  1:3152758/1‑150 (MQ=255)
                                                                                                                                                ttCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAg      >  2:4761127/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGCTCCCGCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCC  >  NC_002947/3881995‑3882292

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