Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 726,960 C→A 5.5% intergenic (‑204/+63) PP_0620 ← / ← PP_0621 GntR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947726,9600CA 5.5% 283.9 / 3.0 91intergenic (‑204/+63)PP_0620/PP_0621GntR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/44);  new base A (3/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GTCAATGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGCCCGGTCGTTGGGCAG  >  NC_002947/726815‑727109
                                                                                                                                                 |                                                                                                                                                     
gTCAATGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTcg                                                                                                                                                   >  2:1384180/1‑150 (MQ=255)
  cAATGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTcgcg                                                                                                                                                 <  1:1810875/150‑1 (MQ=255)
   aaTGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                <  1:349535/150‑1 (MQ=255)
   aaTGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                <  1:2094455/150‑1 (MQ=255)
   aaTGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                <  1:2054250/150‑1 (MQ=255)
        tGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTTGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCCTACATGTTGGTGCCTGCTTCGCGGGCAAg                                                                                                                                           >  1:902048‑M2/1‑149 (MQ=255)
         ggCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgc                                                                                                                                          >  1:1384180‑M2/1‑148 (MQ=255)
         ccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctctgtgggagcgggcacGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGggcatgc                                                                                                                                          >  1:1253263‑M1/92‑147 (MQ=255)
           cGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgccc                                                                                                                                        >  2:231544‑M2/1‑146 (MQ=255)
           aggggggcaGGGCGGGTAAATGGGGTGAAAGGGTCGGGGAGGGATATGTGGCACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgccc                                                                                                                                        >  1:2087241‑M2/2‑146 (MQ=255)
             ggggggAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgcccgc                                                                                                                                      <  1:1842838‑M2/150‑7 (MQ=255)
               gaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccccctgtgggagcgggcacGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTc                                                                                                                                    <  2:234861‑M1/65‑1 (MQ=255)
                gggAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgcccgctcc                                                                                                                                   >  2:971280‑M2/1‑141 (MQ=255)
                  gAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCggg                                                                                                                                               >  1:57709/1‑136 (MQ=255)
                  gAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCggg                                                                                                                                               <  2:57709/136‑1 (MQ=255)
                    ggggTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCTATACACCGTGGTGCCTGCTTCGCGGGCAAGATCGCTATcaca                                                                                                                               <  2:1258115/150‑1 (MQ=38)
                       tggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccccctgtgggagcgggcacGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAgag                                                                                                                            >  1:234861‑M1/78‑149 (MQ=255)
                       ggcGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACCAAATTGCAGATAAAAACTGTGGGAGCGGGCGCGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGAt                                                                                                                            <  2:404937/148‑1 (MQ=255)
                       gTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGAt                                                                                                                            >  2:1366515/1‑150 (MQ=255)
                       gTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCAGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGAGCAAGCCGGCTCCCACAGaa                                                                                                                            >  2:2196977/1‑149 (MQ=255)
                        tCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGTCGTAGTCATACACCATGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGAtt                                                                                                                           >  1:1310384/1‑150 (MQ=38)
                            aaaTGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGa                                                                                                                       >  2:1828852/1‑150 (MQ=255)
                             aaTGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCGCCGAGCACCGGCGAAGCCGGTGCCATACACCGAGGTGCCTGCTTCGCGGGCAAGCGCGCTCCCACAGATTTCgca                                                                                                                      >  1:154273/1‑148 (MQ=37)
                              aTGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAAt                                                                                                                     >  2:2273955/1‑150 (MQ=255)
                              aTGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAAt                                                                                                                     <  2:1553122/150‑1 (MQ=255)
                               tGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAAtt                                                                                                                    <  1:769475/150‑1 (MQ=255)
                                ggggtgAAAGGGTCGGGGAGGTATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTg                                                                                                                   <  2:1977792/147‑1 (MQ=255)
                                 gtgtGAAAGGGTCGGGGAGGGATATGTGGAACACGAATTTGCAGATAAAAACTGTGGGGGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACAGTGGTGCCTGCTTCGCGGGCAAGCCAGCTCACACAGATTTCGAATTgg                                                                                                                  >  1:1888307/1‑150 (MQ=37)
                                  tgggAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGt                                                                                                                 <  2:1894022/147‑1 (MQ=255)
                                      aaaGGGTCGGGGGGGGATATGTGGAACACGAACTTGCAGATAAAAACTGTGGGGACGGGCGTGCCCGAGAACAACGGCGAAGCCGGTGCCATACACCTTGGTGCCTGCTTCGCGGGCAAGCGCGCTCCCACAGATTTCGAATCAGtttga                                                                                                             >  2:2279455/1‑148 (MQ=255)
                                      aaaGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCa                                                                                                             <  1:2273955/150‑1 (MQ=255)
                                       gggggTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAg                                                                                                            >  2:1788358/3‑150 (MQ=255)
                                       aaGGGTCGGGTAGTGATATGTTGAACACGAAATTGCAGATAAAAACTGTGTGAGCCGGCGTTCCCGCTAACACCAGCGAATCCGTTGCTATTCACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAg                                                                                                            <  2:780362/150‑1 (MQ=37)
                                          aagtaacaggtaccacacagtccgcaagcacgccgcgcttccctgtgggagcgggcacGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGgcatg                                                                                                                                           <  2:449511‑M1/58‑3 (MQ=255)
                                          aagtaacaggtaccacacagtccgcaagcacgccgcgcttccctgtgggagcgggcacGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGgcatg                                                                                                                                           >  1:449511‑M1/59‑114 (MQ=255)
                                           gTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGAtt                                                                                                        <  2:2284149/150‑1 (MQ=255)
                                           gTCGGGGAGGGATATGTGGAACAAGACATTGCAGATAAAAACTGGGGGGGCGGGCGTGCCCTCGAACACCGGCGAAGCCGGTGCCATACACCATGGTGCCTGCCTCGCGCGCAAGCCCGCTCCCACAGATTTCTAATTCGTTTCAGGAtt                                                                                                        <  2:1310384/150‑1 (MQ=37)
                                                 gagggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                >  2:1944938/1‑104 (MQ=255)
                                                 gagggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                <  1:1944938/104‑1 (MQ=255)
                                                 gagggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGa                                                                                                          >  2:1828346/1‑142 (MQ=255)
                                                 gagggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGa                                                                                                          <  1:1828346/142‑1 (MQ=255)
                                                  agggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGa                                                                                                 <  1:1828852/150‑1 (MQ=255)
                                                   gggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGAc                                                                                                <  1:2102639/150‑1 (MQ=255)
                                                    ggaTATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACa                                                                                               >  2:1416592/1‑150 (MQ=255)
                                                          gtgGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGCGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGaa                                                                                                                      >  2:1875624/1‑121 (MQ=255)
                                                          gtgGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGCGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGaa                                                                                                                      <  1:1875624/121‑1 (MQ=255)
                                                           tgGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCgtg                                                                                          >  1:2029595/1‑148 (MQ=255)
                                                           tgGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCgtg                                                                                          <  2:2029595/148‑1 (MQ=255)
                                                            ggAACACGAAATTGAAGATAGCAACAGTGGTTGTGGTTGTGCCCGCGATCACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTc                                                                                       >  1:1221171/1‑150 (MQ=37)
                                                              aaCACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                >  1:2087041/1‑91 (MQ=255)
                                                              aaCACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCgg                                                                                                                                                <  2:2087041/91‑1 (MQ=255)
                                                              aaCACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAg                                                                                                            <  2:1868888/127‑1 (MQ=255)
                                                              aaCACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAg                                                                                                            >  1:1868888/1‑127 (MQ=255)
                                                               acacGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGc                                                                                    <  2:289617/150‑1 (MQ=255)
                                                                 acGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCt                                                                                  >  1:767572/1‑150 (MQ=255)
                                                                   gAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtg                                                                                <  2:1187300/150‑1 (MQ=255)
                                                                   gAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtg                                                                                <  2:1187196/150‑1 (MQ=255)
                                                                   gAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtg                                                                                >  1:617571/1‑150 (MQ=255)
                                                                   gAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtg                                                                                <  1:1366515/150‑1 (MQ=255)
                                                                     aaTTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtggt                                                                              <  2:617571/150‑1 (MQ=255)
                                                                             aTAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAggg                                                                      <  1:1707964/150‑1 (MQ=255)
                                                                                       tgGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGa                                                                                                                       <  1:1816473/91‑1 (MQ=255)
                                                                                       tgGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGa                                                                                                                       >  2:1816473/1‑91 (MQ=255)
                                                                                       tgGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGGGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCGCTGCTGCa                                                            >  1:1117834/1‑150 (MQ=255)
                                                                                        gggAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAg                                                           <  1:1477631/150‑1 (MQ=255)
                                                                                        gggAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAg                                                           <  1:656188/150‑1 (MQ=255)
                                                                                           agtgggAGTGCCCGCTATCACCTGCGTAGCCGGTGCGATACCCCTTGGTTCCTGCTTCGCGGGCAAGCCCTCTCCCACAGAGTCCGAATTGCTTTCAGGAGTGAACAGACATCGTGTTCAGCCTTGGTCGCCAGGGCCCTGCTGCAAGTg                                                        <  2:1888307/147‑1 (MQ=255)
                                                                                             cGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgcccgctcccacagtttttatctgcaatttcgtgttccacatatccctccccgaccctttcacaccatttacgacccctccccccc                                                      <  1:1034984‑M2/150‑87 (MQ=255)
                                                                                              gggCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtggtg                                                                             <  2:389073/126‑1 (MQ=255)
                                                                                              gggCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTtggtg                                                                             >  1:389073/1‑126 (MQ=255)
                                                                                                  gTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgcccgctcccacagtttttatctgcaatttcgtgttccacatatccctccccgaccctttcacaccatttacgacccctcccccccgccag                                                 >  2:1034984‑M2/1‑59 (MQ=255)
                                                                                                  gTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCctgc                                                                <  2:1979883/135‑1 (MQ=255)
                                                                                                  gTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCctgc                                                                >  1:1979883/1‑135 (MQ=255)
                                                                                                      ccGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGc                                             >  2:958604/1‑150 (MQ=255)
                                                                                                          gAACACCGGCGAAGCCGGTGCCATACCCCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATa                                         <  1:2196977/150‑1 (MQ=255)
                                                                                                               ccGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAgcccgctcccacagtttttatctgcaatttcgtgttccacatatccctccccgaccctttcacaccatttacgacccctcccccccgccagcattgacccgtcc                                    <  1:1659619‑M2/150‑105 (MQ=255)
                                                                                                                cGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTc                                   >  2:1893352/1‑150 (MQ=255)
                                                                                                                  gCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCAt                                 >  2:1091288/1‑150 (MQ=255)
                                                                                                                    gAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTg                               >  2:1498408/1‑150 (MQ=255)
                                                                                                                           gTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCctgct                        <  2:2331674/150‑1 (MQ=255)
                                                                                                                           gTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCctgct                        >  1:1902995/1‑150 (MQ=255)
                                                                                                                           gTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCct                           >  1:641977/1‑147 (MQ=255)
                                                                                                                           gTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCct                           <  2:641977/147‑1 (MQ=255)
                                                                                                                            tGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGTCCTGCTGCAGGTGGGCCTGGCTGCAATACCACACATTGCCgtgcgg                       <  2:1842838/150‑3 (MQ=255)
                                                                                                                                aTACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAgc                   >  1:497143/1‑150 (MQ=255)
                                                                                                                                    acCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGccc               >  1:981911/1‑150 (MQ=255)
                                                                                                                                        tggtgCCTGCTTCGCGGGCAAGCCCTCTCCCACAGATTTTGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGCCCGGTc           <  2:154273/150‑1 (MQ=255)
                                                                                                                                             ccTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGCCCGGTCGTTgg      >  2:49967/1‑150 (MQ=255)
                                                                                                                                               tGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCa                                                                                                             >  1:74057/1‑45 (MQ=255)
                                                                                                                                               tGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCa                                                                                                             <  2:74057/45‑1 (MQ=255)
                                                                                                                                               tGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGCCCGGTCGTTGGGc    >  1:594180/1‑150 (MQ=255)
                                                                                                                                                 cTTCGCGGGCAAGCCCGCTCCCCCAGATTTCGAATTGGTTTCAGGAGTTATCAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGCCCGGTCGTTGGGCAg  <  1:841106/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTCAATGCTGGCGGGGGGGAGGGGTCGTAAATGGTGTGAAAGGGTCGGGGAGGGATATGTGGAACACGAAATTGCAGATAAAAACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATACACCGTGGTGCCTGCTTCGCGGGCAAGCCCGCTCCCACAGATTTCGAATTGGTTTCAGGATTGAACAGACATCGTGTTCAGCCTTGGTGGCCAGGGCCCTGCTGCAGGTGGGCCTGGCTGCAATACCACTCATTGCCCTGCTGCAGCGCCCGGTCGTTGGGCAG  >  NC_002947/726815‑727109

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