Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,022,851 C→G 6.4% intergenic (‑502/+166) PP_3547 ← / ← emrB NAD(P)‑binding oxidoreductase/multidrug transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,022,8510CG 6.4% 298.7 / 9.0 94intergenic (‑502/+166)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base G (3/3);  total (47/47)
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.57e-01

TTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATC  >  NC_002947/4022702‑4023000
                                                                                                                                                     |                                                                                                                                                     
tttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATcc                                                                                                                                                       <  1:1774889/150‑1 (MQ=255)
 ttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCa                                                                                                                                                      <  1:1516918/150‑1 (MQ=255)
  tGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCAt                                                                                                                                                     <  1:2094422/150‑1 (MQ=255)
    tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGc                                                                                                                                                   <  2:1173378/150‑1 (MQ=255)
    tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGc                                                                                                                                                   >  1:28508/1‑150 (MQ=255)
      aaaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGccc                                                                                                                                                 >  1:813834/1‑150 (MQ=255)
       aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGCGGACCACGCTGCAGACCTCATGCCGGGGCCCGGCGCTGGACGGGCATGGATGCCCg                                                                                                                                                >  2:265188/1‑150 (MQ=17)
         cgctaCAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                              <  2:15473/146‑1 (MQ=255)
           cgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                            <  1:265188/150‑2 (MQ=255)
               aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTcc                                                                                                                                        <  1:2249519/150‑1 (MQ=255)
                  aaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtt                                                                                                                                     >  2:1607269/1‑150 (MQ=255)
                       aCCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                                                <  1:831327/150‑1 (MQ=255)
                         cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGc                                                                                                                              >  1:798751/1‑150 (MQ=255)
                         cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGc                                                                                                                              >  2:2013994/1‑150 (MQ=255)
                          aaaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGc                                                                                                                              >  1:1571331/1‑149 (MQ=255)
                          aaaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGc                                                                                                                              <  2:1571331/149‑1 (MQ=255)
                              aaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTgc                                                                                                                         >  2:1540239/1‑150 (MQ=255)
                                   tGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCCAGCCCTGCGCTCg                                                                                                                    >  1:267599/1‑150 (MQ=255)
                                    ggCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCgg                                                                                                                   <  1:1837342/150‑1 (MQ=39)
                                        ccTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcctc                                                                                                               <  1:1916354/150‑1 (MQ=38)
                                              atcCGGCCCCTGCGCTACGCCCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTccc                                                                                                                                            <  2:646729/115‑3 (MQ=18)
                                              atcCGGCCCCTGCGCTACGCCCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                                                                            >  1:646729/1‑115 (MQ=18)
                                               tcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                                                >  2:67860/1‑126 (MQ=37)
                                               tcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                                                <  1:67860/126‑1 (MQ=37)
                                                ccGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTc                                                                                                       <  2:2115154/150‑1 (MQ=35)
                                                     cccTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAAt                                                                                                  >  2:445688/1‑150 (MQ=35)
                                                         gcgcTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCcgc                                                                                              <  1:445688/150‑1 (MQ=38)
                                                             tACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCgcggcg                                                                                          >  2:495394/1‑150 (MQ=40)
                                                                  tCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTgg                                                                                     <  1:495394/150‑1 (MQ=40)
                                                                   ccGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTggg                                                                                    <  1:2013994/150‑1 (MQ=255)
                                                                          tCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgc                                                                             <  1:1607269/150‑1 (MQ=255)
                                                                            ccTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgc                                                                           >  2:905819/1‑150 (MQ=255)
                                                                            ccTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgc                                                                           <  2:1020232/150‑1 (MQ=255)
                                                                             cTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgcg                                                                          <  2:1793384/150‑1 (MQ=255)
                                                                             cTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgcg                                                                          <  1:367623/150‑1 (MQ=255)
                                                                              tCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCaa                                                                                                   <  1:2401511/124‑1 (MQ=33)
                                                                              tCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCaa                                                                                                   >  2:2401511/1‑124 (MQ=33)
                                                                              tCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCgcgcgc                                                                         >  1:1392053/1‑150 (MQ=255)
                                                                              tCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCgcgcgc                                                                         >  1:2265326/1‑150 (MQ=255)
                                                                                 cTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTa                                                                      >  1:727330/1‑150 (MQ=255)
                                                                                    cGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                                                     <  1:2125064/100‑1 (MQ=34)
                                                                                    cGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                                                     >  2:2125064/1‑100 (MQ=34)
                                                                                        cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaa                                                               >  2:886154/1‑150 (MQ=255)
                                                                                          tGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaaaa                                                             <  2:1912257/150‑1 (MQ=255)
                                                                                            cTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATc                                                           <  2:716905/150‑1 (MQ=255)
                                                                                             tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCg                                                                                                      >  1:1157810/1‑106 (MQ=33)
                                                                                             tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCg                                                                                                      <  2:1157810/106‑1 (MQ=32)
                                                                                               ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaaaa                                                             <  1:1857258/145‑1 (MQ=255)
                                                                                               ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaaaa                                                             >  2:1857258/1‑145 (MQ=255)
                                                                                               ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCgaa                                                        >  1:1474263/1‑150 (MQ=255)
                                                                                               ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaa                                                        >  1:825448/1‑150 (MQ=255)
                                                                                                cGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCgcggc                                                                                           >  1:1003032/1‑114 (MQ=40)
                                                                                                cGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCgcggc                                                                                           <  2:1003032/114‑1 (MQ=40)
                                                                                                cGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                       <  2:2044226/150‑1 (MQ=255)
                                                                                                  ggaggaCCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGc                                                     <  2:813834/150‑1 (MQ=255)
                                                                                                   gaggaCCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGcc                                                    >  1:159324/1‑150 (MQ=255)
                                                                                                        ccGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagc                                               >  2:381880/1‑150 (MQ=255)
                                                                                                        ccGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGACGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagc                                               <  2:2143763/150‑1 (MQ=255)
                                                                                                         cgcgCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcg                                              >  2:1196074/1‑150 (MQ=255)
                                                                                                            gcTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgaga                                           <  2:159324/150‑1 (MQ=255)
                                                                                                            gcTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgaga                                           >  1:2079735/1‑150 (MQ=255)
                                                                                                             cTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagag                                          <  2:526987/150‑1 (MQ=255)
                                                                                                                tgggCCCCAGGATGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGTAAAATCAAAAGCCGAagcgagagcga                                       >  1:2186149/3‑150 (MQ=255)
                                                                                                                 aGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgag                                      >  1:472974/1‑150 (MQ=255)
                                                                                                                   gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagag                                    >  1:677264/1‑150 (MQ=255)
                                                                                                                     cccATCCTGGGGCCCAGCGCTTGCCGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgaga                                     <  2:1508388/147‑1 (MQ=255)
                                                                                                                     cccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcctc                                                                                                               >  2:1784641/1‑73 (MQ=35)
                                                                                                                     cccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcctc                                                                                                               <  1:1784641/73‑1 (MQ=35)
                                                                                                                     cccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcg                                  >  2:1924603/1‑150 (MQ=255)
                                                                                                                     cccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgaga                                     >  1:1508388/1‑147 (MQ=255)
                                                                                                                     ccaaTCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTg                                                                                      >  1:2004509/4‑98 (MQ=255)
                                                                                                                     ccaaTCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTg                                                                                      <  2:2004509/95‑1 (MQ=255)
                                                                                                                      ccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTGAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcga                                 <  2:1392053/150‑1 (MQ=255)
                                                                                                                       cATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGTGCTATCGCGCGCTTAGAAGCAAAATCAAAAGTCGAagcgagagcgagagcgag                                >  2:1238864/1‑150 (MQ=255)
                                                                                                                       cATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgag                                >  1:1297726/1‑150 (MQ=255)
                                                                                                                        aTCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTTGCAATCCGCGGCGCCTGGGCTATCGCGCGCGTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgaga                               >  1:1717384/1‑150 (MQ=255)
                                                                                                                        aTCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgaga                               <  1:381880/150‑1 (MQ=255)
                                                                                                                        aTCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgaga                               >  2:2281825/1‑150 (MQ=255)
                                                                                                                        aTCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgaga                               <  1:73984/150‑1 (MQ=255)
                                                                                                                          ccTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGGAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagc                             >  2:410613/1‑150 (MQ=255)
                                                                                                                            tGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcga                           <  2:727330/150‑1 (MQ=255)
                                                                                                                              gggCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgaga                         >  2:1895711/1‑150 (MQ=255)
                                                                                                                                    aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCACTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgaga                   >  1:1213699/1‑150 (MQ=255)
                                                                                                                                      cgcTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgagagc                 <  1:1533675/150‑1 (MQ=255)
                                                                                                                                          tGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcga                                             <  2:2054043/118‑1 (MQ=255)
                                                                                                                                          tGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcga                                             >  1:2054043/1‑118 (MQ=255)
                                                                                                                                           gACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGcgagagcgagagcgagagcgagagcgagagcgagag            <  1:1196074/150‑1 (MQ=255)
                                                                                                                                            aCGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgaga                         <  2:2375393/136‑1 (MQ=255)
                                                                                                                                            aCGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgaga                         >  1:2375393/1‑136 (MQ=255)
                                                                                                                                              gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCAt         <  2:677264/150‑1 (MQ=255)
                                                                                                                                              gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCAt         >  2:1561765/1‑150 (MQ=255)
                                                                                                                                                gCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCg       <  1:1561765/150‑1 (MQ=255)
                                                                                                                                                  aTCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                       <  1:1343086/100‑1 (MQ=255)
                                                                                                                                                  aTCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                       >  2:1343086/1‑100 (MQ=255)
                                                                                                                                                   tCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGa    <  1:1895711/150‑1 (MQ=255)
                                                                                                                                                    ccATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGAt   <  2:321625/150‑1 (MQ=255)
                                                                                                                                                     cATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATc  <  2:798751/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATC  >  NC_002947/4022702‑4023000

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