Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,455,092 C→G 6.8% C173S (TGC→TCC)  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,0920CG 6.8% 408.1 / 17.3 118C173S (TGC→TCC) 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 = 4.35e-01

Only 100 of 118 total aligned reads displayed.

GTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTAC  >  NC_000913/4454950‑4455236
                                                                                                                                              |                                                                                                                                                
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)
         cACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTaaa                                                                                                                                   <  4:28232/149‑1 (MQ=255)
            cccGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggt                                                                                                                                           <  1:799843/138‑1 (MQ=255)
            cccGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggt                                                                                                                                           >  2:799843/1‑138 (MQ=255)
              cGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGc                                                                                                                              <  3:698075/149‑1 (MQ=255)
                ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCAt                                                                                                                            >  3:263086/1‑149 (MQ=255)
                  cGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATgg                                                                                                                          <  4:34103/149‑1 (MQ=255)
                  cGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATg                                                                                                                           <  1:547113/148‑1 (MQ=255)
                   ggAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGACTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGt                                                                                                                         >  3:717962/1‑149 (MQ=255)
                    gAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGt                                                                                                                         >  1:257469/1‑148 (MQ=255)
                    gAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGACACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGt                                                                                                                         >  1:169917/1‑148 (MQ=255)
                     aaaCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTa                                                                                                                        >  2:357352/1‑148 (MQ=255)
                            aTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGcc                                                                                                                >  4:894520/1‑149 (MQ=255)
                            aTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGc                                                                                                                 <  2:853183/148‑1 (MQ=255)
                             tGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGcc                                                                                                                >  1:786705/1‑148 (MQ=255)
                                    tCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATAcc                                                                                                         <  1:28197/148‑1 (MQ=255)
                                      cTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtg                                                                                                                                          >  3:788823/1‑113 (MQ=255)
                                       ttCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACg                                                                                                      >  2:191397/1‑148 (MQ=255)
                                          gTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAc                                                                                                  >  4:492669/1‑149 (MQ=255)
                                          gTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAc                                                                                                  >  4:699296/1‑149 (MQ=255)
                                           ttGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAc                                                                                                  >  1:604904/1‑148 (MQ=255)
                                            tGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAAcc                                                                                                 <  1:292090/148‑1 (MQ=255)
                                             gATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgtt                                                                                                                                        >  3:711780/1‑106 (MQ=255)
                                             gATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgtt                                                                                                                                        <  4:711780/108‑3 (MQ=255)
                                              aTGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgt                                                                                                                                         >  1:403625/1‑105 (MQ=255)
                                              aTGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGGATGCggtgt                                                                                                                                         <  2:403625/106‑2 (MQ=255)
                                               tGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGcc                                                                                             <  4:138955/149‑1 (MQ=255)
                                               tGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGc                                                                                              <  2:597206/148‑1 (MQ=255)
                                                 gAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGc                                                                                           >  4:657986/1‑149 (MQ=255)
                                                   gTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCaa                                                                                         >  3:512388/1‑149 (MQ=255)
                                                    tAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCaa                                                                                         >  1:905173/1‑148 (MQ=255)
                                                       gTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACGAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                     >  4:593443/1‑149 (MQ=255)
                                                        ttttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGACACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                     <  2:169917/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                                                                                >  4:821903/1‑149 (MQ=255)
                                                            gCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtttg                                                                                 <  1:905181/148‑1 (MQ=255)
                                                             ccTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtttgt                                                                                <  2:422639/148‑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)
                                                                             gcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGaa                                                                <  2:709905/148‑1 (MQ=255)
                                                                               aTCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtc                                                             <  4:834910/149‑1 (MQ=255)
                                                                               aTCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAt                                                              <  2:86246/148‑1 (MQ=255)
                                                                                tCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtct                                                            >  4:537957/1‑149 (MQ=255)
                                                                                 ccGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtct                                                            <  2:430440/148‑1 (MQ=255)
                                                                                  cGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtctt                                                           <  2:905173/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)
                                                                                           cacaGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAc                                                 <  4:114264/149‑1 (MQ=255)
                                                                                            acaGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAc                                                 >  1:165938/1‑148 (MQ=255)
                                                                                                 cAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCg                                            <  1:709398/148‑1 (MQ=255)
                                                                                                 cAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGt                                           <  3:147732/149‑1 (MQ=255)
                                                                                                  aGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAcc                                                <  4:564674/143‑1 (MQ=255)
                                                                                                  aGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTAcc                                                >  3:564674/1‑143 (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:111062/1‑148 (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)
                                                                                                            gAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtg                                                                                     >  1:802249/1‑96 (MQ=255)
                                                                                                              gCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTa                              >  4:104824/1‑149 (MQ=255)
                                                                                                                   ggATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAtctt                                                           >  2:798879/1‑115 (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)
                                                                                                                            gtaaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGAt                <  3:821903/149‑1 (MQ=255)
                                                                                                                             taaAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGAt                <  2:165938/148‑1 (MQ=255)
                                                                                                                                    gTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGaaaa        >  4:546347/1‑149 (MQ=255)
                                                                                                                                    gTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGaaa         <  1:191397/148‑1 (MQ=255)
                                                                                                                                     tGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGga             <  1:11352/143‑1 (MQ=255)
                                                                                                                                     tGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTAACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGGGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAAt       <  3:663589/149‑1 (MQ=255)
                                                                                                                                      gAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATg      <  3:767735/149‑1 (MQ=255)
                                                                                                                                       accacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGt                                           <  3:23194/109‑1 (MQ=255)
                                                                                                                                       accacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGt                                           >  4:23194/3‑111 (MQ=255)
                                                                                                                                       accacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATgg     >  4:429612/3‑149 (MQ=255)
                                                                                                                                        ccacCGGATGCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATgg     >  2:672728/2‑148 (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)
                                                                                                                                              |                                                                                                                                                
GTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTAC  >  NC_000913/4454950‑4455236

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: