Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,021,688 C→T 8.6% P184L (CCG→CTG)  sxy → transcriptional coactivator for CRP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,021,6880CT 8.6% 281.0 / 10.7 93P184L (CCG→CTG) sxytranscriptional coactivator for CRP
Reads supporting (aligned to +/- strand):  ref base C (45/40);  new base T (4/4);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GGTGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAG  >  NC_000913/1021542‑1021836
                                                                                                                                                    |                                                                                                                                                    
ggTGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGt                                                                                                                                                   <  1:595835/150‑1 (MQ=255)
  tGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTgg                                                                                                                                                 <  2:347798/150‑1 (MQ=255)
   ggggATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGc                                                                                                                                                <  1:1197176/150‑1 (MQ=255)
   ggggATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCACGGTGGc                                                                                                                                                >  1:187211/1‑150 (MQ=255)
       aTTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCAcgc                                                                                                                                            >  1:1193404/1‑150 (MQ=255)
        ttAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATCTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCAcgcc                                                                                                                                           <  1:772891/150‑1 (MQ=255)
          aaGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCAcgccgc                                                                                                                                         >  1:92806/1‑150 (MQ=255)
          aaGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCAcgccgc                                                                                                                                         <  1:506345/150‑1 (MQ=255)
            ggATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAGGTTGGTTTTTTGGTTTGAGGTGGCCATTAGGGGCTTACATAAGCTTC‑‑CGTTCCGGTGGCAAgccgcccc                                                                                                                                      >  1:563269/1‑148 (MQ=255)
             gATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAg                                                                                                                                      <  2:1038000/150‑1 (MQ=255)
              aTGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAgg                                                                                                                                     <  2:1108721/150‑1 (MQ=255)
                  ccGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCTGGCGGCGTGCCACCGGGAGcg                                                                                                                                 <  1:574426/149‑2 (MQ=255)
                  aCGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCTGGTGGCGTGCCACCGGGAGcg                                                                                                                                 <  1:267691/150‑2 (MQ=255)
                  aCGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCt                                                                                                                                 >  2:632840/1‑150 (MQ=255)
                   cGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGGTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCTGGCGGCGTGCCACCGGGAGcgc                                                                                                                                <  1:321145/150‑3 (MQ=255)
                     ggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGc                                                                                                                              >  1:110924/1‑150 (MQ=255)
                     ggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGc                                                                                                                              >  1:323073/1‑150 (MQ=255)
                     ggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGc                                                                                                                              >  2:1168275/1‑150 (MQ=255)
                      gCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCa                                                                                                                             <  2:1464939/150‑1 (MQ=255)
                       cGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCAAGCTCCTGGCGGCGTGCCACCGGGAGcgcag                                                                                                                              <  1:1478466/150‑5 (MQ=38)
                        gTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGa                                                                                                                           >  2:800860/1‑150 (MQ=255)
                           aCGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATg                                                                                                                        >  2:469475/1‑150 (MQ=255)
                                  cTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGAc                                                                                                                 >  1:1535582/1‑150 (MQ=255)
                                      gggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActctc                                                                                                             <  1:634638/150‑1 (MQ=255)
                                       ggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActctct                                                                                                            <  1:685804/150‑1 (MQ=255)
                                       ggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActctct                                                                                                            >  2:1379053/1‑150 (MQ=255)
                                        gCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCtt                                                                                                           <  2:1301709/150‑1 (MQ=255)
                                         cAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTa                                                                                                          >  1:463074/1‑150 (MQ=255)
                                         cAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTa                                                                                                          >  2:399335/1‑150 (MQ=255)
                                          aaaaaTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAc                                                                                                         >  1:314554/1‑150 (MQ=255)
                                             aaTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGcc                                                                                                      <  2:1044648/150‑1 (MQ=255)
                                             aaTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGcc                                                                                                      <  2:522502/150‑1 (MQ=255)
                                                gtgtTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGaa                                                                                                   >  1:924017/1‑150 (MQ=255)
                                                  gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                 >  1:784218/1‑150 (MQ=255)
                                                  gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                 >  2:1303813/1‑150 (MQ=255)
                                                  gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                 <  2:56893/150‑1 (MQ=255)
                                                  gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                 <  2:659659/150‑1 (MQ=255)
                                                     ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                              <  1:701206/150‑1 (MQ=255)
                                                     ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                              >  2:1402597/1‑150 (MQ=255)
                                                      gttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGa                                                                                             >  1:104055/1‑150 (MQ=255)
                                                      gttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCAATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGa                                                                                             <  2:910725/150‑1 (MQ=255)
                                                       ttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAg                                                                                            <  2:1296473/150‑1 (MQ=255)
                                                       ttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAg                                                                                            <  2:92806/150‑1 (MQ=255)
                                                       ttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAg                                                                                            <  1:1446106/150‑1 (MQ=255)
                                                              tGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTg                                                                                     <  1:469475/150‑1 (MQ=255)
                                                              tGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTg                                                                                     <  1:1379053/150‑1 (MQ=255)
                                                               gCGGCAGCAAAACAGTCTGGTGACAGAAAAGATGCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGc                                                                                    <  2:1167409/150‑1 (MQ=255)
                                                                 ggCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCgg                                                                                  <  1:1168275/150‑1 (MQ=255)
                                                                     gcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCa                                                                                                                             <  1:796159/103‑1 (MQ=255)
                                                                     gcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCa                                                                                                                             >  2:796159/1‑103 (MQ=255)
                                                                     gcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACt                                                                              >  2:456661/1‑150 (MQ=255)
                                                                       aaaaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                            >  1:580556/1‑150 (MQ=255)
                                                                       aaaaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                            <  2:158268/150‑1 (MQ=255)
                                                                        aaaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGa                                                                           >  2:1228655/1‑150 (MQ=255)
                                                                                 ggTGACAGAAAAGATTCTGTTTATGCTTGTAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTgg                                                                  >  2:576251/1‑150 (MQ=255)
                                                                                      cAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTg                                                             >  1:311974/1‑150 (MQ=255)
                                                                                       aGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGc                                                            >  2:981655/1‑150 (MQ=255)
                                                                                       aGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTACCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTATCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGc                                                            >  1:150709/1‑150 (MQ=255)
                                                                                        gAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCa                                                           >  2:237160/1‑150 (MQ=255)
                                                                                                tCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCa                                                   <  2:563269/150‑1 (MQ=255)
                                                                                                 cTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCaa                                                  >  1:1066479/1‑150 (MQ=255)
                                                                                                    tttATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCCAtct                                               >  1:714608/1‑150 (MQ=255)
                                                                                                     ttATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtctc                                              <  2:62378/150‑1 (MQ=255)
                                                                                                      tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCa                                             >  1:235551/1‑150 (MQ=255)
                                                                                                      tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCa                                             >  2:1449076/1‑150 (MQ=255)
                                                                                                       aTGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAg                                            <  1:1402597/150‑1 (MQ=255)
                                                                                                       aTGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAg                                            >  2:857968/1‑150 (MQ=255)
                                                                                                        tGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAgg                                           <  2:784218/150‑1 (MQ=255)
                                                                                                         gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                          <  1:800860/150‑1 (MQ=255)
                                                                                                         gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                          <  2:509725/150‑1 (MQ=255)
                                                                                                          cTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCa                                         >  1:987640/1‑150 (MQ=255)
                                                                                                                 gTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGc                                  >  1:1398095/1‑150 (MQ=255)
                                                                                                                 gTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGc                                  >  2:1499762/1‑150 (MQ=255)
                                                                                                                    ccATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaat                               <  1:889933/150‑1 (MQ=255)
                                                                                                                          tCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAGGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCaataa                         >  2:1528656/1‑150 (MQ=255)
                                                                                                                          tCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCaataa                         >  2:532364/1‑150 (MQ=255)
                                                                                                                           cGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCaataat                        <  1:199510/150‑1 (MQ=255)
                                                                                                                            ggCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATc                       >  1:880821/1‑150 (MQ=255)
                                                                                                                            ggCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATc                       >  2:254051/1‑150 (MQ=255)
                                                                                                                                  cATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACCCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                 >  1:364631/1‑150 (MQ=255)
                                                                                                                                    tGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGCAGCTTACAGAATGGGCTGACTCTCTTAATCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGAAGACCAATCTCAGCCAGTAATCCAATCAATAATCCAACttgt               >  2:1074571/1‑150 (MQ=255)
                                                                                                                                      aaGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAGGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACttgttg             <  1:1528656/150‑1 (MQ=255)
                                                                                                                                       aGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTtgttgt            >  2:566006/1‑150 (MQ=255)
                                                                                                                                          tGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg         <  2:104055/150‑1 (MQ=255)
                                                                                                                                          tGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg         <  2:924017/150‑1 (MQ=255)
                                                                                                                                            caagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTAGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg         >  2:1414814/4‑150 (MQ=255)
                                                                                                                                             aagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGa        >  2:213634/3‑150 (MQ=255)
                                                                                                                                              agcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAc       >  2:347667/2‑150 (MQ=255)
                                                                                                                                               gcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGACCTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc      >  2:747235/1‑150 (MQ=255)
                                                                                                                                               gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc      <  1:632840/150‑1 (MQ=255)
                                                                                                                                               gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc      <  1:399335/150‑1 (MQ=255)
                                                                                                                                               gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc      >  2:1016731/1‑150 (MQ=255)
                                                                                                                                                 tCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAAGTTGTTGTACGACCag    <  2:187211/150‑1 (MQ=255)
                                                                                                                                                  cccGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCaga   <  2:1535582/150‑1 (MQ=255)
                                                                                                                                                   ccGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCagag  >  2:1090865/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGC‑‑GCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAG  >  NC_000913/1021542‑1021836

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: