Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 1,007,288 C→T F230F (TTC→TTT)  bglP → similar to PTS system, beta‑glucosides specific enzyme IIABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066711,007,2880CT99.9% 3282.7 / NA 829F230F (TTC→TTT) bglPsimilar to PTS system, beta‑glucosides specific enzyme IIABC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (381/447);  total (382/447)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 829 total aligned reads displayed.

GCACTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTCGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTA  >  AM406671/1007142‑1007432
                                                                                                                                                  |                                                                                                                                                
 caCTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCGCAGTACCCACAGCCGGCTATGGTTCAGCTTTAATCCCAATTATTGCCAATACCCCATTTGcc                                                                                                                                               >  2:5483014/1‑148 (MQ=255)
  aCTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAg                                                                                                                                              >  1:5585369/1‑149 (MQ=255)
   cTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGc                                                                                                                                             >  2:9881774/1‑149 (MQ=255)
     tGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCt                                                                                                                                           <  1:8347756/149‑1 (MQ=255)
      gTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCtt                                                                                                                                          >  1:4543575/1‑149 (MQ=255)
        ttATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCtttt                                                                                                                                        <  1:8198090/149‑1 (MQ=255)
        ttATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCtttt                                                                                                                                        <  1:9017996/149‑1 (MQ=255)
         tATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttttt                                                                                                                                       <  1:8873543/149‑1 (MQ=255)
         tATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttttt                                                                                                                                       <  2:3378990/149‑1 (MQ=255)
           tCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAg                                                                                                                                     >  2:4335727/1‑149 (MQ=255)
            ccAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGGTCATCTGTAATGCCAATTATTGCCATTAAAGCATTTGCAGCCTTTTTAGa                                                                                                                                    >  1:1669430/1‑149 (MQ=255)
             cAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                   <  2:1386968/149‑1 (MQ=255)
             cAGGATTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                   >  2:10079490/1‑149 (MQ=255)
               ggTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAaca                                                                                                                                 >  2:511786/1‑149 (MQ=255)
                gTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                <  2:4006349/149‑1 (MQ=255)
                gTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                <  2:5534456/149‑1 (MQ=255)
                gTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                <  1:6850388/149‑1 (MQ=255)
                gTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                <  1:6469540/149‑1 (MQ=255)
                gTTGCGTAGATGGAGCTTTCGCAAAAACATTTGCAGAAAGGGGTGAGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                <  2:7690224/149‑1 (MQ=255)
                 tttCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACAcc                                                                                                                               >  2:3628132/1‑149 (MQ=255)
                  ttCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCa                                                                                                                              <  1:2031536/149‑1 (MQ=255)
                   tCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAg                                                                                                                             >  1:6184773/1‑149 (MQ=255)
                   tCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAg                                                                                                                             <  1:5087880/149‑1 (MQ=255)
                   tCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAg                                                                                                                             <  1:4993608/149‑1 (MQ=255)
                    cGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGt                                                                                                                            <  1:10110825/149‑1 (MQ=255)
                    cGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGt                                                                                                                            <  2:8573141/149‑1 (MQ=255)
                     gTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACATCATTTGTAGCCTTTTTAGAACAGCAGtt                                                                                                                           >  1:5918374/1‑149 (MQ=255)
                      tAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTa                                                                                                                          >  1:7253631/1‑149 (MQ=255)
                      tAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTa                                                                                                                          >  2:4564669/1‑149 (MQ=255)
                       aGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaa                                                                                                                         >  2:6682238/1‑149 (MQ=255)
                        gATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCCGTTaaa                                                                                                                        >  2:1561236/1‑149 (MQ=255)
                        gATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaaa                                                                                                                        >  2:6810081/1‑149 (MQ=255)
                         aTGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAg                                                                                                                       <  1:8663619/149‑1 (MQ=255)
                             aaCTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaaaa                                                                                                                   <  2:631808/149‑1 (MQ=255)
                                   cGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTc                                                                                                             <  1:3815361/149‑1 (MQ=255)
                                     cAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCt                                                                                                           <  2:9623495/149‑1 (MQ=255)
                                     cAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCt                                                                                                           <  1:10002682/149‑1 (MQ=255)
                                       aaaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGa                                                                                                         >  1:1929589/1‑149 (MQ=255)
                                       aaaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGa                                                                                                         >  1:219403/1‑149 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                        >  2:5011044/1‑149 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                        >  1:2591631/1‑149 (MQ=255)
                                             tttGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGt                                                                                                   <  1:7419726/149‑1 (MQ=255)
                                              ttGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTc                                                                                                  >  1:7670489/1‑149 (MQ=255)
                                                 cAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCaaa                                                                                               >  2:180355/1‑149 (MQ=255)
                                                  aGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAAc                                                                                              <  2:417716/149‑1 (MQ=255)
                                                   gAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACt                                                                                             <  1:9742832/149‑1 (MQ=255)
                                                    aaaGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACtt                                                                                            <  1:9882657/149‑1 (MQ=255)
                                                     aaGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACttt                                                                                           <  1:9353377/149‑1 (MQ=255)
                                                       gtggtgGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTc                                                                                         <  2:7437146/149‑1 (MQ=255)
                                                       gtggtgGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTc                                                                                         <  2:2886211/149‑1 (MQ=255)
                                                        tggtggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCt                                                                                        <  2:8328549/149‑1 (MQ=255)
                                                        tggtggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCt                                                                                        <  1:10312024/149‑1 (MQ=255)
                                                          gtggGGTCCGCCTCTTTGGGATACCGTTCTCAGTACCAACCGCTGGTTATGGTCCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTATCTTg                                                                                      >  2:5628008/1‑149 (MQ=255)
                                                          gtggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTg                                                                                      <  2:9754631/149‑1 (MQ=255)
                                                           tggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGa                                                                                     >  2:5932300/1‑149 (MQ=255)
                                                                   gCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTc                                                                             <  1:2720847/149‑1 (MQ=255)
                                                                   gCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTc                                                                             <  2:10467106/149‑1 (MQ=255)
                                                                   gCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTc                                                                             <  1:2930837/149‑1 (MQ=255)
                                                                    ccTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCt                                                                            >  1:7847510/1‑149 (MQ=255)
                                                                       ctTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTa                                                                         >  2:2013800/1‑149 (MQ=255)
                                                                           gggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGc                                                                     >  1:4747310/1‑149 (MQ=255)
                                                                            ggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCg                                                                    <  2:7882026/149‑1 (MQ=255)
                                                                            ggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCg                                                                    <  2:3133355/149‑1 (MQ=255)
                                                                              aTACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGAt                                                                  >  1:2287840/1‑149 (MQ=255)
                                                                                 ccGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTAt                                                               <  2:2273786/149‑1 (MQ=255)
                                                                                      ctcAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCaa                                                          <  2:2455101/149‑1 (MQ=255)
                                                                                      ctcAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCaa                                                          <  2:9236639/149‑1 (MQ=255)
                                                                                       tcAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAAt                                                         <  1:5650410/149‑1 (MQ=255)
                                                                                        cAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGTCACCATTCTTTACGGCGATTATCTCAAtt                                                        >  2:6136353/1‑149 (MQ=255)
                                                                                         aGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTc                                                       <  1:1574836/149‑1 (MQ=255)
                                                                                           tACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCct                                                     >  1:5211747/1‑149 (MQ=255)
                                                                                                aCAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGTCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTggg                                                <  2:39001/149‑1 (MQ=255)
                                                                                                aCAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTggg                                                <  2:87277/149‑1 (MQ=255)
                                                                                                 catctagTTATGGTTAATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCATTTTTATAACACCAGTTAAACAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACCGCGATTATCTCAATTCCTCTTGGc                                               <  2:6460065/143‑1 (MQ=255)
                                                                                                  aGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCt                                              >  2:10273355/1‑149 (MQ=255)
                                                                                                     tGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTcc                                           <  1:1240612/149‑1 (MQ=255)
                                                                                                        ttATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAg                                        >  1:3961088/1‑149 (MQ=255)
                                                                                                         tATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCGTTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGt                                       <  2:7081447/149‑1 (MQ=255)
                                                                                                         tATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGt                                       <  1:9571330/149‑1 (MQ=255)
                                                                                                           tGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTa                                     <  1:3652530/149‑1 (MQ=255)
                                                                                                            ggTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTAt                                    >  1:357217/1‑149 (MQ=255)
                                                                                                              ttCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTg                                  >  1:317169/1‑149 (MQ=255)
                                                                                                               tCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTgg                                 >  1:6665991/1‑149 (MQ=255)
                                                                                                                cATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGa                                >  2:8147138/1‑149 (MQ=255)
                                                                                                                   cTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACCCCAGTTAAAGAAAATTATTCCTGATGTAGCCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCTTAGTTATTGGACCa                             >  1:5476465/1‑149 (MQ=255)
                                                                                                                     gTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGt                           >  2:1269445/1‑149 (MQ=255)
                                                                                                                      tAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCa                             >  2:219517/1‑146 (MQ=255)
                                                                                                                      tAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTa                          <  1:7077413/149‑1 (MQ=255)
                                                                                                                       aaTGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTaa                         <  1:4356762/149‑1 (MQ=255)
                                                                                                                         tGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATg                       <  2:2607632/149‑1 (MQ=255)
                                                                                                                              attattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCt                  <  1:8055457/149‑1 (MQ=255)
                                                                                                                               ttattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCtt                 >  1:9169156/1‑149 (MQ=255)
                                                                                                                               ttattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCtt                 <  1:301414/149‑1 (MQ=255)
                                                                                                                                tagtGCCATTACAGCATTTGCAGCCGTTTTAGAGCACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTg                <  2:4755039/146‑1 (MQ=255)
                                                                                                                                   tGCCATTACAGCATTTGCAGCCTTTCTAGAACACCAGTTAAAGAAAATTGTTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCGCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCtt             >  2:5722951/1‑149 (MQ=255)
                                                                                                                                       aTTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTAGACCAGTAATGAACCTTGCTTCAGa         <  1:3349899/149‑1 (MQ=255)
                                                                                                                                          aCAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGc      <  1:10455264/149‑1 (MQ=255)
                                                                                                                                           cAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCCATTCCTCTGGGCTTCCTCGTTATTGGACCAGTAATGCACCTTGCTTCAGATGCg     >  2:9123663/1‑149 (MQ=255)
                                                                                                                                             gCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGtt   <  1:4074645/149‑1 (MQ=255)
                                                                                                                                              cATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGTCACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAATCTTGCTTCAGATGTGttt  >  2:9429783/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCACTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTCGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTA  >  AM406671/1007142‑1007432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑