Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 970,845 C→A A78D (GCT→GAT)  llmg_1005 → sugar kinase and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM406671970,8450CA99.8% 3276.8 / NA 846A78D (GCT→GAT) llmg_1005sugar kinase and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (386/456);  minor base G (2/0);  total (390/456)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 859 total aligned reads displayed.

TGGTACTCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGCTCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGTCCTTTTA  >  AM406671/970699‑970990
                                                                                                                                                  |                                                                                                                                                 
 ggTACTCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGatct                                                                                                                                                >  2:3460584/1‑149 (MQ=255)
 ggTACTCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGatct                                                                                                                                                <  1:274609/149‑1 (MQ=255)
    aatCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGTAGTAGTTGAGCAATTTAAACCAAATTCTGATGCTTGTGGAGTACCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTccc                                                                                                                                             <  2:6943866/147‑1 (MQ=255)
    aCTCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTccc                                                                                                                                             <  1:2963774/149‑1 (MQ=255)
     cTCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCt                                                                                                                                            <  1:1558320/149‑1 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           <  2:9269268/149‑1 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           >  2:215322/1‑149 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           <  1:9182067/149‑1 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           <  1:7058906/149‑1 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           <  1:2035633/149‑1 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           >  1:5190221/1‑149 (MQ=255)
      tCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCtt                                                                                                                                           >  2:5702347/1‑149 (MQ=255)
           aTAATTTAGAACAGTTTTATCAGTCCCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTACTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCACTTATAtt                                                                                                                                      >  2:5238798/1‑149 (MQ=255)
            tAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTc                                                                                                                                     >  2:6694417/1‑149 (MQ=255)
              aTTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAACATTCTGATGTTTGTGAAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTCGCGTGCTTGAAGGAGCAAGTGATCTCCCTTATATTCAt                                                                                                                                   >  2:2968134/1‑149 (MQ=255)
               tttAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGANGGAGCAAGTGATCTCCCTTATATTCATg                                                                                                                                  <  1:5303533/149‑1 (MQ=255)
               tttAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATg                                                                                                                                  <  1:2530100/149‑1 (MQ=255)
                ttAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGa                                                                                                                                 >  2:2709980/1‑149 (MQ=255)
                 tAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGAt                                                                                                                                <  2:4275284/149‑1 (MQ=255)
                 tAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGAt                                                                                                                                >  2:9143260/1‑149 (MQ=255)
                  aGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGAtt                                                                                                                               <  1:1309197/149‑1 (MQ=255)
                     aCAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTg                                                                                                                            >  1:4225452/1‑149 (MQ=255)
                     aCAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTg                                                                                                                            <  2:6082088/149‑1 (MQ=255)
                       aGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGaa                                                                                                                          >  1:4328752/1‑149 (MQ=255)
                        gTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGaaa                                                                                                                         >  1:3938549/1‑149 (MQ=255)
                          tttATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAAtt                                                                                                                       <  1:8358456/149‑1 (MQ=255)
                           ttATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTc                                                                                                                      >  1:4922797/1‑149 (MQ=255)
                             aTCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCAt                                                                                                                    >  2:7385415/1‑149 (MQ=255)
                             aTCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCAt                                                                                                                    <  2:8350020/149‑1 (MQ=255)
                               cAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGc                                                                                                                  <  2:3312116/149‑1 (MQ=255)
                               cAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGc                                                                                                                  <  1:7511331/149‑1 (MQ=255)
                                aGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCt                                                                                                                 >  2:3714494/1‑149 (MQ=255)
                                  tctcTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGa                                                                                                               >  2:8950635/1‑149 (MQ=255)
                                    tctAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAAc                                                                                                             >  2:2421877/1‑149 (MQ=255)
                                         cGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGaa                                                                                                        <  1:8953631/149‑1 (MQ=255)
                                           gCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGaaaa                                                                                                      <  2:100341/149‑1 (MQ=255)
                                             agtagtTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCCGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGaaaaaa                                                                                                    >  1:5761334/1‑148 (MQ=255)
                                             agtagtTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAAc                                                                                                    >  1:7904651/1‑149 (MQ=255)
                                               tagtTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGc                                                                                                  >  1:6864855/1‑149 (MQ=255)
                                               tagtTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGc                                                                                                  <  1:5940665/149‑1 (MQ=255)
                                                agtTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCt                                                                                                 >  2:1788795/1‑149 (MQ=255)
                                                  ttGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCttt                                                                                               >  1:7311978/1‑149 (MQ=255)
                                                  ttGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCttt                                                                                               <  2:552939/149‑1 (MQ=255)
                                                   tGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTg                                                                                              >  1:767172/1‑149 (MQ=255)
                                                     aGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGt                                                                                            <  1:6215568/149‑1 (MQ=255)
                                                     aGCAATTTAAAGAAAAGTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGt                                                                                            <  1:6983661/149‑1 (MQ=255)
                                                      gCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTc                                                                                           >  2:7495562/1‑149 (MQ=255)
                                                      gCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTc                                                                                           <  1:3848614/149‑1 (MQ=255)
                                                         aTTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTAc                                                                                        >  1:1656959/1‑149 (MQ=255)
                                                          tttAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTAcc                                                                                       >  2:7248722/1‑149 (MQ=255)
                                                             aaaGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGt                                                                                    >  2:4632628/1‑149 (MQ=255)
                                                             aaaGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGt                                                                                    >  2:9492260/1‑149 (MQ=255)
                                                             aaaGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGt                                                                                    >  1:2361201/1‑149 (MQ=255)
                                                                gAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATc                                                                                 >  1:54878/1‑149 (MQ=255)
                                                                  aaaTTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGa                                                                               >  2:1292285/1‑149 (MQ=255)
                                                                   aaTTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATACATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGAt                                                                              <  2:6492672/149‑1 (MQ=255)
                                                                     ttCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTg                                                                            >  1:824483/1‑149 (MQ=255)
                                                                      tCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGa                                                                           >  2:7799012/1‑149 (MQ=255)
                                                                      tCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGa                                                                           <  1:3464702/149‑1 (MQ=255)
                                                                       cTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGaa                                                                          <  2:855499/149‑1 (MQ=255)
                                                                          aTGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAat                                                                       <  2:1450836/149‑1 (MQ=255)
                                                                             tttGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAatgat                                                                    >  1:3180412/1‑149 (MQ=255)
                                                                              ttGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAatgatg                                                                   <  1:9371852/149‑1 (MQ=255)
                                                                                 tgGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCa                                                                <  2:5075228/149‑1 (MQ=255)
                                                                                    aGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAAtt                                                             <  1:3705233/149‑1 (MQ=255)
                                                                                    aGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAAtt                                                             <  2:5962678/149‑1 (MQ=255)
                                                                                      tAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTgc                                                           <  1:5613175/149‑1 (MQ=255)
                                                                                      tAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTgc                                                           >  1:7990452/1‑149 (MQ=255)
                                                                                       aGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTgcg                                                          <  1:4683973/149‑1 (MQ=255)
                                                                                        gCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGa                                                                               <  1:3307595/127‑1 (MQ=255)
                                                                                         cAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGcgcg                                                        <  2:7935159/149‑1 (MQ=255)
                                                                                         cAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGcgcg                                                        >  2:3155191/1‑149 (MQ=255)
                                                                                             tAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGccc                                                    <  2:2021614/149‑1 (MQ=255)
                                                                                                ctctCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCg                                                 <  1:658415/149‑1 (MQ=255)
                                                                                                  ctcCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCt                                               >  1:4455424/1‑149 (MQ=255)
                                                                                                   tcCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTg                                              <  1:8793432/149‑1 (MQ=255)
                                                                                                      tGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAgg                                           <  1:6444881/149‑1 (MQ=255)
                                                                                                        gAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTa                                         >  1:6355095/1‑149 (MQ=255)
                                                                                                         aGCAGTTAACACGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTaa                                        >  1:2104317/1‑149 (MQ=255)
                                                                                                          gCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTGGAAAAACGCTTTGGTCTACCAGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTaaa                                       <  2:4905174/149‑1 (MQ=255)
                                                                                                            agctAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTGCCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAAt                                     <  2:5390797/146‑1 (MQ=255)
                                                                                                              ttAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAAttt                                   <  1:8532387/149‑1 (MQ=255)
                                                                                                              ttAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAAttt                                   >  1:1486617/1‑149 (MQ=255)
                                                                                                               tAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTg                                  >  1:5090058/1‑149 (MQ=255)
                                                                                                               tAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTg                                  >  2:7134389/1‑149 (MQ=255)
                                                                                                                     ggCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCgg                            >  2:2590927/1‑149 (MQ=255)
                                                                                                                      gCTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGc                           >  1:1230075/1‑149 (MQ=255)
                                                                                                                       cTTCTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGGTGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCa                          >  2:6306326/1‑149 (MQ=255)
                                                                                                                          cTGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTTCGGCaaaa                       >  2:1653742/1‑149 (MQ=255)
                                                                                                                           tGGCGTGATTGAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGACAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAg                      >  2:806557/1‑149 (MQ=255)
                                                                                                                                     gAAGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGa            <  2:8808571/149‑1 (MQ=255)
                                                                                                                                      aaGGAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGAt           >  2:386117/1‑149 (MQ=255)
                                                                                                                                        ggAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGt         >  2:7689661/1‑149 (MQ=255)
                                                                                                                                        ggAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGt         >  2:8197483/1‑149 (MQ=255)
                                                                                                                                        ggAGCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGt         <  2:413445/149‑1 (MQ=255)
                                                                                                                                           gCAAGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGTCCt      <  1:1788795/149‑1 (MQ=255)
                                                                                                                                            cAAGTGATCTCCCTTATATTCATGCTTTTGAAATTCATGCTGAACTCGAAAAACGCTTTGGTCTACCCGTACCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAACTTCGTGCGGCAAAAGGTCTCTCAGGTGTCCtt     >  2:4496139/1‑149 (MQ=255)
                                                                                                                                             aaGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGTCCttt    >  2:2244406/1‑149 (MQ=255)
                                                                                                                                             aaGTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGTCCttt    >  1:1077987/1‑149 (MQ=255)
                                                                                                                                               gTGATCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGTCCTTTTa  >  2:3992563/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGTACTCCTGATAATTTAGAACAGTTTTATCAGTCTCTAACGGCAGTAGTTGAGCAATTTAAAGAAAATTCTGATGTTTGTGGAGTAGCAATTAGCTCTCCTGGAGCAGTTAACAAGGCTTCTGGCGTGATTGAAGGAGCAAGTGCTCTCCCTTATATTCATGATTTTGAAATTCATGCTGAACTTGAAAAACGCTTTGGTCTACCCGTATCGATTGAAAATGATGCCAATTGCGCGGCCCTCGCTGAGGTAAAATTTGGTGCGGCAAAAGGTCTCTCAGATGTCCTTTTA  >  AM406671/970699‑970990

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: