Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,326,752 G→C 7.8% S215T (AGT→ACT)  atoB → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,7520GC 7.8% 400.0 / 21.7 115S215T (AGT→ACT) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (57/49);  new base C (3/6);  total (60/55)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

Only 100 of 115 total aligned reads displayed.

AAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTG  >  NC_000913/2326607‑2326898
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                                          <  1:834917/148‑1 (MQ=255)
        cGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaa                                                                                                                                          <  2:405642/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)
                            aGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACg                                                                                                                     >  3:289225/1‑149 (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:379716/1‑148 (MQ=255)
                                    cTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAg                                                                                                              <  2:805646/148‑1 (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)
                                           tACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGt                                                                                                      >  4:212844/1‑149 (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)
                                                   caGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCAt                                                                                                  <  1:804725/145‑1 (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)
                                                                   ccGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                               >  2:857899/1‑148 (MQ=255)
                                                                    cGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTGACTGAAGACGaag                                                                                                                                         >  3:877299/1‑88 (MQ=255)
                                                                    cGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaaa                                                                              <  2:901166/148‑1 (MQ=255)
                                                                     gCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGc                                                                            >  3:343326/1‑149 (MQ=255)
                                                                      cAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGc                                                                            >  2:763221/1‑148 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                               <  4:348980/142‑1 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                               >  3:348980/1‑142 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAgg                                                                         <  1:763221/148‑1 (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:192321/149‑1 (MQ=255)
                                                                                cGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCa                                                                 <  3:837360/149‑1 (MQ=255)
                                                                                 ggTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                               >  2:658340/1‑134 (MQ=255)
                                                                                 ggTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                               <  1:658340/134‑1 (MQ=255)
                                                                                      ttttACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCt                                                            <  2:453849/148‑1 (MQ=255)
                                                                                             gCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                     <  2:581573/148‑1 (MQ=255)
                                                                                             gCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                     >  1:451143/1‑148 (MQ=255)
                                                                                                 aaaTCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAATAGTCACCGCTGGGAACGCGTc                                                 >  2:710170/1‑148 (MQ=255)
                                                                                                 aaaTCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTc                                                 >  2:709790/1‑148 (MQ=255)
                                                                                                     cGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTa                                            <  3:278236/149‑1 (MQ=255)
                                                                                                       tCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAAtt                                                                                                                                        <  1:274924/55‑1 (MQ=255)
                                                                                                       tCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAAtt                                                                                                                                        >  2:274924/1‑55 (MQ=255)
                                                                                                         ccGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTaa                                        <  3:725890/149‑1 (MQ=255)
                                                                                                                 tgttgtCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAg                                                                          >  3:443782/1‑107 (MQ=255)
                                                                                                                 tgttgtCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAg                                                                          <  4:443782/107‑1 (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                         >  2:340788/1‑148 (MQ=255)
                                                                                                                             aaagaaaACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTgg                    <  4:289225/149‑1 (MQ=255)
                                                                                                                                gaaaACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTAGTATTAACGACGGTGCTGCCGCTCTGGTGa                 >  4:274086/1‑149 (MQ=255)
                                                                                                                                 aaaaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGa                 >  2:635796/1‑148 (MQ=255)
                                                                                                                                 aaaaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAt                >  3:40850/1‑149 (MQ=255)
                                                                                                                                   aaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAtt               >  1:852793/1‑148 (MQ=255)
                                                                                                                                    aCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCTAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTa              <  1:451202/148‑1 (MQ=255)
                                                                                                                                    aCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTa              >  1:309583/1‑148 (MQ=255)
                                                                                                                                     ccTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAAcg                                                     <  1:594858/108‑1 (MQ=255)
                                                                                                                                      cTTCGTCTTCCGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGTTTATgg           >  3:492165/1‑149 (MQ=255)
                                                                                                                                      attCGTCTTGACTGAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGa                 >  3:350956/2‑143 (MQ=255)
                                                                                                                                      attCGTCTTGACTGAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGa                 <  4:350956/142‑1 (MQ=255)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGct                        >  1:727564/1‑135 (MQ=255)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGct                        <  2:727564/135‑1 (MQ=255)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGga          >  4:603495/1‑149 (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)
                                                                                                                                           tCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGgaaga       <  1:537373/148‑1 (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)
                                                                                                                                                aGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTg  <  2:672665/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTG  >  NC_000913/2326607‑2326898

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