Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,326,750 C→G 7.9% F214L (TTC→TTG atoB → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,7500CG 7.9% 396.8 / 21.0 114F214L (TTC→TTGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (57/48);  new base G (3/6);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 114 total aligned reads displayed.

AAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATC  >  NC_000913/2326607‑2326896
                                                                                                                                               |                                                                                                                                                  
aaaGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTc                                                                                                                                                >  1:896233/1‑148 (MQ=255)
    agTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGa                                                                                                                                            >  2:733305/1‑148 (MQ=255)
    agTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGAc                                                                                                                                           <  4:510986/149‑1 (MQ=255)
      tACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACg                                                                                                                                          >  1:453849/1‑148 (MQ=255)
      tACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGa                                                                                                                                         >  3:639318/1‑149 (MQ=255)
        cGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaag                                                                                                                                       <  3:196440/149‑2 (MQ=255)
        cGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaag                                                                                                                                       <  4:488011/149‑2 (MQ=255)
        cGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaa                                                                                                                                        <  2:405642/148‑1 (MQ=255)
        cGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaa                                                                                                                                        <  1:834917/148‑1 (MQ=255)
         ggAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAAtt                                                                                                                                      <  3:353096/149‑1 (MQ=255)
             ttACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTccc                                                                                                                                   >  2:595113/1‑148 (MQ=255)
                 ccGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAg                                                                                                                              >  4:725890/1‑149 (MQ=255)
                  cGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAg                                                                                                                              >  1:664702/1‑148 (MQ=255)
                      aaaTGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAATACGAATTCCCGAAAGCGAAt                                                                                                                         <  4:581644/149‑1 (MQ=255)
                       aaTGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAAtt                                                                                                                        >  3:21554/1‑149 (MQ=255)
                          gCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCaa                                                                                                                     <  3:479992/149‑1 (MQ=255)
                          gCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCa                                                                                                                      >  2:537373/1‑148 (MQ=255)
                          gCAGGATGAACTGGCGCTACATTCACAGCGCAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCa                                                                                                                      <  1:586894/148‑1 (MQ=255)
                           cAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAAc                                                                                                                    <  3:3841/149‑1 (MQ=255)
                              gATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGc                                                                                                                 >  4:782662/1‑149 (MQ=255)
                               aTGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGc                                                                                                                 >  1:901166/1‑148 (MQ=255)
                                  aaCTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGaa                                                                                                             >  4:573032/1‑149 (MQ=255)
                                    cTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAg                                                                                                            <  2:805646/148‑1 (MQ=255)
                                    cTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAg                                                                                                            >  2:379716/1‑148 (MQ=255)
                                    cTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGc                                                                                                           <  4:599621/149‑1 (MQ=255)
                                     tGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAg                                                                                                                                             >  2:644293/1‑114 (MQ=255)
                                     tGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAg                                                                                                                                             <  1:644293/114‑1 (MQ=255)
                                      ggCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGt                                                                                                         >  4:187406/1‑149 (MQ=255)
                                      ggCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGt                                                                                                         <  4:835087/149‑1 (MQ=255)
                                        cgcTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGtt                                                                                                        <  2:775421/148‑1 (MQ=255)
                                         gcTACATTCACAGCGAAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAg                                                                                                      >  4:278236/1‑149 (MQ=255)
                                          cTACAATCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAg                                                                                                      >  2:646800/1‑148 (MQ=255)
                                           tACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAgg                                                                                                     <  1:646800/148‑1 (MQ=255)
                                           tACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGt                                                                                                    <  4:472785/149‑1 (MQ=255)
                                               ttCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCAt                                                                                                >  3:235334/1‑149 (MQ=255)
                                               ttCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCAt                                                                                                >  4:837360/1‑149 (MQ=255)
                                                tCACAGCGTAAAGCGGCAGCCGCAATTGAGGTCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCAt                                                                                                <  1:162861/148‑1 (MQ=255)
                                                 cacaGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAg                                                                                                                                             >  1:151104/1‑102 (MQ=255)
                                                 cacaGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAg                                                                                                                                             <  2:151104/102‑1 (MQ=255)
                                                   caGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCAt                                                                                                <  1:804725/145‑1 (MQ=255)
                                                   caGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCAt                                                                                                >  2:804725/1‑145 (MQ=255)
                                                          aaGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGc                                                                                      >  1:672665/1‑148 (MQ=255)
                                                           aGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCt                                                                                    <  3:212844/149‑1 (MQ=255)
                                                             cGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTc                                                                                  <  3:292970/149‑1 (MQ=255)
                                                              ggCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCGTTGCGCCCGGCCTTc                                                                                  <  1:595113/148‑1 (MQ=255)
                                                               gCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTgcgc                                                                                           <  4:342123/138‑1 (MQ=255)
                                                               gCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTgcgc                                                                                           >  3:342123/1‑138 (MQ=255)
                                                                   ccGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                             >  2:857899/1‑148 (MQ=255)
                                                                    cGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaag                                                                                                                                       >  3:877299/1‑88 (MQ=255)
                                                                    cGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaaa                                                                            <  2:664702/148‑1 (MQ=255)
                                                                     gCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGc                                                                          >  3:343326/1‑149 (MQ=255)
                                                                      cAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGc                                                                          >  2:763221/1‑148 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                             >  3:348980/1‑142 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                             <  4:348980/142‑1 (MQ=255)
                                                                           gAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAAc                                                                    >  3:45294/1‑149 (MQ=255)
                                                                            aGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAAc                                                                    <  1:733305/148‑1 (MQ=255)
                                                                              tCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAAAGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGt                                                                 <  3:187406/149‑1 (MQ=255)
                                                                               ccGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGt                                                                 >  1:706297/1‑148 (MQ=255)
                                                                                cGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCa                                                               <  3:837360/149‑1 (MQ=255)
                                                                                cGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCa                                                               <  3:192321/149‑1 (MQ=255)
                                                                                 ggTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                             >  2:658340/1‑134 (MQ=255)
                                                                                      ttttACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCt                                                          <  2:453849/148‑1 (MQ=255)
                                                                                        ttACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTggg                                                       <  3:797201/149‑1 (MQ=255)
                                                                                             gCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                   >  1:451143/1‑148 (MQ=255)
                                                                                             gCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                   <  2:581573/148‑1 (MQ=255)
                                                                                                 aaaTCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAATAGTCACCGCTGGGAACGCGTc                                               >  2:710170/1‑148 (MQ=255)
                                                                                                 aaaTCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTc                                               >  2:709790/1‑148 (MQ=255)
                                                                                                     cGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTa                                          <  3:278236/149‑1 (MQ=255)
                                                                                                       tCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAAtt                                                                                                                                      >  2:274924/1‑55 (MQ=255)
                                                                                                       tCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAAtt                                                                                                                                      <  1:274924/55‑1 (MQ=255)
                                                                                                                 tgttgtCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAg                                                                        >  3:443782/1‑107 (MQ=255)
                                                                                                                  gttgtCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGtgc                             >  3:53167/1‑149 (MQ=255)
                                                                                                                     gtCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGtgctgc                          >  3:698520/1‑149 (MQ=255)
                                                                                                                      tCACTTGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGtgctgc                          >  2:646904/1‑148 (MQ=255)
                                                                                                                       cACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGcc                         >  1:559565/1‑148 (MQ=255)
                                                                                                                         cTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGc                       >  1:531060/1‑148 (MQ=255)
                                                                                                                             aaagaaaACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTgg                  <  4:289225/149‑1 (MQ=255)
                                                                                                                              aagaaaACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTgg                  >  1:735242/1‑148 (MQ=255)
                                                                                                                               agaaaACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTg                >  3:202371/1‑149 (MQ=255)
                                                                                                                                gaaaACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTAGTATTAACGACGGTGCTGCCGCTCTGGTGa               >  4:274086/1‑149 (MQ=255)
                                                                                                                                  aaaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAtt             <  4:639318/149‑1 (MQ=255)
                                                                                                                                   aaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAtt             >  1:852793/1‑148 (MQ=255)
                                                                                                                                    aCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCTAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTa            <  1:451202/148‑1 (MQ=255)
                                                                                                                                    aCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTa            >  1:309583/1‑148 (MQ=255)
                                                                                                                                     ccTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                   >  2:594858/1‑108 (MQ=255)
                                                                                                                                     ccTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                   <  1:594858/108‑1 (MQ=255)
                                                                                                                                     ccTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATg          <  4:21554/149‑1 (MQ=255)
                                                                                                                                      cTTCGTCTTCCGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGTTTATgg         >  3:492165/1‑149 (MQ=255)
                                                                                                                                      attCGTCTTGACTGAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGa               <  4:350956/142‑1 (MQ=255)
                                                                                                                                      attCGTCTTGACTGAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGa               >  3:350956/2‑143 (MQ=255)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGct                      >  1:727564/1‑135 (MQ=255)
                                                                                                                                        tCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGgaa       <  4:698520/149‑1 (MQ=255)
                                                                                                                                        tCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGga        <  2:896233/148‑1 (MQ=255)
                                                                                                                                          gTCTTGACTGAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGgaaga     >  3:75912/1‑149 (MQ=255)
                                                                                                                                          gTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGgaaga     <  4:492165/149‑1 (MQ=255)
                                                                                                                                          gTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGgaaga     >  4:774691/1‑149 (MQ=255)
                                                                                                                                            cTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGgaagaa    <  1:635796/148‑1 (MQ=255)
                                                                                                                                             ttCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAAt   >  1:753589/1‑148 (MQ=255)
                                                                                                                                             ttCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATc  >  4:147972/1‑149 (MQ=255)
                                                                                                                                             ttCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATc  <  3:826102/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
AAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATC  >  NC_000913/2326607‑2326896

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