Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,071,095 A→T 6.5% A50A (GCT→GCA fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,071,0950AT 6.5% 319.1 / 11.2 95A50A (GCT→GCAfbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base A (41/46);  new base T (3/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-02

GTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCTTAGACA  >  NC_000913/3070947‑3071243
                                                                                                                                                    |                                                                                                                                                    
gTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTagca                                                                                                                                                      >  4:402883/1‑149 (MQ=255)
      ccATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGtt                                                                                                                                                >  1:277115/1‑149 (MQ=255)
      ccATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGtt                                                                                                                                                >  3:122703/1‑149 (MQ=255)
        aTCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTc                                                                                                                                              >  3:666243/1‑149 (MQ=255)
          cTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCa                                                                                                                                            <  2:141311/149‑1 (MQ=255)
          cTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCa                                                                                                                                            <  1:446376/149‑1 (MQ=255)
          cTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCa                                                                                                                                            <  1:729227/149‑1 (MQ=255)
          cTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCa                                                                                                                                            <  3:622709/149‑1 (MQ=255)
           tGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAg                                                                                                                                           <  2:490582/149‑1 (MQ=255)
                aaCGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACgg                                                                                                                                      >  4:556431/1‑149 (MQ=255)
                 aCGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACgac                                                                                                                                     >  4:687889/1‑147 (MQ=255)
                 aCGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGc                                                                                                                                     >  3:489283/1‑149 (MQ=255)
                   gTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGt                                                                                                                                   >  1:814930/1‑149 (MQ=255)
                    tGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGtt                                                                                                                                  <  3:27171/149‑1 (MQ=255)
                     gATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTg                                                                                                                                 >  4:598516/1‑149 (MQ=255)
                       tGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGAt                                                                                                                               <  2:407276/149‑1 (MQ=255)
                       tGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGAt                                                                                                                               <  3:634117/149‑1 (MQ=255)
                       tGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGAt                                                                                                                               <  1:5460/149‑1 (MQ=255)
                            aCCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGt                                                                                                                          <  3:687889/149‑1 (MQ=255)
                            aCCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGt                                                                                                                          >  3:781899/1‑149 (MQ=255)
                              cAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCa                                                                                                                        >  1:547788/1‑149 (MQ=255)
                               agagATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAg                                                                                                                       <  1:277781/149‑1 (MQ=255)
                               agagATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAg                                                                                                                       <  4:122703/149‑1 (MQ=255)
                               agagATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAg                                                                                                                       <  2:277115/149‑1 (MQ=255)
                                  gaTCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTAc                                                                                                                    >  3:277636/1‑149 (MQ=255)
                                          ccAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCTGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAg                                                                                                            >  2:571712/1‑149 (MQ=255)
                                          ccAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCTGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAg                                                                                                            >  4:779020/1‑149 (MQ=255)
                                          ccAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCTGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAg                                                                                                            >  2:164225/1‑149 (MQ=255)
                                              gATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCGGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGt                                                                                                           >  1:46865/1‑146 (MQ=255)
                                              gATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCGGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGt                                                                                                           <  2:46865/146‑1 (MQ=255)
                                              gATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTa                                                                                                        >  4:631383/1‑149 (MQ=255)
                                              gATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTa                                                                                                        <  2:598257/149‑1 (MQ=255)
                                                    agcaCCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCAGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTAct                                                                                                      <  2:671192/145‑1 (MQ=255)
                                                    agcaCCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCAGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTAct                                                                                                      >  1:671192/1‑145 (MQ=255)
                                                         ccTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGt                                                                                             >  1:793676/1‑149 (MQ=255)
                                                            gCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCg                                                                                          <  2:793676/149‑1 (MQ=255)
                                                              ggAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGaa                                                                                        <  2:404810/149‑1 (MQ=255)
                                                              ggAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCAGGCAGTGCGaa                                                                                        <  2:730394/149‑1 (MQ=255)
                                                              ggAACGTCAGAGTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGaa                                                                                        <  2:304980/149‑1 (MQ=255)
                                                               gAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAg                                                                                       >  2:1257/1‑149 (MQ=255)
                                                                   gTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAgttgt                                                                                   >  1:10764/1‑149 (MQ=255)
                                                                   gTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAgttgt                                                                                   <  3:402883/149‑1 (MQ=255)
                                                                   gTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAgttgt                                                                                   >  3:458170/1‑149 (MQ=255)
                                                                     cAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGttt                                                                                 >  3:350861/1‑149 (MQ=255)
                                                                      aGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGtttt                                                                                <  4:210636/149‑1 (MQ=255)
                                                                      aGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGtttt                                                                                >  2:526822/1‑149 (MQ=255)
                                                                      aGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGtttt                                                                                >  4:764641/1‑149 (MQ=255)
                                                                         tttCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCtt                                                                             <  3:608875/149‑1 (MQ=255)
                                                                         tttCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCtt                                                                             <  2:10764/149‑1 (MQ=255)
                                                                          ttCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCttt                                                                            <  2:150021/149‑1 (MQ=255)
                                                                           tCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCtttt                                                                           <  3:764641/149‑1 (MQ=255)
                                                                            cACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTg                                                                          <  3:598516/149‑1 (MQ=255)
                                                                             aCGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGa                                                                                                                            >  3:250255/1‑98 (MQ=255)
                                                                             aCGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGa                                                                                                                            <  4:250255/98‑1 (MQ=255)
                                                                             aCGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGc                                                                         >  1:399928/1‑149 (MQ=255)
                                                                             aCGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGc                                                                         >  2:746891/1‑149 (MQ=255)
                                                                               gCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTa                                                                       >  4:543943/1‑149 (MQ=255)
                                                                               gCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTa                                                                       >  4:707293/1‑149 (MQ=255)
                                                                                 cTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTAcc                                                                     <  1:1257/149‑1 (MQ=255)
                                                                                  tttACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCt                                                                              <  2:718218/139‑1 (MQ=255)
                                                                                  tttACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCt                                                                              >  1:718218/1‑139 (MQ=255)
                                                                                  tttACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCt                                                                    <  1:6111/149‑1 (MQ=255)
                                                                                       cAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGaaa                                                               >  4:362581/1‑149 (MQ=255)
                                                                                        aGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGaaaa                                                              <  4:489283/149‑1 (MQ=255)
                                                                                             tAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTc                                                         <  4:665282/149‑1 (MQ=255)
                                                                                              aaaGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCt                                                        <  2:814930/149‑1 (MQ=255)
                                                                                               aaGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTg                                                       >  2:392731/1‑149 (MQ=255)
                                                                                                 ggAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTa                                                     >  3:712582/1‑149 (MQ=255)
                                                                                                  gAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTACCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTAc                                                    <  4:535204/149‑1 (MQ=255)
                                                                                                       accaccGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTActgc                                                                                                    <  1:537768/96‑1 (MQ=255)
                                                                                                       accaccGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTActgc                                                                                                    >  2:537768/1‑96 (MQ=255)
                                                                                                       accaccGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGtcat                                               >  1:627439/1‑149 (MQ=255)
                                                                                                        ccaccGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGtcatc                                              <  3:631383/149‑1 (MQ=255)
                                                                                                         caccGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGtcatca                                             <  2:168810/149‑1 (MQ=255)
                                                                                                          accGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCAc                                            >  4:179776/1‑149 (MQ=255)
                                                                                                           ccGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCAcc                                           >  4:310246/1‑149 (MQ=255)
                                                                                                            cGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCa                                          <  2:547788/149‑1 (MQ=255)
                                                                                                             gTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAg                                         >  1:113978/1‑149 (MQ=255)
                                                                                                              ttGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGt                                        <  4:712582/149‑1 (MQ=255)
                                                                                                                 gagaACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGAt                                     >  1:797846/1‑149 (MQ=255)
                                                                                                                    aaCTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTAc                                  <  4:666243/149‑1 (MQ=255)
                                                                                                                        gAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCa                              <  4:667396/149‑1 (MQ=255)
                                                                                                                               aaCCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTa                       >  3:361050/1‑149 (MQ=255)
                                                                                                                                aCCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTAc                      <  3:362581/149‑1 (MQ=255)
                                                                                                                                  cGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGa                    <  4:277636/149‑1 (MQ=255)
                                                                                                                                     cgcTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAAt                 <  2:399928/149‑1 (MQ=255)
                                                                                                                                        tttAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAAt                 <  1:112565/146‑1 (MQ=255)
                                                                                                                                        tttAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAAt                 >  2:112565/1‑146 (MQ=255)
                                                                                                                                        tttAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCaa              >  2:488660/1‑149 (MQ=255)
                                                                                                                                         ttAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCaaa             >  2:644835/1‑149 (MQ=255)
                                                                                                                                           aaCTTTAGCTGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCaaaaa           <  3:779020/149‑1 (MQ=255)
                                                                                                                                           aaCTTTAGCTGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCaaaaa           <  1:571712/149‑1 (MQ=255)
                                                                                                                                           aaCTTTAGCTGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCaaaaa           <  1:164225/149‑1 (MQ=255)
                                                                                                                                           aaCTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCaaaaa           <  4:781899/149‑1 (MQ=255)
                                                                                                                                             cTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATc         >  4:76689/1‑149 (MQ=255)
                                                                                                                                              tttAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCt        >  4:785754/1‑149 (MQ=255)
                                                                                                                                                  gcagcGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCTTAGa    <  3:707293/149‑1 (MQ=255)
                                                                                                                                                  gcagcGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCTTAGa    <  1:392731/149‑1 (MQ=255)
                                                                                                                                                  gcagcGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTATCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCTTAGa    <  3:543943/149‑1 (MQ=255)
                                                                                                                                                    agcGGTTTCCTGTACGGCGTTGATGGAGTCAGTACCAACGAAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTTCCTGGAAAACTTTCTGTACGTCATTACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCTTAGACa  >  3:446326/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGTTTCCAGTACGGCGTTGATGGAGTCAGTACCGACGCAGTTTACTGCTGGCAGTGCGAAGTTGTTTTCTTTTGCTACCTGGAAAACTTTCTGTACGTCATCACCAGTGATTACGCCAGGTTTTACGAAATCAAAAATCTTAGACA  >  NC_000913/3070947‑3071243

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: