Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,320 A→G 8.5% 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 8.5% 263.5 / 9.0 93intergenic (+348/‑172)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (45/40);  new base G (6/2);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTT  >  NC_002947/1899174‑1899468
                                                                                                                                                   |                                                                                                                                                    
tGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                                                                                   >  2:150691/1‑150 (MQ=255)
tGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                                                                                   <  2:958911/150‑1 (MQ=255)
tGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                                                                                   <  2:1099405/150‑1 (MQ=255)
         ttGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAg                                                                                                                                          <  1:1110947/150‑1 (MQ=255)
          gatctcAGCAAAGCGCCTGGGAGCATGAGGTTGACATGAATGGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAgg                                                                                                                                         <  1:1143163/145‑1 (MQ=38)
             cccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                                                                                      >  1:781260/1‑150 (MQ=255)
             cccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                                                                                      >  2:1391905/1‑150 (MQ=255)
               cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                                                                                     >  2:808435/1‑149 (MQ=255)
               cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                                                                                     <  1:808435/149‑1 (MQ=255)
                aGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCCTGCCGCGTAGCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGt                                                                                                                                   <  2:1403982/150‑1 (MQ=255)
                aGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGTAGCGGGCATTCCCGCTAAAgggagatc                                                                                                                                   >  2:1143163/1‑145 (MQ=39)
                 gCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGta                                                                                                                                  >  1:516764/1‑150 (MQ=255)
                        gcCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTgg                                                                                                                           <  2:628427/150‑1 (MQ=255)
                           tGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGccc                                                                                                                        >  2:428168/1‑150 (MQ=255)
                            gggAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                                                                                   <  1:502939/122‑1 (MQ=255)
                            gggAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                                                                                   >  2:502939/1‑122 (MQ=255)
                              gAGCATGAGGTTGCCACCGCATTGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCACGCTCCCACAGGTACTTCCCCGCCCTCAGGCCTTGCAATATCACTGTGTTAGCGGGAATGCCCGCGAAAGGGACCGTGTAGCTGACCCATa                                                                                                                     >  2:346656/1‑150 (MQ=255)
                                gCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATaaa                                                                                                                   <  1:139303/150‑1 (MQ=255)
                                   tGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTg                                                                                                                >  1:1353443/1‑150 (MQ=255)
                                    gAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                                                                                           >  2:818213/1‑122 (MQ=255)
                                    gAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                                                                                           <  1:818213/122‑1 (MQ=255)
                                     aGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                                                                              <  1:43270/150‑1 (MQ=255)
                                       gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                                                                                            <  2:634995/150‑1 (MQ=255)
                                       gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                                                                                            <  1:493317/150‑1 (MQ=255)
                                       gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                                                                                            >  2:493317/1‑150 (MQ=255)
                                       gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCAGGTATCGCTGGCCCATAAATTGCCGa                                                                                                            <  2:1157334/150‑1 (MQ=255)
                                         tGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATc                                                                                                          >  2:970239/1‑150 (MQ=255)
                                            cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                                                                                       >  2:1300942/1‑150 (MQ=255)
                                            cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                                                                                       >  2:208737/1‑150 (MQ=255)
                                              ccGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCt                                                                                                     >  2:1281238/1‑150 (MQ=255)
                                                            aaCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACTGCGCCGCGCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGt                                                                                        >  1:355134/1‑150 (MQ=255)
                                                             aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                                                                                      >  1:669775/1‑103 (MQ=255)
                                                             aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                                                                                      <  2:669775/103‑1 (MQ=255)
                                                               cTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCaa                                                                                     >  1:457690/1‑150 (MQ=255)
                                                                      ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCCCCGTCGAc                                                                              >  1:922920/1‑150 (MQ=255)
                                                                      ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGAc                                                                              <  1:150691/150‑1 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACa                                                                             >  2:701827/1‑150 (MQ=255)
                                                                         ctctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                                                                              >  2:1052790/1‑115 (MQ=255)
                                                                         ctctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                                                                              <  1:1052790/115‑1 (MQ=255)
                                                                          tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATGCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                          >  1:1243524/1‑150 (MQ=255)
                                                                          tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGcc                                                                                                                                                  >  2:1259795/1‑78 (MQ=255)
                                                                          tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGcc                                                                                                                                                  <  1:1259795/78‑1 (MQ=255)
                                                                          tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                          <  1:939116/150‑1 (MQ=255)
                                                                          tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                          >  1:955013/1‑150 (MQ=255)
                                                                          tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                          <  1:208737/150‑1 (MQ=255)
                                                                              cgcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTca                                                                      <  2:516764/150‑1 (MQ=255)
                                                                               gcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTGGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcac                                                                     <  1:1284444/150‑1 (MQ=255)
                                                                                     aCGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGTCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTc                                                               >  1:249182/1‑150 (MQ=255)
                                                                                             tcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                       >  2:505496‑M2/55‑150 (MQ=255)
                                                                                             tcccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                       >  2:855685‑M2/55‑150 (MQ=255)
                                                                                             TCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAgcccg                                                                        >  1:1147881‑M1/1‑132 (MQ=255)
                                                                                             TCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAgcccg                                                                        <  2:1147881‑M1/133‑2 (MQ=255)
                                                                                             TCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                       <  1:1019109‑M1/150‑1 (MQ=255)
                                                                                              cccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAg                                                      >  2:1421273‑M2/54‑150 (MQ=255)
                                                                                              CCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAg                                                      >  2:955013‑M1/1‑150 (MQ=255)
                                                                                               ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                     >  1:1153826‑M2/53‑150 (MQ=255)
                                                                                               ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                     >  2:308284‑M2/53‑150 (MQ=255)
                                                                                               ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                              >  2:106841‑M2/53‑141 (MQ=255)
                                                                                               ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                              <  1:106841‑M2/89‑1 (MQ=255)
                                                                                               ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                     >  2:1119965‑M2/53‑150 (MQ=255)
                                                                                               ccacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                     >  2:1001149‑M2/53‑150 (MQ=255)
                                                                                               CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                     <  1:412811‑M1/150‑1 (MQ=255)
                                                                                               CCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAgg                                                     <  1:57793‑M1/150‑1 (MQ=255)
                                                                                                cacagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgcg                                                   >  1:187775‑M2/52‑150 (MQ=255)
                                                                                                 acagggatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTGCCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgcg                                                   >  2:406547‑M2/51‑150 (MQ=255)
                                                                                                   AGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAAc                                                                                    <  1:1246335‑M1/115‑1 (MQ=255)
                                                                                                   AGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAAc                                                                                    >  2:1246335‑M1/1‑115 (MQ=255)
                                                                                                        cTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGtt                                            >  2:93056/1‑150 (MQ=255)
                                                                                                            gccgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGccc                                        >  1:1239510/1‑150 (MQ=255)
                                                                                                             caaggcctgagggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGcccc                                       >  1:286496‑M2/39‑150 (MQ=255)
                                                                                                              aaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCg                                      >  1:831982‑M2/38‑150 (MQ=255)
                                                                                                              aaggcctgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCg                                      >  2:1160777‑M2/38‑150 (MQ=255)
                                                                                                                 ccTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCtt                                   <  2:781260/150‑1 (MQ=255)
                                                                                                                  cTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTg                                  >  2:368548/1‑150 (MQ=255)
                                                                                                                    tgagggcggcgcagtacctgtgggagcgggcGTGCCCGCGAAGAGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCa                                <  2:1153826‑M2/119‑1 (MQ=255)
                                                                                                                     aGGCGTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTCGACAGCCCTCACCACTTCGCGACAAAGGCGCCGTGTTGCCCCGCTTGCac                               >  2:94844/1‑149 (MQ=255)
                                                                                                                       gggcggcgcagtacctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCa                             <  2:286496‑M2/122‑1 (MQ=255)
                                                                                                                               atatCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCg                     <  1:1281238/150‑1 (MQ=255)
                                                                                                                                  tacctgtgggagcgggcGTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCgggg                  <  1:308284‑M2/133‑1 (MQ=255)
                                                                                                                                   cccTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGtgttgccccggcatgctgatgcaagcggggc                 <  2:1353443‑M1/150‑32 (MQ=255)
                                                                                                                                   acctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGc                 <  1:1001149‑M2/134‑1 (MQ=255)
                                                                                                                                    cctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGcc                                                 >  2:121934‑M2/16‑117 (MQ=255)
                                                                                                                                    cctgtgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGcc                                                 <  1:121934‑M2/102‑1 (MQ=255)
                                                                                                                                        tgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCg            <  1:141410‑M2/139‑1 (MQ=255)
                                                                                                                                        tgggagcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCg            >  1:1318580‑M2/12‑150 (MQ=255)
                                                                                                                                         tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                             <  2:364561‑M2/90‑1 (MQ=255)
                                                                                                                                         tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                             >  1:364561‑M2/11‑100 (MQ=255)
                                                                                                                                         tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCggg                   <  2:1338091‑M2/132‑1 (MQ=255)
                                                                                                                                         tttcgcgggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCggg                   >  1:1338091‑M2/11‑142 (MQ=255)
                                                                                                                                          ggAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGtgttgccctggcatgctgatgcaagcggggcaacacgg          <  1:970239‑M1/150‑39 (MQ=255)
                                                                                                                                          ggAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGgctgccccggcatgctgatgcaagcggggcaacacggc          <  2:355134‑M1/150‑39 (MQ=255)
                                                                                                                                                ggcATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttt    >  2:585951‑M2/4‑150 (MQ=255)
                                                                                                                                                  cATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttttt  <  1:1300942‑M2/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTT  >  NC_002947/1899174‑1899468

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