Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,455,095 C→G 6.8% G172A (GGA→GCA)  fbp ← fructose‑1,6‑bisphosphatase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,455,0950CG 6.8% 407.0 / 15.1 118G172A (GGA→GCA) fbpfructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base C (55/55);  new base G (5/3);  total (60/58)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 119 total aligned reads displayed.

GTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGA  >  NC_000913/4454950‑4455243
                                                                                                                                                 |                                                                                                                                                    
gTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg                                                                                                                                                   >  4:114292/1‑149 (MQ=255)
gTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCg                                                                                                                                                    >  1:709905/1‑148 (MQ=255)
gTACTTCTTCACCCAGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg                                                                                                                                                   >  4:887132/1‑149 (MQ=255)
    ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtggt                                                                                                                                               >  3:138955/1‑149 (MQ=255)
    ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtgg                                                                                                                                                <  2:577528/148‑1 (MQ=255)
     tcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtggt                                                                                                                                               >  1:316356/1‑148 (MQ=255)
       ttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGgtgt                                                                                                                                             >  2:709398/1‑148 (MQ=255)
           ccccGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAAcc                                                                                                                                        <  4:642583/149‑1 (MQ=255)
            cccGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggt                                                                                                                                                  <  1:799843/138‑1 (MQ=255)
            cccGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggt                                                                                                                                                  >  2:799843/1‑138 (MQ=255)
              cGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAg                                                                                                                                      >  1:597206/1‑148 (MQ=255)
              cGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGc                                                                                                                                     <  3:698075/149‑1 (MQ=255)
                ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCAt                                                                                                                                   >  3:263086/1‑149 (MQ=255)
                  cGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATgg                                                                                                                                 <  4:34103/149‑1 (MQ=255)
                   ggAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGACTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGt                                                                                                                                >  3:717962/1‑149 (MQ=255)
                    gAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGt                                                                                                                                >  1:257469/1‑148 (MQ=255)
                    gAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGACACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGt                                                                                                                                >  1:169917/1‑148 (MQ=255)
                    gAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCCGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTa                                                                                                                               <  4:536370/149‑1 (MQ=255)
                     aaaCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTa                                                                                                                               >  2:357352/1‑148 (MQ=255)
                       aCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTagag                                                                                                                            <  4:263086/149‑1 (MQ=255)
                            aTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGcc                                                                                                                       >  4:894520/1‑149 (MQ=255)
                            aTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGc                                                                                                                        <  2:853183/148‑1 (MQ=255)
                                    tCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATAcc                                                                                                                <  1:28197/148‑1 (MQ=255)
                                      cTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtg                                                                                                                                                 >  3:788823/1‑113 (MQ=255)
                                          gTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAc                                                                                                         >  4:699296/1‑149 (MQ=255)
                                          gTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAc                                                                                                         >  4:492669/1‑149 (MQ=255)
                                           ttGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAc                                                                                                         >  1:604904/1‑148 (MQ=255)
                                            tGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAcc                                                                                                        <  1:292090/148‑1 (MQ=255)
                                             gATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgtt                                                                                                                                               <  4:711780/108‑3 (MQ=255)
                                              aTGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgt                                                                                                                                                <  2:403625/106‑2 (MQ=255)
                                              aTGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgt                                                                                                                                                >  1:403625/1‑105 (MQ=255)
                                               tGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGcc                                                                                                    <  4:138955/149‑1 (MQ=255)
                                               tGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGc                                                                                                     <  2:597206/148‑1 (MQ=255)
                                                    tAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCaa                                                                                                >  1:905173/1‑148 (MQ=255)
                                                       gTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACGAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                            >  4:593443/1‑149 (MQ=255)
                                                        ttttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGACACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                            <  2:169917/148‑1 (MQ=255)
                                                         tttGTCTTTCTCCGGGAAGCGCATCCGTTCCTGGCAAAGGCAGAAACCGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCAGCCGGTGGTGTAAACCAGCGTGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgt                                                                                           <  2:251885/148‑1 (MQ=255)
                                                         tttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtt                                                                                          <  3:807312/149‑1 (MQ=255)
                                                         tttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgt                                                                                           <  2:257469/148‑1 (MQ=255)
                                                          ttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgttt                                                                                         <  3:114292/149‑1 (MQ=255)
                                                            gCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtttgt                                                                                       >  3:77574/1‑149 (MQ=255)
                                                            gCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtttgt                                                                                       >  4:821903/1‑149 (MQ=255)
                                                            gCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtttg                                                                                        <  1:905181/148‑1 (MQ=255)
                                                             ccTTTCTCCGGGAAGCGCATCCGTCCCTGGCACAGGCAGAAAACGCAGAGCGAAGGATCGTAAGTAAAGGCGTGAACCCCACATCCGGTGGTGTAAACCAGCTTGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtttgt                                                                                       <  2:490588/148‑1 (MQ=255)
                                                              cTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTa                                                                                     <  3:107590/149‑1 (MQ=255)
                                                              cTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTa                                                                                     >  4:709170/1‑149 (MQ=255)
                                                               tttCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTAc                                                                                    <  3:657986/149‑1 (MQ=255)
                                                                           gcgcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTggagg                                                                         >  1:829214/1‑148 (MQ=255)
                                                                            cgcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTggagga                                                                        >  2:792477/1‑148 (MQ=255)
                                                                             gcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGgaagaag                                                                      >  3:93174/1‑148 (MQ=255)
                                                                             gcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGaa                                                                       <  2:709905/148‑1 (MQ=255)
                                                                               aTCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtc                                                                    <  4:834910/149‑1 (MQ=255)
                                                                               aTCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAt                                                                     <  2:86246/148‑1 (MQ=255)
                                                                                 ccGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtct                                                                   <  2:430440/148‑1 (MQ=255)
                                                                                    ttCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATcttctt                                                               <  3:139169/149‑1 (MQ=255)
                                                                                     tCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATcttctt                                                               <  2:316356/148‑1 (MQ=255)
                                                                                      ccTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGCTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTcc                                                             <  4:512388/149‑1 (MQ=255)
                                                                                        tGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCg                                                            >  1:73285/1‑148 (MQ=255)
                                                                                          gCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGtt                                                          >  1:13131/1‑148 (MQ=255)
                                                                                            acaGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAc                                                        >  1:165938/1‑148 (MQ=255)
                                                                                                gCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCg                                                   >  3:691912/1‑149 (MQ=255)
                                                                                                 cAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCg                                                   <  1:709398/148‑1 (MQ=255)
                                                                                                  aGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAcc                                                       >  3:564674/1‑143 (MQ=255)
                                                                                                  aGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAcc                                                       <  4:564674/143‑1 (MQ=255)
                                                                                                  aGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTg                                                 >  4:89958/1‑149 (MQ=255)
                                                                                                   gAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTg                                                 >  2:880407/1‑148 (MQ=255)
                                                                                                     aaaCGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCa                                              <  3:894520/149‑1 (MQ=255)
                                                                                                      aaCGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCaa                                             >  4:584172/1‑149 (MQ=255)
                                                                                                      aaCGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCa                                              >  2:667065/1‑148 (MQ=255)
                                                                                                         gCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACa                                           <  2:786705/148‑1 (MQ=255)
                                                                                                         gCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAg                                          >  3:843805/1‑149 (MQ=255)
                                                                                                           cGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTcc                                                             >  4:388494/1‑128 (MQ=255)
                                                                                                           cGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTcc                                                             <  3:388494/128‑1 (MQ=255)
                                                                                                            gAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                            >  1:802249/1‑96 (MQ=255)
                                                                                                            gAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                            <  2:802249/96‑1 (MQ=255)
                                                                                                              gCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTa                                     >  4:104824/1‑149 (MQ=255)
                                                                                                                   ggATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtctt                                                                  <  1:798879/115‑1 (MQ=255)
                                                                                                                    gATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACgcgg                               >  4:767735/1‑149 (MQ=255)
                                                                                                                        gtaagtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGTAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGt                           >  3:452595/1‑149 (MQ=255)
                                                                                                                         taagtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTa                          >  4:663589/1‑149 (MQ=255)
                                                                                                                           agtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAg                         >  1:159546/1‑148 (MQ=255)
                                                                                                                           agtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAg                         <  2:13131/148‑1 (MQ=255)
                                                                                                                           agtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGa                        <  3:104921/149‑1 (MQ=255)
                                                                                                                            gtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGAt                       <  3:89958/149‑1 (MQ=255)
                                                                                                                             taaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGAt                       <  2:165938/148‑1 (MQ=255)
                                                                                                                                   cGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGaa                 >  2:708289/1‑148 (MQ=255)
                                                                                                                                    gTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGaaaa               >  4:546347/1‑149 (MQ=255)
                                                                                                                                    gTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGaaa                <  1:191397/148‑1 (MQ=255)
                                                                                                                                     tGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGga                    <  1:11352/143‑1 (MQ=255)
                                                                                                                                     tGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGga                    >  2:11352/1‑143 (MQ=255)
                                                                                                                                     tGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTAACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGGGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAAt              <  3:663589/149‑1 (MQ=255)
                                                                                                                                      gAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATg             <  3:767735/149‑1 (MQ=255)
                                                                                                                                       accacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGt                                                  >  4:23194/3‑111 (MQ=255)
                                                                                                                                       accacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGt                                                  <  3:23194/109‑1 (MQ=255)
                                                                                                                                       accacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATgg            >  4:429612/3‑149 (MQ=255)
                                                                                                                                        acacCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATgg            >  2:125767/1‑148 (MQ=255)
                                                                                                                                         cacCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGt           >  2:674306/1‑148 (MQ=255)
                                                                                                                                         cacCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGACGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTa          <  3:699296/149‑1 (MQ=255)
                                                                                                                                          acCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTAc         <  3:492669/149‑1 (MQ=255)
                                                                                                                                                 ccGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACaga  >  4:209199/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGA  >  NC_000913/4454950‑4455243

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: