Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,021,696 C→T 5.2% R187C (CGC→TGC)  sxy → transcriptional coactivator for CRP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,021,6960CT 5.2% 308.5 / 3.9 96R187C (CGC→TGC) sxytranscriptional coactivator for CRP
Reads supporting (aligned to +/- strand):  ref base C (48/43);  new base T (3/2);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 102 total aligned reads displayed.

GGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGC  >  NC_000913/1021547‑1021845
                                                                                                                                                     |                                                                                                                                                     
agATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTCACATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAc                                                                                                                                                       >  1:1346551/2‑150 (MQ=255)
   ttAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgcc                                                                                                                                                    >  1:477410/1‑150 (MQ=255)
   ttAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgcc                                                                                                                                                    >  2:1321973/1‑150 (MQ=255)
   ttAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgcc                                                                                                                                                    >  2:1188431/1‑150 (MQ=255)
       ggATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCa                                                                                                                                                >  1:113609/1‑150 (MQ=255)
        gATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAg                                                                                                                                               <  1:667154/150‑1 (MQ=255)
        gATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAg                                                                                                                                               >  1:286779/1‑150 (MQ=255)
        gATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAg                                                                                                                                               <  1:391237/150‑1 (MQ=255)
          tGTACGGGCGTTACGTATCCTTGGGGCAAATAGGTGTTGGTTGCGTCTGCGGCAGCAAAAGAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGa                                                                                                                                             <  2:924674/150‑1 (MQ=255)
           gTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCTGGCGGCGTGCCACCGGGAg                                                                                                                                            <  2:151014/150‑1 (MQ=255)
             aCGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCt                                                                                                                                          <  1:295078/150‑1 (MQ=255)
                ggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGc                                                                                                                                       >  1:1161540/1‑150 (MQ=255)
                ggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGc                                                                                                                                       >  1:446900/1‑150 (MQ=255)
                ggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGc                                                                                                                                       >  1:1155640/1‑150 (MQ=255)
                 gCGTGACGTATACTTGGGGCAAAAATGGGTTGGGTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCTGGCGGCGTGCCACCGGGAGcgcagc                                                                                                                                      <  1:320185/150‑6 (MQ=255)
                  cGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAg                                                                                                                                     >  1:821747/1‑150 (MQ=255)
                     tACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAAt                                                                                                                                  >  2:526262/1‑150 (MQ=255)
                      aCGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATg                                                                                                                                 <  1:526262/150‑1 (MQ=255)
                      aCGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATg                                                                                                                                 >  1:1479354/1‑150 (MQ=255)
                       cGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                                >  2:496693/1‑150 (MQ=255)
                         tataCTTGGGGCAAAAATGTGTTGGTTGCGAGTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATTCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTTGCACTCCGCCAGGAGCTTGCAGAATGGGc                                                                                                                              <  1:278076/150‑1 (MQ=255)
                          ataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCt                                                                                                                             <  1:662675/150‑1 (MQ=255)
                           taCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTg                                                                                                                            >  2:824563/1‑150 (MQ=255)
                             cTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGAc                                                                                                                          >  1:393584/1‑150 (MQ=255)
                             cTTGGGGCAAAAATGGGTTGGTTGCGACTGCGGCAGCAAAACAGTCTCGGGACAGAAAAGATTCTGTTTATGCGTGAAGGTCCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGAc                                                                                                                          <  2:1170113/150‑1 (MQ=255)
                               tGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                                >  1:795951/1‑142 (MQ=255)
                               tGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                                <  2:795951/142‑1 (MQ=255)
                                ggggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActct                                                                                                                       <  2:203955/150‑1 (MQ=255)
                                 gggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActctc                                                                                                                      <  2:1321687/150‑1 (MQ=255)
                                  ggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACTGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActctct                                                                                                                     <  2:642380/150‑1 (MQ=255)
                                  ggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActctct                                                                                                                     >  2:1159999/1‑150 (MQ=255)
                                     aaaaaTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAc                                                                                                                  <  2:483523/150‑1 (MQ=255)
                                            tgtTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCGAGGAGCTTGCAGAATGGGCTGACTGTCTTTCGCCGaaa                                                                                                           <  2:1125922/150‑1 (MQ=255)
                                                ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                                       <  2:62797/150‑1 (MQ=255)
                                                  ttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAg                                                                                                     <  1:1318392/150‑1 (MQ=255)
                                                      gACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTc                                                                                                 >  1:1095828/1‑150 (MQ=255)
                                                          gCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGc                                                                                             <  1:330422/150‑1 (MQ=255)
                                                           cGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCg                                                                                            >  2:305535/1‑150 (MQ=255)
                                                            ggCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCgg                                                                                           >  1:426795/1‑150 (MQ=255)
                                                             gcagcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGa                                                                                          >  1:449053/1‑150 (MQ=255)
                                                              cagcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGaa                                                                                         >  1:468670/1‑150 (MQ=255)
                                                               agcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAAc                                                                                        <  2:286779/150‑1 (MQ=255)
                                                                gcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACt                                                                                       >  1:1317034/1‑150 (MQ=255)
                                                                  aaaaCAGTCTGGTGACAGAAATGTTTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                                     <  2:1346551/150‑1 (MQ=255)
                                                                    aaCAGTCTGGTGACAGAAAATATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGTTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                                   <  2:1553328/150‑1 (MQ=255)
                                                                    aaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                                   <  2:1256925/150‑1 (MQ=255)
                                                                    aaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                                   <  2:514841/150‑1 (MQ=255)
                                                                     aCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGt                                                                                  <  1:276487/150‑1 (MQ=255)
                                                                              tGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCg                                                                         <  1:164085/150‑1 (MQ=255)
                                                                              tGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCg                                                                         >  2:1502344/1‑150 (MQ=255)
                                                                                 cAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTg                                                                      >  1:1621873/1‑150 (MQ=255)
                                                                                  aGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGc                                                                     <  2:1155640/150‑1 (MQ=255)
                                                                                  aGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGc                                                                     <  2:1161540/150‑1 (MQ=255)
                                                                                   gAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCa                                                                    <  2:763420/150‑1 (MQ=255)
                                                                                   gAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCa                                                                    <  1:305535/150‑1 (MQ=255)
                                                                                   gAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCa                                                                    >  2:748335/1‑150 (MQ=255)
                                                                                        gATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAgacga                                                               >  2:818546/1‑150 (MQ=255)
                                                                                         aTTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGAAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAgacgac                                                              >  2:264727/1‑150 (MQ=255)
                                                                                          ttatGTTTATGCTTGAACGTAACATTATCGTCAATAATGAATCTGCGCACCAGGTGTCACGCCTCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGAcc                                                             >  2:1469199/4‑150 (MQ=255)
                                                                                           tCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCa                                                            <  1:1321973/150‑1 (MQ=255)
                                                                                            cTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCaa                                                           <  2:821747/150‑1 (MQ=255)
                                                                                              gTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtc                                                         <  1:496693/150‑1 (MQ=255)
                                                                                                  aTGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAg                                                     <  2:1479354/150‑1 (MQ=255)
                                                                                                    gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCTAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                   <  2:1425245/150‑1 (MQ=255)
                                                                                                    gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                   <  1:70917/150‑1 (MQ=255)
                                                                                                    gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                   <  1:1502344/150‑1 (MQ=255)
                                                                                                    gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                   >  1:813828/1‑150 (MQ=255)
                                                                                                    gCTAGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                                   <  2:1241179/150‑1 (MQ=255)
                                                                                                     cTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCa                                                  >  1:1221484/1‑150 (MQ=255)
                                                                                                     cTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCa                                                  >  1:778160/1‑150 (MQ=255)
                                                                                                           ggTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATg                                            >  2:1619163/1‑150 (MQ=255)
                                                                                                            gTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGc                                           <  1:162509/150‑1 (MQ=255)
                                                                                                              gCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaa                                         >  1:1205340/1‑150 (MQ=255)
                                                                                                                cATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaatc                                       <  2:1621873/150‑1 (MQ=255)
                                                                                                                 aTTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaatca                                      <  2:728849/150‑1 (MQ=255)
                                                                                                                 aTTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaatca                                      >  1:1447409/1‑150 (MQ=255)
                                                                                                                       ggCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATc                                >  1:395144/1‑150 (MQ=255)
                                                                                                                             cATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                          >  1:475721/1‑150 (MQ=255)
                                                                                                                               tGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACttgt                        >  2:452787/1‑150 (MQ=255)
                                                                                                                                gAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtctc                                                       <  2:749475/118‑1 (MQ=255)
                                                                                                                                gAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtctc                                                       >  1:749475/1‑118 (MQ=255)
                                                                                                                                 aaGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACttgttg                      <  1:1159999/150‑1 (MQ=255)
                                                                                                                                  ctgcaagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTtgttgt                     >  2:94528/7‑150 (MQ=255)
                                                                                                                                    cTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTAc                   >  2:454266/1‑150 (MQ=255)
                                                                                                                                     tGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAACCAGGAGTTTCCTGCGGCACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCGAAAATAGAACTTGTTGTACg                  >  1:1222977/1‑150 (MQ=255)
                                                                                                                                     tGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg                  >  2:474882/1‑150 (MQ=255)
                                                                                                                                      gcgcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGa                 >  2:506607/1‑150 (MQ=255)
                                                                                                                                      aagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGa                 >  2:1112953/3‑150 (MQ=255)
                                                                                                                                        gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc               >  2:630436/1‑150 (MQ=255)
                                                                                                                                        gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc               <  2:446900/150‑1 (MQ=255)
                                                                                                                                        gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc               >  2:193642/1‑150 (MQ=255)
                                                                                                                                          tCCCGGTGGCATGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCag             >  1:935690/1‑150 (MQ=255)
                                                                                                                                            ccGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCagag           >  1:1026811/1‑150 (MQ=255)
                                                                                                                                            ccGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGATGTACGACCagag           >  1:538902/1‑150 (MQ=255)
                                                                                                                                             cGGTGGCACGCCTCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAgg          >  2:330304/1‑150 (MQ=255)
                                                                                                                                              ggTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGATGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGt         >  1:1631197/1‑150 (MQ=255)
                                                                                                                                                   cACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCttt    >  1:1347305/1‑150 (MQ=255)
                                                                                                                                                   cACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCttt    >  2:712115/1‑150 (MQ=255)
                                                                                                                                                     cgccgcCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGc  <  1:824563/150‑1 (MQ=255)
                                                                                                                                                     cgccgcCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGc  <  2:477410/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGC  >  NC_000913/1021547‑1021845

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: