Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,320 A→G 7.0% intergenic (+348/‑172) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,3200AG 7.0% 547.9 / 21.2 188intergenic (+348/‑172)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (70/100);  new base G (7/6);  total (81/107)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

Only 100 of 188 total aligned reads displayed.

ATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTT  >  NC_002947/1899173‑1899468
                                                                                                                                                   |                                                                                                                                                     
         gggctgtcgacgttgacagcccgggatgaggtagatcagcaatgtatgggcccgctataccggccctttcgcgggcatGCCCGCTCCCCCAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                                                                                            >  2:5293311‑M1/79‑150 (MQ=255)
         aTTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCCCGCCCGCTCCCACAGGTACGGCGCCGCCCTCAGGCCGGGCGATATCCCGGTGGGAGCGGGCCTGCCCGCGaaa                                                                                                                                            >  1:268441/1‑150 (MQ=39)
         aTTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                                                                                            >  2:1420178/1‑150 (MQ=255)
               ccAGCAAAGCGCCTGGGAGCATGTGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCACTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                                                                                      <  1:306038/150‑1 (MQ=255)
                cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCCTGCCCGCGAAAGGGCCgg                                                                                                                                     >  2:2307825/1‑150 (MQ=255)
                     aaGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACCGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAg                                                                                                                                <  1:2014907/150‑1 (MQ=255)
                           cTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGcc                                                                                                                         >  2:1888404/1‑150 (MQ=255)
                                 gCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATaaa                                                                                                                   <  1:1327383/150‑1 (MQ=255)
                                   aTGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAAtt                                                                                                                 >  2:1153493/1‑150 (MQ=255)
                                    tGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTg                                                                                                                <  2:4867778/150‑1 (MQ=255)
                                      tggTTGCCACCGCATGGTGCTTAACCTGTGCAGGCCTCTTCGCGTGCACGCCCGCTCCCAATGTTACTGCGCCGCCATAAGGCCTTGCGATCACCCTGTGTGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                                                                              <  1:4275184/149‑1 (MQ=37)
                                      aGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                                                                              <  1:1241712/150‑1 (MQ=255)
                                      aGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                                                                              <  1:2212012/150‑1 (MQ=255)
                                        gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCTATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                                                                                            <  1:1165524/150‑1 (MQ=255)
                                        gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                                                                                            <  1:3915293/150‑1 (MQ=255)
                                        gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                                                                                            >  2:3536777/1‑150 (MQ=255)
                                        gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCCCGCCCCCAGGTACTGCGCCCCCCGCAGGCCGTGCGCTATACCGGGGGGGGCCGGCATGCCCCCCCAAGGGCCCGGGTAGCAGGGCCCAGAAAttctgg                                                                                                             >  1:638181/1‑146 (MQ=255)
                                         ttGCCACCGCATGGTGCGTAACCTGTGCAGGCCCCTTCGCGGGCACGCCCCCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCGGTGGGAGCGGGCATGCCCGCGCAAGGGCCGGTATAGCTGGCCCATAAATTGCTgtt                                                                                                           >  2:3796858/1‑148 (MQ=255)
                                             cACCGCTTGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                                                                                       <  1:25274/150‑1 (MQ=255)
                                             cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                                                                                       >  1:4438836/1‑150 (MQ=255)
                                             cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                                                                                       <  2:3021401/150‑1 (MQ=255)
                                             cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCGTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAAATGCTGATCTAc                                                                                                       >  2:258807/1‑150 (MQ=255)
                                              aCCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAcc                                                                                                      >  2:2881337/1‑150 (MQ=255)
                                                      gTGGGTTACCTGTGCAGGGGTCTTCGGGGGCTGGGCCGCTCCCACAGGTAGTGGGCCGCCCTCAGGCGTTGCGATATCGCTGTGGGAGCGGGCTTGCCCGCGAAAGGGCCGGTATAGCTGGCCCTTTAATTGCTGATCTTCCTCATCCCg                                                                                              <  1:1165391/150‑1 (MQ=255)
                                                       tGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCGTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCACATCCCgg                                                                                             >  1:3970837/1‑150 (MQ=255)
                                                        gCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCggg                                                                                            <  1:865126/150‑1 (MQ=255)
                                                            aaCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCGGTGGGAGCGGGCAtgccc                                                                                                                                                  >  1:1271785‑M2/1‑88 (MQ=255)
                                                             cccTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACTGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCCGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTc                                                                                       <  1:1773478/149‑1 (MQ=255)
                                                             aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTc                                                                                       <  2:2432511/150‑1 (MQ=255)
                                                             aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTc                                                                                       <  2:5160797/150‑1 (MQ=255)
                                                             aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTc                                                                                       <  1:2362361/150‑1 (MQ=255)
                                                             aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTc                                                                                       <  2:4438836/150‑1 (MQ=255)
                                                               cTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACTGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCaa                                                                                     <  1:4993757/150‑1 (MQ=255)
                                                                tgtgCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTTTCTGCGCCGCCGTCAGGCCTTGCTATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAAc                                                                                    <  1:4644251/150‑1 (MQ=255)
                                                                   gCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTc                                                                                 <  1:5127039/150‑1 (MQ=255)
                                                                   gCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTc                                                                                 <  1:2195249/150‑1 (MQ=255)
                                                                     aGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCGTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGa                                                                               >  2:2722657/1‑150 (MQ=255)
                                                                      gccgtctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATTTCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGAc                                                                              <  1:1859453/146‑1 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                                                                                            >  2:1540639/1‑88 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                                                                                            <  1:1540639/88‑1 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCAGGGCTGTCGACGTTGACa                                                                             <  2:3479868/150‑1 (MQ=255)
                                                                            ttCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCt                                                                        <  1:979561/150‑1 (MQ=255)
                                                                              cgcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGCTATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCGGTCAACGTCGACAGCCCTca                                                                      >  2:4698138/1‑150 (MQ=255)
                                                                                 gggcatgcccgctcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcacca                                                                   >  2:1251212‑M2/67‑150 (MQ=255)
                                                                                  ggCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcaccac                                                                  >  2:1254821/1‑150 (MQ=255)
                                                                                       GCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                             <  2:2909966‑M1/150‑1 (MQ=255)
                                                                                        CCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACg                                                            <  2:2515592‑M1/150‑1 (MQ=255)
                                                                                        CCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACATCACg                                                            <  1:3729294‑M1/150‑1 (MQ=255)
                                                                                         ccgctcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGa                                                           >  1:3082833‑M2/59‑150 (MQ=255)
                                                                                         ccGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCAtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagca                                                            <  1:5344665‑M2/150‑92 (MQ=255)
                                                                                           gCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCAtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcat                                                                                                    >  1:2545266‑M2/1‑57 (MQ=255)
                                                                                           gCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCAtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcat                                                                                                    <  2:2545266‑M2/108‑52 (MQ=255)
                                                                                             tcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                       >  1:3456339‑M2/55‑150 (MQ=255)
                                                                                             tcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                       >  1:4076311‑M2/55‑150 (MQ=255)
                                                                                             tCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCAtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcctgagcactgc                                                        <  1:4864440‑M2/150‑96 (MQ=255)
                                                                                              cccACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCGtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagcactgcc                                                       <  1:1345713‑M2/150‑97 (MQ=255)
                                                                                               ccACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCGtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagcactgccc                                                      <  1:2476208‑M2/150‑98 (MQ=255)
                                                                                               CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTACACGACAAAgg                                                     <  1:2916062‑M1/150‑1 (MQ=255)
                                                                                               CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCACACCACTTCACGACAAAgg                                                     <  1:186220‑M1/150‑1 (MQ=255)
                                                                                               CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCACACCACATCACGACAAagc                                                     <  1:4222144‑M1/150‑2 (MQ=255)
                                                                                                 ACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGAAAAAGgcg                                                   <  2:4586288‑M1/150‑1 (MQ=255)
                                                                                                  CAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgcgc                                                  >  2:5218582‑M1/1‑150 (MQ=255)
                                                                                                  CAGGTACTGCGCCGCCCTCAGGCCGTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgcgc                                                  >  2:773143‑M1/1‑150 (MQ=255)
                                                                                                   agggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGcc                                                 >  1:4951569‑M2/49‑150 (MQ=255)
                                                                                                   aGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCGtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagcactgccctccc                                                  <  1:4164910‑M2/150‑102 (MQ=255)
                                                                                                   AGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGt                                                                                        >  1:5200692‑M1/1‑111 (MQ=255)
                                                                                                   AGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGcc                                                 >  1:1420178‑M1/1‑150 (MQ=255)
                                                                                                   AGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCAAAAATTGCTGATCTAACTCATCCCGGGGTGGCGACGTCGAGCGCCCTCACCACATCACGACAAAGGCGcc                                                 <  1:4770742‑M1/150‑1 (MQ=255)
                                                                                                     ggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCgt                                               >  2:3696885‑M2/47‑150 (MQ=255)
                                                                                                     GTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCGCACCACTTCACGACAAAGGCGCcgg                                               <  1:178357‑M1/150‑2 (MQ=255)
                                                                                                       aCTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCgtgt                                             <  1:1378661/150‑1 (MQ=255)
                                                                                                        cTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGtt                                            <  1:1254748/150‑1 (MQ=255)
                                                                                                        cTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGtt                                            <  2:5228748/150‑1 (MQ=255)
                                                                                                            gccgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaa                                                        >  2:4404528/1‑134 (MQ=255)
                                                                                                            gccgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaa                                                        <  1:4404528/134‑1 (MQ=255)
                                                                                                            gccgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGccc                                        >  2:5380845/1‑150 (MQ=255)
                                                                                                              cgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCAtgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagcactgccctcccatgtcccagga                                       >  2:2983645‑M2/1‑38 (MQ=255)
                                                                                                                cccTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCt                                    <  1:3419838/150‑1 (MQ=255)
                                                                                                                cccTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCt                                    <  1:3109215/150‑1 (MQ=255)
                                                                                                                 gcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAAt                                                                                                                  <  2:233814‑M2/37‑1 (MQ=255)
                                                                                                                  cTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTg                                  >  2:5016763/1‑150 (MQ=255)
                                                                                                                         cTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGc                           <  2:4507002/150‑1 (MQ=255)
                                                                                                                          ttGCGATATCCCTGTGGGAGCGGGCCTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCGGTCAACGTCGACAGCCCTCACCACGTCACGACAAAGGCGCCGGGTTGCCCCGCTTgcagcagca                          >  1:106011/1‑150 (MQ=255)
                                                                                                                          ttGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGACGACAGCCCACACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTAGCATCAGCa                          >  1:5324622/1‑150 (MQ=255)
                                                                                                                               atatCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCg                     <  1:2850586/150‑1 (MQ=255)
                                                                                                                                 gtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCggg                   <  1:3696885‑M2/132‑1 (MQ=255)
                                                                                                                                  tCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCgggg                  <  2:966244/150‑1 (MQ=255)
                                                                                                                                  tCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCGCACCACGTCACGACAAAGGCGCCGGGTTGCCCCGCGAGCATCAGCATGCCgggg                  >  1:3006402/1‑150 (MQ=255)
                                                                                                                                      tgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCgt                                               >  2:825508‑M2/14‑117 (MQ=255)
                                                                                                                                      tgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGc              >  2:3769583‑M2/14‑150 (MQ=255)
                                                                                                                                       gtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGAcagcccggg                                                                      >  1:3234837‑M2/13‑90 (MQ=255)
                                                                                                                                         tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCgg           >  2:1336165‑M2/11‑150 (MQ=255)
                                                                                                                                         gggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCgg           <  2:3970837‑M2/140‑1 (MQ=255)
                                                                                                                                          ttcgcggGCAcgcccgcgaagaggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagcactgccctcccatgtcccaggagatgccggtgatcgaggacgggtttgag           <  2:602176‑M2/147‑141 (MQ=255)
                                                                                                                                          ggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGgc          <  2:3236153‑M2/141‑1 (MQ=255)
                                                                                                                                           gagcgggcATGCCCGCGAAAGGGCCGGTATAGATGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGgcg         >  2:3330626‑M2/9‑150 (MQ=255)
                                                                                                                                              cgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGa                                                           >  2:3320699‑M2/6‑97 (MQ=255)
                                                                                                                                              cgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGa                                                           <  1:3320699‑M2/92‑1 (MQ=255)
                                                                                                                                              cgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGa                                                           >  2:595332‑M2/6‑97 (MQ=255)
                                                                                                                                              cGGGCAtgcccgcgaaagggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctggggcaatttccagcatgagcactgccctcccatgtcccaggagatgccggtgatcgaggacgggtttgaggtgc       >  2:3110885‑M2/1‑6 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTT  >  NC_002947/1899173‑1899468

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