Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,021,703 A→G 8.3% Q189R (CAG→CGG)  sxy → transcriptional coactivator for CRP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,021,7030AG 8.3% 493.7 / 25.5 172Q189R (CAG→CGG) sxytranscriptional coactivator for CRP
Reads supporting (aligned to +/- strand):  ref base A (72/84);  new base G (10/4);  total (82/90)
Fisher's exact test for biased strand distribution p-value = 9.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 176 total aligned reads displayed.

GGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGC  >  NC_000913/1021554‑1021852
                                                                                                                                                     |                                                                                                                                                     
ggATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCa                                                                                                                                                       <  2:2338749/150‑1 (MQ=255)
 gATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAg                                                                                                                                                      >  2:802515/1‑150 (MQ=255)
    gTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCTGGCGGCGTGCCACCGGGAg                                                                                                                                                   <  1:628730/150‑1 (MQ=255)
    gTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGACGCCAGGAg                                                                                                                                                   <  1:465987/150‑1 (MQ=255)
     tACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCTGGCGGCGTGCCACCGGGAGc                                                                                                                                                  <  1:2145883/150‑1 (MQ=255)
      aCGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCt                                                                                                                                                 >  2:2024554/1‑150 (MQ=255)
        gggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTg                                                                                                                                               <  1:53024/150‑1 (MQ=255)
        gggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGGTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTg                                                                                                                                               <  1:2393980/150‑1 (MQ=255)
               aCGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATg                                                                                                                                        <  1:968833/150‑1 (MQ=255)
               aCGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATg                                                                                                                                        <  2:1287636/150‑1 (MQ=255)
                cgcataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                                       <  2:314635/147‑1 (MQ=255)
                cGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                                       <  1:2110701/150‑1 (MQ=255)
                cGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                                       <  2:1074923/150‑1 (MQ=255)
                 gTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATggg                                                                                                                                      >  1:670437/1‑150 (MQ=255)
                  tataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGc                                                                                                                                     <  1:105437/150‑1 (MQ=255)
                  tataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGc                                                                                                                                     <  2:346181/150‑1 (MQ=255)
                      cTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGAc                                                                                                                                 >  2:286369/1‑150 (MQ=255)
                         ggggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActct                                                                                                                              >  2:2074138/1‑150 (MQ=255)
                         ggggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActct                                                                                                                              >  2:405935/1‑150 (MQ=255)
                             cAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTa                                                                                                                          <  1:360274/150‑1 (MQ=255)
                              aaaaaTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAc                                                                                                                         >  1:730817/1‑150 (MQ=255)
                              aaaaaTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAc                                                                                                                         <  1:1280143/150‑1 (MQ=255)
                                  aTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCg                                                                                                                     <  1:1144512/150‑1 (MQ=255)
                                  aTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCg                                                                                                                     <  2:878384/150‑1 (MQ=255)
                                     tgtTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGaaa                                                                                                                  >  1:757555/1‑150 (MQ=255)
                                     tgtTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAAGCCggat                                                                                                                  <  1:2360732/150‑4 (MQ=255)
                                      gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                                 <  1:759339/150‑1 (MQ=255)
                                      gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                                 <  1:184623/150‑1 (MQ=255)
                                         ggttgCGACTGCGGCAGCAAATCAGTCTGGTGACAGAAAACATTCAGTTTATGCTTGAAGGTGCCAGTTTCGGCAGTCATGAAGCTGCGCTCCCGGTGGCACGCCGGCAGGAGCTTGCAGAATGGCCTGACTCTCTTACGCCGAAACAgg                                                                                                              >  2:2126532/1‑150 (MQ=255)
                                         ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCATAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                                              <  1:1582486/150‑1 (MQ=255)
                                         ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                                              <  2:260987/150‑1 (MQ=255)
                                          gttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGa                                                                                                             >  2:935104/1‑150 (MQ=255)
                                            tgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGt                                                                                                           >  1:1593144/1‑150 (MQ=255)
                                                aCTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTcc                                                                                                       <  1:182184/150‑1 (MQ=255)
                                                       cagcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGaa                                                                                                >  2:1012746/1‑150 (MQ=255)
                                                          caAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACtt                                                                                             >  2:2194666/1‑150 (MQ=255)
                                                           aaaaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                                            >  2:78724/1‑150 (MQ=255)
                                                           aaaaCAGTCTGGTGACAGAAAACATTCTGTTTATGCTTGAAGGGGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                                            <  2:451110/150‑1 (MQ=255)
                                                            aaaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGa                                                                                           >  2:2419218/1‑150 (MQ=255)
                                                             aaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                                          <  2:1131313/150‑1 (MQ=255)
                                                             aaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                                          <  2:108414/150‑1 (MQ=255)
                                                              aCAGTCTGGTGACAGAAAAGATTCTGGTTATGCCTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGt                                                                                         <  2:1585527/150‑1 (MQ=255)
                                                                     ggTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTgg                                                                                  >  1:2078179/1‑150 (MQ=255)
                                                                       ggACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCg                                                                                <  1:2468422/149‑1 (MQ=255)
                                                                            gAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCa                                                                           >  2:1501321/1‑150 (MQ=255)
                                                                              aaaGATTCTGTTTATGCTTGAAGGTGCCATTATGGGCATTCATGAAGCTGCGCTCCGGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAga                                                                         <  2:1512307/150‑1 (MQ=255)
                                                                                 gATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAgacga                                                                      <  1:84083/150‑1 (MQ=255)
                                                                                  aTTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAgacgac                                                                     <  2:778561/150‑1 (MQ=255)
                                                                                   ttCTGTTTATGCTTGATGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCGTGCAGAATGGGCTGACGCTCTTACGCCGAAACAGGAGGTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGccc                                                                    <  1:2212208/150‑1 (MQ=255)
                                                                                   ttCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGAcc                                                                    <  1:2220867/150‑1 (MQ=255)
                                                                                       gTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGATCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtc                                                                <  1:2369852/150‑1 (MQ=255)
                                                                                       gTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtc                                                                >  2:1951327/1‑150 (MQ=255)
                                                                                       gTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtc                                                                >  2:1764694/1‑150 (MQ=255)
                                                                                         ttATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGGGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGATTAATCTGGCGTTGCAGACGACCAAtctc                                                              <  2:1730172/150‑1 (MQ=255)
                                                                                          tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCa                                                             <  1:1959469/150‑1 (MQ=255)
                                                                                          tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATATGGCGTTGCAGACGACCAATCTCa                                                             >  2:2056093/1‑150 (MQ=255)
                                                                                            tGCTTGAAGGTGCCATTATCGCCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTCCGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGGTGCAGACGACCAATCTCAgg                                                           <  2:1656386/150‑1 (MQ=255)
                                                                                             gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                          >  2:678668/1‑150 (MQ=255)
                                                                                             gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                          <  1:1064290/150‑1 (MQ=255)
                                                                                                tGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCATAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGt                                                       <  2:2246279/150‑1 (MQ=255)
                                                                                                 gAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTa                                                      <  2:462144/150‑1 (MQ=255)
                                                                                                   aGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAAt                                                    <  1:802515/150‑1 (MQ=255)
                                                                                                   aGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAAt                                                    >  2:1939354/1‑150 (MQ=255)
                                                                                                     gTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGc                                                  <  2:34219/150‑1 (MQ=255)
                                                                                                      tGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGca                                                 >  1:1290327/1‑150 (MQ=255)
                                                                                                        ccATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaat                                               >  2:1128805/1‑150 (MQ=255)
                                                                                                         cATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaatc                                              <  2:773176/150‑1 (MQ=255)
                                                                                                         cATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaatc                                              >  1:528132/1‑150 (MQ=255)
                                                                                                               cGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCaataat                                        <  1:241713/150‑1 (MQ=255)
                                                                                                                ggCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATc                                       >  1:315412/1‑150 (MQ=255)
                                                                                                                  cATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCa                                     <  1:2024554/150‑1 (MQ=255)
                                                                                                                  cATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCa                                     >  1:1128218/1‑150 (MQ=255)
                                                                                                                    ttCATGAAGCGGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAAc                                   <  1:188010/150‑1 (MQ=255)
                                                                                                                      cATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                                 >  2:1809325/1‑150 (MQ=255)
                                                                                                                      cATGAAGCTGCGCTCCCGGTGGCACGCCGACAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                                 >  2:1057561/1‑150 (MQ=255)
                                                                                                                      cATGAAGCTGCGCTCCCGGTGGCACCCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                                 >  1:1813101/1‑150 (MQ=255)
                                                                                                                       aTGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACttg                                >  1:881460/1‑150 (MQ=255)
                                                                                                                           ctgcaagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTtgttgt                            >  2:1466920/7‑150 (MQ=255)
                                                                                                                           aGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTtgttgt                            >  2:1898621/1‑150 (MQ=255)
                                                                                                                           aGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTtgttgt                            <  2:5480/150‑1 (MQ=255)
                                                                                                                             gcaagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg                         >  1:180884/5‑150 (MQ=255)
                                                                                                                             cTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAAAGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCTTTCCAGACGACCAATCTCAGGCAGTAATCCAATCAATAATCCAACTTGATGTAc                          >  2:2273730/1‑150 (MQ=255)
                                                                                                                              tGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg                         >  1:412797/1‑150 (MQ=255)
                                                                                                                               gcgcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGa                        >  2:2108773/1‑150 (MQ=255)
                                                                                                                                cgcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAc                       <  2:2071936/150‑1 (MQ=255)
                                                                                                                                agcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAc                       >  2:1275089/2‑150 (MQ=255)
                                                                                                                                 gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc                      <  1:1764694/150‑1 (MQ=255)
                                                                                                                                  cTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCa                     >  1:1420427/1‑150 (MQ=255)
                                                                                                                                  cTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCa                     >  2:629839/1‑150 (MQ=255)
                                                                                                                                    cccGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCaga                   >  1:1862961/1‑150 (MQ=255)
                                                                                                                                            cACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCttt           >  2:1410308/1‑150 (MQ=255)
                                                                                                                                            cACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCttt           <  1:2194666/150‑1 (MQ=255)
                                                                                                                                             aCGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTg          <  1:377090/150‑1 (MQ=255)
                                                                                                                                              cgccgcCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGc         <  2:341975/150‑1 (MQ=255)
                                                                                                                                              cgccgcCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTGCTTTGc         <  1:1484776/150‑1 (MQ=255)
                                                                                                                                               gccgccAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCt        <  2:2209430/150‑1 (MQ=255)
                                                                                                                                                ccgccAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTg       <  1:167987/150‑1 (MQ=255)
                                                                                                                                                ccgccAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTg       <  2:821929/150‑1 (MQ=255)
                                                                                                                                                  gccAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTc     <  2:663455/150‑1 (MQ=255)
                                                                                                                                                     aGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGc  <  1:2074138/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGC  >  NC_000913/1021554‑1021852

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