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

Only 100 of 843 total aligned reads displayed.

GCACTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTCGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTA  >  AM406671/1007142‑1007432
                                                                                                                                                  |                                                                                                                                                
gcaCTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGc                                                                                                                                                >  2:948004/1‑149 (MQ=255)
    ttGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGcc                                                                                                                                            <  1:8070972/149‑1 (MQ=255)
     tGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCt                                                                                                                                           >  1:7108174/1‑149 (MQ=255)
      gTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCtt                                                                                                                                          >  1:3830490/1‑149 (MQ=255)
       tttATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttt                                                                                                                                         <  2:1232295/149‑1 (MQ=255)
        ttATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCtttt                                                                                                                                        <  1:546244/149‑1 (MQ=255)
         tATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttttt                                                                                                                                       <  2:3944244/149‑1 (MQ=255)
         tATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTAGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCttttt                                                                                                                                       <  2:7800053/149‑1 (MQ=255)
          aTCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTa                                                                                                                                      <  1:948443/149‑1 (MQ=255)
             cAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                   <  1:6572868/149‑1 (MQ=255)
             cAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGaa                                                                                                                                   <  1:8845774/149‑1 (MQ=255)
               ggTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAaca                                                                                                                                 <  2:3239430/149‑1 (MQ=255)
               ggTTACGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAaca                                                                                                                                 <  2:2879132/149‑1 (MQ=255)
                gTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAacac                                                                                                                                <  1:4269304/149‑1 (MQ=255)
                   tCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCTTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAg                                                                                                                             <  2:3006190/149‑1 (MQ=255)
                      tAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTa                                                                                                                          <  2:4157807/149‑1 (MQ=255)
                       aGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaa                                                                                                                         >  2:7962686/1‑149 (MQ=255)
                        gATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATTCCGTTCTCAGTACCAACAGCTGGTTATGGTTTATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTaaa                                                                                                                        <  2:4559143/149‑1 (MQ=255)
                         aTGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAg                                                                                                                       >  1:8958239/1‑149 (MQ=255)
                         aTGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAg                                                                                                                       >  1:8601292/1‑149 (MQ=255)
                         aTGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAg                                                                                                                       <  1:3547340/149‑1 (MQ=255)
                          tGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACTCCAGTTAAAGa                                                                                                                      <  1:4895545/149‑1 (MQ=255)
                          tGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGa                                                                                                                      <  1:8541746/149‑1 (MQ=255)
                          tGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGa                                                                                                                      <  2:6884445/149‑1 (MQ=255)
                           ggAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaa                                                                                                                     >  2:7556771/1‑149 (MQ=255)
                           ggAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaa                                                                                                                     <  2:6720806/149‑1 (MQ=255)
                           ggAACTGTCGCAAAAACATTTGCAGAAAGAGGTCGGGCAAGCCTCTTTCGGACACCGCAATCAGTACCAACAGCTGGCCATGGCTCCTCTGTCCTCCCCATCATTGCCATTACAGCCTTTGCCGCCTTTTTAGAACACCAGTTAAAGaa                                                                                                                     >  2:8441942/1‑149 (MQ=255)
                            gAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGaaa                                                                                                                    >  1:1163770/1‑149 (MQ=255)
                                 gTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAattat                                                                                                               >  1:1772700/1‑149 (MQ=255)
                                  tCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAattatt                                                                                                              <  1:6701816/149‑1 (MQ=255)
                                  tCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAattatt                                                                                                              <  2:3449111/149‑1 (MQ=255)
                                    gCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTcc                                                                                                            <  2:9156091/149‑1 (MQ=255)
                                     cAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCt                                                                                                           <  1:781023/149‑1 (MQ=255)
                                       aaaaCATTTGCAGAAAGTGGTGGTGTAAGCCTCTTTGGGATACCGTTCTCAGTACCATCGGCTGGTTATGGTGCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAGGAAAATTATTCCTGa                                                                                                         <  2:2388216/149‑1 (MQ=255)
                                       aaaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAatta                                                                                                                <  1:6073116/142‑1 (MQ=255)
                                        aagcATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                        <  2:3130313/146‑1 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                        >  2:1669388/1‑149 (MQ=255)
                                        aaaCATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGAt                                                                                                        >  1:6970940/1‑149 (MQ=255)
                                           cATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTCTTGCCATTACAGCATTTGCAGCGTTTTTAGACCACCAGTGAAAGAAAATTATTCCTGATGTa                                                                                                     <  2:2039534/149‑1 (MQ=255)
                                             tttGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGt                                                                                                   >  1:148847/1‑149 (MQ=255)
                                             tttGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGt                                                                                                   >  1:2391452/1‑149 (MQ=255)
                                                   gAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACt                                                                                             <  1:9319271/149‑1 (MQ=255)
                                                    aaaGTGGTGGGGTAAGCCTCTTTGTGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACtt                                                                                            >  1:3447911/1‑149 (MQ=255)
                                                      aGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACtttt                                                                                          >  1:248672/1‑149 (MQ=255)
                                                      aGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACtttt                                                                                          <  2:6598540/149‑1 (MQ=255)
                                                        tggtggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCt                                                                                        >  2:7409664/1‑149 (MQ=255)
                                                           tggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACGGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGGAGCCTTTTTATAACACGAGTTAAAGAAAATGATTCCTGATGTAGTCAAACTTTACTTGa                                                                                     >  2:1824588/1‑149 (MQ=255)
                                                           tggGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGa                                                                                     >  1:4024233/1‑149 (MQ=255)
                                                             gggTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGaca                                                                                   <  1:1000964/149‑1 (MQ=255)
                                                              ggTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGacac                                                                                  <  1:7285222/149‑1 (MQ=255)
                                                               gTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACAcc                                                                                 <  1:3034117/149‑1 (MQ=255)
                                                                 aaGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCAt                                                                               >  2:7116436/1‑149 (MQ=255)
                                                                   gCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTc                                                                             >  1:3177238/1‑149 (MQ=255)
                                                                    ccTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCt                                                                            <  1:2782455/149‑1 (MQ=255)
                                                                    ccTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATCCCTGATGTAGTCAAACTTTTCTTGACACCATTCt                                                                            >  1:2933929/1‑149 (MQ=255)
                                                                      tctTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCttt                                                                          <  2:5318349/149‑1 (MQ=255)
                                                                       ctTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTa                                                                         >  2:547708/1‑149 (MQ=255)
                                                                        tttGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTAc                                                                        <  2:8220850/149‑1 (MQ=255)
                                                                        tttGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTAc                                                                        <  2:8682995/149‑1 (MQ=255)
                                                                         ttGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACg                                                                       >  2:5673781/1‑149 (MQ=255)
                                                                            ggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACgg                                                                      >  2:1152787/1‑147 (MQ=255)
                                                                            ggATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCg                                                                    >  1:4659599/1‑149 (MQ=255)
                                                                             gATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGa                                                                   <  1:6191556/149‑1 (MQ=255)
                                                                             gATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGa                                                                   >  2:5070677/1‑149 (MQ=255)
                                                                              aTACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGAt                                                                  <  2:8094995/149‑1 (MQ=255)
                                                                                aCCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTa                                                                <  2:4777892/149‑1 (MQ=255)
                                                                                   gTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTCTTCGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCACACTCCTCTTGACACCATTCTTTACCCCGatcatcc                                                             >  2:1555376/1‑148 (MQ=255)
                                                                                         aGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTc                                                       >  1:131913/1‑149 (MQ=255)
                                                                                            aCCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCctt                                                    >  2:375484/1‑148 (MQ=255)
                                                                                            aCCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCctc                                                    <  1:7440373/149‑1 (MQ=255)
                                                                                             ccAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCctct                                                   <  1:1627690/149‑1 (MQ=255)
                                                                                                  aGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCt                                              >  1:5162659/1‑149 (MQ=255)
                                                                                                    cTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTc                                            >  2:2552748/1‑149 (MQ=255)
                                                                                                      ggTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCt                                          <  2:3513351/149‑1 (MQ=255)
                                                                                                         tATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGt                                       <  2:2213684/149‑1 (MQ=255)
                                                                                                            ggTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTAt                                    <  2:4326700/149‑1 (MQ=255)
                                                                                                            ggTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTAt                                    <  2:9414700/149‑1 (MQ=255)
                                                                                                             gTTCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTAtt                                   >  1:577306/1‑149 (MQ=255)
                                                                                                              ttCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTg                                  >  1:2646551/1‑149 (MQ=255)
                                                                                                              ttCATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATGATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTg                                  >  2:6052338/1‑149 (MQ=255)
                                                                                                                cATCTGTAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTgaa                                >  2:5844908/1‑147 (MQ=255)
                                                                                                                      tAATGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTa                          >  1:6803288/1‑149 (MQ=255)
                                                                                                                      tAATGCCAATTATTGCCATTACAGCATTTGCAGCCTCTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTa                          >  2:6426934/1‑149 (MQ=255)
                                                                                                                        aTGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAAt                        >  1:2213200/1‑149 (MQ=255)
                                                                                                                         tGCCAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATg                       >  1:5915032/1‑149 (MQ=255)
                                                                                                                            cAATTATTGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAAc                    >  1:7735892/1‑149 (MQ=255)
                                                                                                                               ttattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCtt                 <  1:9188590/149‑1 (MQ=255)
                                                                                                                               ttattGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCtt                 <  1:7586801/149‑1 (MQ=255)
                                                                                                                                   tGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTNGTTATTGGACCAGTAATGAACCTTGCtt             <  1:5358835/149‑1 (MQ=255)
                                                                                                                                   tGCCATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCtt             >  2:6175013/1‑149 (MQ=255)
                                                                                                                                     ccATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTa                                     <  1:6425845/123‑1 (MQ=255)
                                                                                                                                     ccATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCa           <  1:200260/149‑1 (MQ=255)
                                                                                                                                     ccATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCa           >  2:4412955/1‑149 (MQ=255)
                                                                                                                                      cATTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAg          <  2:5210005/149‑1 (MQ=255)
                                                                                                                                       aTTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGa         <  1:842499/149‑1 (MQ=255)
                                                                                                                                       aTTACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGa         >  1:3332283/1‑149 (MQ=255)
                                                                                                                                        ttACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGAt        >  2:2716325/1‑149 (MQ=255)
                                                                                                                                         tACAGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATg       <  2:8955455/149‑1 (MQ=255)
                                                                                                                                           cAGCATTTGCAGCCTTTTTTGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCg     <  1:4751872/149‑1 (MQ=255)
                                                                                                                                            aGCATTTGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGt    >  1:4385686/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCACTTGTTTATCCAGGTTTCGTAGATGGAACTGTCGCAAAAACATTTGCAGAAAGTGGTGGGGTAAGCCTCTTTGGGATACCGTTCTCAGTACCAACAGCTGGTTATGGTTCATCTGTAATGCCAATTATTGCCATTACAGCATTCGCAGCCTTTTTAGAACACCAGTTAAAGAAAATTATTCCTGATGTAGTCAAACTTTTCTTGACACCATTCTTTACGGCGATTATCTCAATTCCTCTGGGCTTCCTAGTTATTGGACCAGTAATGAACCTTGCTTCAGATGCGTTA  >  AM406671/1007142‑1007432

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