Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,326,754 C→G 7.8% Q216E (CAA→GAA)  atoB → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,7540CG 7.8% 404.8 / 21.5 116Q216E (CAA→GAA) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (58/49);  new base G (3/6);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 118 total aligned reads displayed.

AAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGG  >  NC_000913/2326607‑2326901
                                                                                                                                                   |                                                                                                                                                   
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)
                          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                                                                                                                                                  <  1:644293/114‑1 (MQ=255)
                                     tGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAg                                                                                                                                                  >  2:644293/1‑114 (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:212844/1‑149 (MQ=255)
                                           tACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGt                                                                                                         <  4:472785/149‑1 (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                                                                                                     >  2:804725/1‑145 (MQ=255)
                                                           aGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCt                                                                                         <  3:212844/149‑1 (MQ=255)
                                                             cGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTc                                                                                       <  3:292970/149‑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)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                                  >  3:348980/1‑142 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATaa                                                                                  <  4:348980/142‑1 (MQ=255)
                                                                         ttGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAgg                                                                            <  1:763221/148‑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                                                                                  <  1:658340/134‑1 (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                                                                                                                                           <  1:274924/55‑1 (MQ=255)
                                                                                                         ccGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTaa                                           <  3:725890/149‑1 (MQ=255)
                                                                                                                 tgttgtCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAg                                                                             <  4:443782/107‑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)
                                                                                                                         cTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGc                            >  2:340788/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)
                                                                                                                                 aaaaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGa                    >  2:635796/1‑148 (MQ=255)
                                                                                                                                 aaaaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAt                   >  3:40850/1‑149 (MQ=255)
                                                                                                                                  aaaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAtt                  <  4:639318/149‑1 (MQ=255)
                                                                                                                                   aaCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGAtt                  >  1:852793/1‑148 (MQ=255)
                                                                                                                                    aCCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTa                 >  1:309583/1‑148 (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)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGct                           <  2:727564/135‑1 (MQ=255)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGct                           >  1:727564/1‑135 (MQ=255)
                                                                                                                                       ttCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGga             >  4:603495/1‑149 (MQ=255)
                                                                                                                                        tCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGga             <  2:896233/148‑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       <  3:826102/149‑1 (MQ=255)
                                                                                                                                                aGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTg     <  2:672665/148‑1 (MQ=255)
                                                                                                                                                  tCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGCAGAATCTgcgg  >  3:585466/1‑149 (MQ=255)
                                                                                                                                                  tCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTgcgg  >  4:127160/1‑149 (MQ=255)
                                                                                                                                                   cAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTgcgg  >  1:243035/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACAGCCGAAATCGTCCCGGTAAATGTTGTCACTCGAAAGAAAACCTTCGTCTTCAGTCAAGACGAATTCCCGAAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGG  >  NC_000913/2326607‑2326901

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: