Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,021,688 C→T 8.1% 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.1% 516.3 / 24.5 174P184L (CCG→CTG) sxytranscriptional coactivator for CRP
Reads supporting (aligned to +/- strand):  ref base C (71/85);  new base T (9/5);  total (81/93)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

Only 100 of 176 total aligned reads displayed.

AGGTGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAG  >  NC_000913/1021541‑1021836
                                                                                                                                                   |                                                                                                                                                    
aGGTGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCTGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCgg                                                                                                                                                    <  2:1283694/150‑1 (MQ=255)
   tGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTgg                                                                                                                                                 <  1:1259975/150‑1 (MQ=255)
    ggggATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGc                                                                                                                                                <  1:2222234/150‑1 (MQ=255)
      ggATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAc                                                                                                                                              <  2:1175822/150‑1 (MQ=255)
      ggATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAc                                                                                                                                              <  2:1560570/150‑1 (MQ=255)
         ttAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgcc                                                                                                                                           >  1:746001/1‑150 (MQ=255)
         ttAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGGGACAGAAAAGATGCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgcc                                                                                                                                           <  2:1732522/150‑1 (MQ=255)
           aaGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgccgc                                                                                                                                         >  2:167987/1‑150 (MQ=255)
           aaGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgccgc                                                                                                                                         >  2:240204/1‑150 (MQ=255)
           aaGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAcgccgc                                                                                                                                         <  1:925407/150‑1 (MQ=255)
            aGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACgccgcc                                                                                                                                        <  2:245012/150‑1 (MQ=255)
                 gTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGACGCCAGGAg                                                                                                                                   <  1:465987/150‑1 (MQ=255)
                     gggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTg                                                                                                                               <  1:53024/150‑1 (MQ=255)
                     gggCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGGTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTg                                                                                                                               <  1:2393980/150‑1 (MQ=255)
                        cGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCTGGCGGCGTGCCACCGGGAGcgcagct                                                                                                                            <  1:488150/150‑7 (MQ=255)
                        cGTTACGTATACTTGGGGCAAAAATGCGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAg                                                                                                                            <  2:2168867/150‑1 (MQ=255)
                            aCGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATg                                                                                                                        <  1:968833/150‑1 (MQ=255)
                             cgcataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                       <  2:314635/147‑1 (MQ=255)
                             cGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATgg                                                                                                                       <  2:1074923/150‑1 (MQ=255)
                              gTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATggg                                                                                                                      >  1:670437/1‑150 (MQ=255)
                               tataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGc                                                                                                                     <  2:346181/150‑1 (MQ=255)
                               tataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGc                                                                                                                     <  1:622059/150‑1 (MQ=255)
                               tataCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGc                                                                                                                     <  1:105437/150‑1 (MQ=255)
                                  aCTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGa                                                                                                                  >  1:34219/1‑150 (MQ=255)
                                   cTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGAc                                                                                                                 >  1:77654/1‑150 (MQ=255)
                                     tGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGAATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGATTGCAGAATGGGCTGActa                                                                                                               >  2:780349/1‑149 (MQ=255)
                                      ggggCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGActct                                                                                                              >  2:405935/1‑150 (MQ=255)
                                               aTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCg                                                                                                     <  1:1144512/150‑1 (MQ=255)
                                               aTGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCg                                                                                                     <  2:878384/150‑1 (MQ=255)
                                                  tgtTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTAAGCCggat                                                                                                  <  1:2360732/150‑4 (MQ=255)
                                                   gttggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAAc                                                                                                 <  1:184623/150‑1 (MQ=255)
                                                      ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCATAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                              <  1:1582486/150‑1 (MQ=255)
                                                      ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                              <  2:2032313/150‑1 (MQ=255)
                                                      ggttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAgg                                                                                              >  1:245884/1‑150 (MQ=255)
                                                       gttgCGACTGCGGCATCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCTGAATTCATGAAGCTGCGCTCCCGTTGGCACGCCGCGATGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGa                                                                                             <  2:2087675/150‑1 (MQ=255)
                                                       gttgCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGa                                                                                             >  2:935104/1‑150 (MQ=255)
                                                             aCTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTcc                                                                                       <  1:182184/150‑1 (MQ=255)
                                                              cTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCTGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCt                                                                                      <  2:77654/150‑1 (MQ=255)
                                                                gCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGc                                                                                    <  1:2103637/150‑1 (MQ=255)
                                                                 cGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCg                                                                                   >  1:757735/1‑150 (MQ=255)
                                                                  ggCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCgg                                                                                  <  2:2256422/150‑1 (MQ=255)
                                                                  ggCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCgg                                                                                  >  2:470106/1‑150 (MQ=255)
                                                                   gcagcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGa                                                                                 >  2:1724343/1‑150 (MQ=255)
                                                                    cagcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGaa                                                                                >  2:1339650/1‑150 (MQ=255)
                                                                    cagcaAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGaa                                                                                >  2:1012746/1‑150 (MQ=255)
                                                                       caAAAGAGTCTGGTGACAGAAAAGATTCTGTTTATGGTTGAAGGTGGCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCAGGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACtt                                                                             <  2:1593144/150‑1 (MQ=255)
                                                                       caAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACtt                                                                             >  2:2194666/1‑150 (MQ=255)
                                                                        aaaaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                            >  2:78724/1‑150 (MQ=255)
                                                                        aaaaCAGTCTGGTGACAGAAAACATTCTGTTTATGCTTGAAGGGGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTg                                                                            <  2:451110/150‑1 (MQ=255)
                                                                          aaCAGTGTGGTTTCAGAAAAGATTGTGTTTATGGTTGAAGGTGCCTTTATGGGCATTGATGAAGGTGCGCTCCGGGTGGCACGCCGCCAGGTGGTTGCAGAATGGTCTGTCTCTCTTACGCCTTAACATGAGTTTCCTGCGGAACTTGAg                                                                          <  2:2300397/150‑1 (MQ=255)
                                                                          aaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                          <  2:108414/150‑1 (MQ=255)
                                                                          aaCAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAg                                                                          <  2:1131313/150‑1 (MQ=255)
                                                                                cTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCt                                                                    >  2:2404586/1‑150 (MQ=255)
                                                                                  ggTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTgg                                                                  >  2:2121097/1‑150 (MQ=255)
                                                                                  ggTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTgg                                                                  >  1:2078179/1‑150 (MQ=255)
                                                                                  ggTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTgg                                                                  <  1:1535033/150‑1 (MQ=255)
                                                                                    tGACAGAAAAGATTCTGTTTATGCTTGGAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCg                                                                <  1:1257921/150‑1 (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:309301/150‑1 (MQ=255)
                                                                                                ttCTGTTTATGCTTGATGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCGTGCAGAATGGGCTGACGCTCTTACGCCGAAACAGGAGGTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGccc                                                    <  1:2212208/150‑1 (MQ=255)
                                                                                                ttCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGAcc                                                    >  1:2167524/1‑150 (MQ=255)
                                                                                                    gTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGATCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtc                                                <  1:2369852/150‑1 (MQ=255)
                                                                                                    gTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCCGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAAtc                                                >  2:662459/1‑150 (MQ=255)
                                                                                                       tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCa                                             <  1:1638778/150‑1 (MQ=255)
                                                                                                       tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCa                                             <  1:1959469/150‑1 (MQ=255)
                                                                                                       tATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATATGGCGTTGCAGACGACCAATCTCa                                             >  2:2056093/1‑150 (MQ=255)
                                                                                                          gCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGc                                          <  1:1064290/150‑1 (MQ=255)
                                                                                                                aGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAAt                                    >  2:1939354/1‑150 (MQ=255)
                                                                                                                aGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAAt                                    <  1:802515/150‑1 (MQ=255)
                                                                                                                 ggTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATg                                   >  2:808958/1‑150 (MQ=255)
                                                                                                                  gTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGc                                  <  2:34219/150‑1 (MQ=255)
                                                                                                                   tGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGca                                 >  1:1290327/1‑150 (MQ=255)
                                                                                                                   tGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGAAGAATGGGCTGCCTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGca                                 >  2:2006086/1‑150 (MQ=255)
                                                                                                                     ccATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaat                               >  2:1128805/1‑150 (MQ=255)
                                                                                                                     ccATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaat                               >  2:273306/1‑150 (MQ=255)
                                                                                                                     ccATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGcaat                               >  1:952536/1‑150 (MQ=255)
                                                                                                                          aTCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCaata                          >  2:709580/1‑150 (MQ=255)
                                                                                                                             ggCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAACCTCAGGCAGTAATGCAATCAATAATc                       >  2:1885099/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)
                                                                                                                                   cATGAAGCTGCGCTCCCGGTGGCACGCCGACAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                 >  2:1057561/1‑150 (MQ=255)
                                                                                                                                   cATGAAGCTGCGCTCCCGGTGGCACCCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACtt                 >  1:1813101/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)
                                                                                                                                          gcaagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATGCAACTTGTTGTAc          >  1:2118715/5‑150 (MQ=255)
                                                                                                                                          gcaagcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAAT‑GGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg         >  1:180884/5‑150 (MQ=255)
                                                                                                                                           tGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACg         >  1:412797/1‑150 (MQ=255)
                                                                                                                                             cgcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACTCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACTACCAATGTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAc       <  2:1775417/150‑1 (MQ=255)
                                                                                                                                             cgcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAc       <  2:2071936/150‑1 (MQ=255)
                                                                                                                                             agcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAc       >  2:1275089/2‑150 (MQ=255)
                                                                                                                                              gcTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc      >  1:686561/1‑150 (MQ=255)
                                                                                                                                              gcTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGAcc      <  1:1764694/150‑1 (MQ=255)
                                                                                                                                               cTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCa     >  2:629839/1‑150 (MQ=255)
                                                                                                                                               cTCCTGGCGGCGTGCCACCGGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCa     >  1:1420427/1‑150 (MQ=255)
                                                                                                                                                 cccGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCaga   >  1:1862961/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGTGGGGATTAAGGATGTACGGGCGTTACGTATACTTGGGGCAAAAATGTGTTGGTTGCGACTGCGGCAGCAAAACAGTCTGGTGACAGAAAAGATTCTGTTTATGCTTGAAGGTGCCATTATCGGCATTCATGAAGCTGCGCTCCCGGTGGCACGCCGCCAGGAGCTTGCAGAATGGGCTGACTCTCTTACGCCGAAACAGGAGTTTCCTGCGGAACTTGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAG  >  NC_000913/1021541‑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: ‑