Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 329,703 A→T 5.9% K81* (AAG→TAG)  betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913329,7030AT 5.9% 244.6 / 3.7 80K81* (AAG→TAG) betTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base T (2/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

AACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGTATATGCAGCCGCCGGAAGGCGCGGGAC  >  NC_000913/329555‑329850
                                                                                                                                                    |                                                                                                                                                   
aaCGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGa                                                                                                                                                     <  6:25042/149‑1 (MQ=255)
    aTCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCt                                                                                                                                                 >  7:353599/1‑149 (MQ=255)
     tCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTc                                                                                                                                                >  1:324854/1‑149 (MQ=255)
     tCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTc                                                                                                                                                >  7:59958/1‑149 (MQ=255)
        tGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCggg                                                                                                                                             >  6:11302/1‑149 (MQ=255)
        tGTTTCGCGACTTCTCGGCCCTGTGGAATGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCggg                                                                                                                                             >  1:274201/1‑149 (MQ=255)
            tCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCa                                                                                                                                          <  2:324854/148‑1 (MQ=255)
              gcgACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGaa                                                                                                                                       >  8:82007/1‑149 (MQ=255)
                gACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAAc                                                                                                                                      >  6:126226/1‑148 (MQ=255)
                     ctcGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATcc                                                                                                                                 >  2:68077/1‑148 (MQ=255)
                      tcGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCa                                                                                                                                >  2:341927/1‑148 (MQ=255)
                                 gATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAAtt                                                                                                                      >  1:194502/1‑147 (MQ=255)
                                     ggCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGcc                                                                                                                 <  5:11302/148‑1 (MQ=255)
                                      gCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATc                                                                                                                                  >  1:343433/1‑130 (MQ=255)
                                      gCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATc                                                                                                                                  <  2:343433/130‑1 (MQ=255)
                                       ccGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCctgc                                                                                                              >  4:124901/1‑149 (MQ=255)
                                       ccGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCctgc                                                                                                              <  8:59958/149‑1 (MQ=255)
                                        cGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAACctgct                                                                                                             <  5:60758/149‑1 (MQ=255)
                                               tGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTACTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGaa                                                                                                                                       >  5:270313/1‑116 (MQ=255)
                                               tGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTACTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGaa                                                                                                                                       <  6:270313/116‑1 (MQ=255)
                                                 gACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTggg                                                                                                     <  1:341927/148‑1 (MQ=255)
                                                      ggTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTAAATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGa                                                                                               <  5:313938/149‑1 (MQ=255)
                                                       gTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGAt                                                                                              <  5:126226/149‑1 (MQ=255)
                                                        tttCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATg                                                                                             <  3:124901/149‑1 (MQ=255)
                                                          tCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGc                                                                                            <  7:82007/148‑1 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAACCCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGt                                                                                         >  2:384209/1‑149 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGc                                                                                                                  <  8:223184/124‑1 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGc                                                                                                                  >  7:223184/1‑124 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGAt                                                                                              >  4:347921/1‑144 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGAt                                                                                              <  3:347921/144‑1 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGt                                                                                         >  7:220872/1‑149 (MQ=255)
                                                            tAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGt                                                                                         >  2:247557/1‑149 (MQ=255)
                                                             aaaaCCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGtt                                                                                        >  4:274722/1‑149 (MQ=255)
                                                              aaaCCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGttt                                                                                       >  3:354646/1‑149 (MQ=255)
                                                                aCCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTtgc                                                                                     >  3:80042/1‑149 (MQ=255)
                                                                  cTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTtgctg                                                                                   >  2:149719/1‑149 (MQ=255)
                                                                     cGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGcc                                                                                 <  2:194502/148‑1 (MQ=255)
                                                                      ggTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCgg                                                                               <  1:149719/149‑1 (MQ=255)
                                                                      ggTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCgg                                                                               >  7:322086/1‑149 (MQ=255)
                                                                        ttGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGTAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGa                                                                             >  4:212493/1‑149 (MQ=255)
                                                                        ttGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGGGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGa                                                                             <  1:197795/149‑1 (MQ=255)
                                                                          ggTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGAt                                                                            <  8:353599/148‑1 (MQ=255)
                                                                           gTACTATCTGCTGGCGGCAACGCTCTATATTGTCTGTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCg                                                                          >  1:227269/1‑149 (MQ=255)
                                                                             aCTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAACCCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGATTGCTGCCGGGATCgg                                                                         >  3:269922/1‑148 (MQ=255)
                                                                               tATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATc                                                                                                                                  <  7:18973/89‑1 (MQ=255)
                                                                               tATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATc                                                                                                                                  <  7:378229/89‑1 (MQ=255)
                                                                               tATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATc                                                                                                                                  >  8:378229/1‑89 (MQ=255)
                                                                                   tgcgtgcggcAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGGGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATc                                                                     <  2:244477/141‑1 (MQ=255)
                                                                                      tgGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCt                                                                >  7:296069/1‑148 (MQ=255)
                                                                                       ggcggcAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGa                                                              >  6:182937/1‑149 (MQ=255)
                                                                                        gcttgaaCGCTCTATTTTGTCTTTGTGGTCTTTATCGCTTGTTCGCGTTTTGGTTCGGTGACGCTCGGGCCTGAACATTCGAAACCGGATTTCTGCCTGCTGAGTTGGGCGGCGATGCTGGTTGCTTCCGGGATCGGTATCGACCTGAt                                                             <  4:269922/144‑1 (MQ=255)
                                                                                        gcggcAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGAt                                                             <  1:68077/149‑1 (MQ=255)
                                                                                         cggcAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATg                                                            <  5:182937/149‑1 (MQ=255)
                                                                                            cAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGTAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGttc                                                         <  3:212493/149‑1 (MQ=255)
                                                                                            cAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGttc                                                         <  4:80042/149‑1 (MQ=255)
                                                                                            cAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGttc                                                         <  1:134834/149‑1 (MQ=255)
                                                                                            cAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTGGTTCGCGTTGTGGTTCGGTGTAGCTCGGGCCAGAACATTCCTAACCGGATTTCTGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGttc                                                         <  3:279953/149‑1 (MQ=255)
                                                                                                gCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCgg                                                                               <  8:288343/123‑1 (MQ=255)
                                                                                                gCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCgg                                                                               >  7:288343/1‑123 (MQ=255)
                                                                                                gCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTtcttct                                                     >  4:315850/1‑149 (MQ=255)
                                                                                                gCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTCTGCTGCCGGGATCGGTATCGACCTGATGTtcttct                                                     <  8:220872/149‑1 (MQ=255)
                                                                                                 ctctATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTtctc                                                    >  6:350110/1‑149 (MQ=255)
                                                                                                      tatTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGTAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAg                                               >  2:357822/1‑149 (MQ=255)
                                                                                                      tatTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCGAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAg                                               >  8:208791/1‑149 (MQ=255)
                                                                                                      tatTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAATCCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAg                                               >  3:56764/1‑149 (MQ=255)
                                                                                                      tatTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTa                                                <  8:296069/148‑1 (MQ=255)
                                                                                                          gTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAACCCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGa                                           >  6:246055/1‑149 (MQ=255)
                                                                                                               tgtgGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAAGAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCGGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCg                                       <  2:274201/148‑1 (MQ=255)
                                                                                                                     cTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACg                                 >  5:95308/1‑148 (MQ=255)
                                                                                                                        tATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGt                             <  2:227269/149‑1 (MQ=255)
                                                                                                                          tCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGCCCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGta                            >  3:320003/1‑148 (MQ=255)
                                                                                                                          tCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGta                            <  1:247557/148‑1 (MQ=255)
                                                                                                                          tCGCTTGTTCGCGTTTTGGTTCGGTCAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTAGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGta                            >  4:189857/1‑148 (MQ=255)
                                                                                                                                    tcgTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGTATATGCAgccgc                 >  5:256119/2‑149 (MQ=255)
                                                                                                                                    gcgTTTTGGTTCGGTGTAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGTATATGCAgccgc                 <  1:357822/149‑1 (MQ=255)
                                                                                                                                       ttttGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGttc                                                         >  7:308025/1‑106 (MQ=255)
                                                                                                                                       ttttGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGttc                                                         <  8:308025/106‑1 (MQ=255)
                                                                                                                                               cGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAAc                                         <  3:342400/114‑1 (MQ=255)
                                                                                                                                               cGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAAc                                         >  4:342400/1‑114 (MQ=255)
                                                                                                                                                ggTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGGTCTTCTCCGTAGCCGAACCGGTAACGCAGTATATGCAGCCGCCGGGAGGCGCgg     <  3:274722/149‑1 (MQ=255)
                                                                                                                                                   gAAGCTCGGGCCAGAACAATCCACACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTAGCGACCTGATGTTCTTCTCCGTAGCCGAACCGGGAACGCAGTATATGCAGCCGCCGGACGGCGCGGGAc  >  8:5399/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTGTATCGCTTGTTCGCGTTTTGGTTCGGTGAAGCTCGGGCCAGAACAATCCAAACCGGAATTCAGCCTGCTGAGTTGGGCGGCGATGCTGTTTGCTGCCGGGATCGGTATCGACCTGATGTTCTTCTCCGTAGCCGAACCGGTAACGCAGTATATGCAGCCGCCGGAAGGCGCGGGAC  >  NC_000913/329555‑329850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: