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% 3505.4 / NA 892F230F (TTC→TTTbglPsimilar to PTS system, beta‑glucosides specific enzyme IIABC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (443/447);  minor base G (1/0);  total (445/447)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 892 total aligned reads displayed.

GCACTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTCGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTAG  >  AM406671/1007142‑1007433
                                                                                                                                                  |                                                                                                                                                 
gcaCTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTGAGCATCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGc                                                                                                                                                 >  2:5801313/1‑149 (MQ=255)
gcaCTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGc                                                                                                                                                 <  2:8246677/149‑1 (MQ=255)
     tGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTAGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCt                                                                                                                                            >  1:4715998/1‑149 (MQ=255)
      gTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCtt                                                                                                                                           >  1:4326589/1‑149 (MQ=255)
       tttATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttt                                                                                                                                          >  1:8974886/1‑149 (MQ=255)
       tttATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttt                                                                                                                                          <  1:2726875/149‑1 (MQ=255)
         tATCCAGGTTTCGTAGATGGAACTGTCGNAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATGGCCATTACAGCATTTGCAGCCttttt                                                                                                                                        >  1:5136812/1‑149 (MQ=255)
          aTCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTa                                                                                                                                       <  1:9227979/149‑1 (MQ=255)
            ccAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGa                                                                                                                                     <  1:1824509/149‑1 (MQ=255)
            ccAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGa                                                                                                                                     <  1:5835199/149‑1 (MQ=255)
             cAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                    <  2:8740397/149‑1 (MQ=255)
             cAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                    >  1:8564220/1‑149 (MQ=255)
             cAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                    >  1:602936/1‑149 (MQ=255)
                gTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                 <  2:1688567/149‑1 (MQ=255)
                  ttCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCa                                                                                                                               <  2:7034912/149‑1 (MQ=255)
                    cGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGt                                                                                                                             <  2:7819840/149‑1 (MQ=255)
                      tAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTa                                                                                                                           >  1:6449525/1‑149 (MQ=255)
                      tAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTa                                                                                                                           <  1:2653389/149‑1 (MQ=255)
                       aGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaa                                                                                                                          >  2:1149132/1‑149 (MQ=255)
                       aGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaa                                                                                                                          >  2:1920816/1‑149 (MQ=255)
                        gATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaaa                                                                                                                         >  2:2816753/1‑149 (MQ=255)
                          tGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGa                                                                                                                       <  2:3633290/149‑1 (MQ=255)
                           ggAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaa                                                                                                                      >  2:6513069/1‑149 (MQ=255)
                            gAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTATCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCGGCTCAAGaaa                                                                                                                     >  2:4615056/1‑149 (MQ=255)
                             aaCTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaaaa                                                                                                                    <  1:4895241/149‑1 (MQ=255)
                             aaCTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaaaa                                                                                                                    >  1:7236745/1‑149 (MQ=255)
                             aaCTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaaaa                                                                                                                    <  1:4041699/149‑1 (MQ=255)
                                tGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAatta                                                                                                                 <  2:6004550/149‑1 (MQ=255)
                                cgTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGATATGGTTCATCTGTAATGCCAATGATTGCCATTACAGCATTGGCAGCCTTTTTAGAACACCAGTTAAAGAAAatta                                                                                                                 >  2:2448105/2‑149 (MQ=255)
                                 gTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAattat                                                                                                                >  1:5361821/1‑149 (MQ=255)
                                  tCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAattatt                                                                                                               >  2:1809420/1‑149 (MQ=255)
                                    gCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTCGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTcc                                                                                                             <  1:5655001/149‑1 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                         >  1:1597171/1‑149 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                         >  2:5578223/1‑149 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTACTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACGCCAGTTAAAGAAAATTATTCCTGAt                                                                                                         <  2:5032814/149‑1 (MQ=255)
                                         aaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATg                                                                                                        >  1:7901262/1‑149 (MQ=255)
                                          aCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttt                                                                                                                                          <  2:3043706/114‑1 (MQ=255)
                                          aCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGt                                                                                                       <  2:4112796/149‑1 (MQ=255)
                                              ttGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTc                                                                                                   >  2:26257/1‑149 (MQ=255)
                                                      aGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACtttt                                                                                           >  1:6998407/1‑149 (MQ=255)
                                                      aGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACtttt                                                                                           <  2:7481969/149‑1 (MQ=255)
                                                        tggtggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCt                                                                                         <  1:504994/149‑1 (MQ=255)
                                                              ggTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGacac                                                                                   <  1:3043925/149‑1 (MQ=255)
                                                               gTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACAcc                                                                                  <  1:6815309/149‑1 (MQ=255)
                                                                tAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCa                                                                                 >  2:1307579/1‑149 (MQ=255)
                                                                   gCCTCTTTGGGATACCGGTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTc                                                                              >  1:3527621/1‑149 (MQ=255)
                                                                    ccTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCt                                                                             >  2:8560991/1‑149 (MQ=255)
                                                                    ccTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCt                                                                             <  1:8376295/149‑1 (MQ=255)
                                                                     ctctTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCtt                                                                            <  1:1323399/149‑1 (MQ=255)
                                                                       ttttGGGAGACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTa                                                                          <  2:3401347/148‑1 (MQ=255)
                                                                       ctTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTa                                                                          >  2:8347019/1‑149 (MQ=255)
                                                                         ttGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACg                                                                        >  1:7098701/1‑149 (MQ=255)
                                                                         ttGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACg                                                                        >  1:7411255/1‑149 (MQ=255)
                                                                          tGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACgg                                                                       <  2:3797387/149‑1 (MQ=255)
                                                                          tGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACgg                                                                       <  1:6586643/149‑1 (MQ=255)
                                                                           gggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGc                                                                      >  2:4239623/1‑149 (MQ=255)
                                                                            ggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCg                                                                     >  2:1052943/1‑149 (MQ=255)
                                                                            ggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCg                                                                     <  1:1042595/149‑1 (MQ=255)
                                                                                 ccGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTAt                                                                <  1:8674146/149‑1 (MQ=255)
                                                                                    ttCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTAtctc                                                             >  1:4570873/1‑149 (MQ=255)
                                                                                     tctcAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCa                                                            >  1:3754347/1‑149 (MQ=255)
                                                                                     tctcAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCa                                                            >  2:8421899/1‑149 (MQ=255)
                                                                                     tctcAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAGGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCa                                                            <  2:7642290/149‑1 (MQ=255)
                                                                                          gTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTcc                                                       >  1:7649884/1‑149 (MQ=255)
                                                                                           tACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCTATTCct                                                      >  2:7118800/1‑149 (MQ=255)
                                                                                           tACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCct                                                      <  1:7801736/149‑1 (MQ=255)
                                                                                           tACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCct                                                      >  2:8922284/1‑149 (MQ=255)
                                                                                              cAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTg                                                   >  1:8155573/1‑149 (MQ=255)
                                                                                                aCAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTggg                                                 <  1:6239896/149‑1 (MQ=255)
                                                                                                 cAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGc                                                <  2:525930/149‑1 (MQ=255)
                                                                                                      ggTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCt                                           <  1:9027105/149‑1 (MQ=255)
                                                                                                       gTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTa                                          <  2:8172768/149‑1 (MQ=255)
                                                                                                        ttATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAg                                         >  1:2069519/1‑149 (MQ=255)
                                                                                                        ttATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAg                                         <  2:7299399/149‑1 (MQ=255)
                                                                                                           tGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTATTTTATTTGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTa                                      >  2:2104028/1‑149 (MQ=255)
                                                                                                             gTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTAtt                                    <  1:5469768/149‑1 (MQ=255)
                                                                                                             gTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTAtt                                    >  2:6250235/1‑149 (MQ=255)
                                                                                                              ttCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAATATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTg                                   <  2:6112538/149‑1 (MQ=255)
                                                                                                              ttCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTg                                   >  2:4838759/1‑149 (MQ=255)
                                                                                                              ttCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTg                                   <  2:917252/149‑1 (MQ=255)
                                                                                                               tCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGAGACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCGTCCTAGTTATTgg                                  >  1:4480873/1‑149 (MQ=255)
                                                                                                               tCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTgg                                  <  1:1895377/149‑1 (MQ=255)
                                                                                                                cATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGa                                 <  2:4404352/149‑1 (MQ=255)
                                                                                                                 aTCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGAc                                >  1:2208869/1‑149 (MQ=255)
                                                                                                                     gTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGt                            <  2:218447/149‑1 (MQ=255)
                                                                                                                             aaTTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAAcc                    <  1:8756475/149‑1 (MQ=255)
                                                                                                                              attattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCt                   >  2:7952423/1‑149 (MQ=255)
                                                                                                                              attattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCt                   >  2:807155/1‑149 (MQ=255)
                                                                                                                              attattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCt                   >  1:3385119/1‑149 (MQ=255)
                                                                                                                              attattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCt                   <  1:8259039/149‑1 (MQ=255)
                                                                                                                                tattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTg                 >  1:2413918/1‑149 (MQ=255)
                                                                                                                                   tGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCtt              <  2:6857181/149‑1 (MQ=255)
                                                                                                                                   tGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCtt              <  2:6358754/149‑1 (MQ=255)
                                                                                                                                      cATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAg           <  1:8066424/149‑1 (MQ=255)
                                                                                                                                       aTTACAGCATTTGCAGCCTTTTTATAATACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGa          >  2:6687400/1‑149 (MQ=255)
                                                                                                                                        ttACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATATAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTCTTGGACCAGTGATGAACCTTGCTTCAGAt         >  2:5365939/1‑149 (MQ=255)
                                                                                                                                          aCAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGc       >  1:2857216/1‑149 (MQ=255)
                                                                                                                                             gCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGtt    >  1:5258716/1‑149 (MQ=255)
                                                                                                                                              cATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTa   <  1:2936458/149‑1 (MQ=255)
                                                                                                                                              cATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTa   <  2:3167049/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCACTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTCGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTAG  >  AM406671/1007142‑1007433

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