Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 1,783,142 G→A A603V (GCA→GTA)  llmg_1803 ← putative secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066711,783,1420GA99.9% 2843.5 / NA 745A603V (GCA→GTA) llmg_1803putative secreted protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (348/394);  minor base C (0/1);  total (349/396)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 745 total aligned reads displayed.

GGATTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTGCCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCATTAGTGATAGACGT  >  AM406671/1782996‑1783287
                                                                                                                                                  |                                                                                                                                                 
 gATTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCg                                                                                                                                                <  2:4306653/149‑1 (MQ=255)
 gATTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCAATAGTGGAAGGTTACCg                                                                                                                                                >  2:3645486/1‑149 (MQ=255)
  aTTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGt                                                                                                                                               >  2:1478097/1‑149 (MQ=255)
  aTTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGt                                                                                                                                               <  2:7707088/149‑1 (MQ=255)
  aTTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGt                                                                                                                                               >  2:6835589/1‑149 (MQ=255)
   ttttATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGtt                                                                                                                                              <  1:1135768/149‑1 (MQ=255)
    tttATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGttt                                                                                                                                             >  1:2561462/1‑149 (MQ=255)
       attaGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGACGGTTACCGTTtgtg                                                                                                                                          >  2:1965223/1‑149 (MQ=255)
       attaGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTtgtg                                                                                                                                          >  1:7611143/1‑149 (MQ=255)
       attaGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTtgtg                                                                                                                                          <  1:422533/149‑1 (MQ=255)
        ttaGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGtgc                                                                                                                                         >  2:3926106/1‑149 (MQ=255)
           gCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGtgctgc                                                                                                                                      <  1:5437889/149‑1 (MQ=255)
           gCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGtgctgc                                                                                                                                      >  2:5252195/1‑149 (MQ=255)
            ccACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGtgctgct                                                                                                                                     <  2:4535354/149‑1 (MQ=255)
            ccACTGTCGCTGTGGATGCGGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGCGACGTAGCCCTAGTTGAAGGTTACCGTTTGtgctgct                                                                                                                                     >  2:2075363/1‑149 (MQ=255)
             cACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGtgctgctg                                                                                                                                    <  1:533963/149‑1 (MQ=255)
              gcTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGGCCCCCTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATCGTTGAAGGTTACCGTTTGtgctgctgc                                                                                                                                   <  2:3334169/148‑1 (MQ=255)
              aCTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGtgctgctgc                                                                                                                                   >  2:3117584/1‑149 (MQ=255)
               cTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTgctgctgct                                                                                                                                  >  1:8515597/1‑149 (MQ=255)
                tGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGctgctgctg                                                                                                                                 >  1:994742/1‑149 (MQ=255)
                  tCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTc                                                                                                                               <  2:7257908/149‑1 (MQ=255)
                    gCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAg                                                                                                                             <  1:4792378/149‑1 (MQ=255)
                     cTGTGGATGCTGTATTCCAGCCAGTCGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGc                                                                                                                            >  2:5050946/1‑149 (MQ=255)
                      tgtgGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCg                                                                                                                           >  2:2584918/1‑149 (MQ=255)
                      tgtgGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACGCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCg                                                                                                                           <  1:3575391/149‑1 (MQ=255)
                         ggATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGa                                                                                                                        <  2:7440766/149‑1 (MQ=255)
                         ggATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGa                                                                                                                        >  1:4355910/1‑149 (MQ=255)
                         ggATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGa                                                                                                                        >  2:9069028/1‑149 (MQ=255)
                         ggATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGa                                                                                                                        <  1:3823017/149‑1 (MQ=255)
                          gATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGAc                                                                                                                       >  1:4938414/1‑149 (MQ=255)
                            tGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCg                                                                                                                     >  1:2727298/1‑149 (MQ=255)
                             gCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGt                                                                                                                    >  1:1910059/1‑149 (MQ=255)
                             gCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGt                                                                                                                    <  1:6213514/149‑1 (MQ=255)
                             gCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGt                                                                                                                    >  1:7460538/1‑149 (MQ=255)
                              cTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGta                                                                                                                   <  1:7257869/149‑1 (MQ=255)
                               tGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGtat                                                                                                                  <  1:874211/149‑1 (MQ=255)
                               tGTATTCCAGCCACTTGGTGAAGTAACCTTTGCCCCTGCCGTCCGACTCGCTCATAAGCTAACAGGTATAGCAGCCGTGCTGAATGCTCCTGAGCCGTAGCAACAGTTGACAGTTACCGTTTGTGCTGCTGCTGTCACCGTGTCCGtat                                                                                                                  >  2:227455/1‑149 (MQ=255)
                                 tATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAg                                                                                                                <  2:6342320/149‑1 (MQ=255)
                                  aTTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGt                                                                                                               >  1:4427497/1‑149 (MQ=255)
                                      cagccagTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTAtt                                                                                                           >  1:2972169/1‑149 (MQ=255)
                                        gccagTTGGTGGTGTGACTTTTGCTCCTGTCTTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCTTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGt                                                                                                         <  2:5536842/149‑1 (MQ=255)
                                         ccagTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGtt                                                                                                        >  1:8415971/1‑149 (MQ=255)
                                            gTTGGTGGTGTGACATTTGCTCCTGTCGTCCCGCTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGc                                                                                                     <  1:6954673/149‑1 (MQ=255)
                                                 ttgtgtGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTggg                                                                                                <  2:3496941/147‑1 (MQ=255)
                                                 tggtgTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTggg                                                                                                <  1:3261311/149‑1 (MQ=255)
                                                  ggtgTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGt                                                                                               <  1:9114451/149‑1 (MQ=255)
                                                  ggtgTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGATTGGGt                                                                                               >  2:3219576/1‑149 (MQ=255)
                                                   gtgtGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGACGGTTACCGTTTGTGCTACTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGtt                                                                                              >  1:101678/1‑149 (MQ=255)
                                                       gACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGGGCTCCTGAGTCGTAGCCATCGTTGAAGGTTACCGTTTGTGGTGCGGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTgcgc                                                                                          <  2:6932407/149‑1 (MQ=255)
                                                          aTTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTaa                                                                                       <  1:7753159/149‑1 (MQ=255)
                                                                    gtcgtcCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTgttgggtt                                                                             >  2:2850132/1‑149 (MQ=255)
                                                                      cgtcCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTgttgggttgg                                                                           >  1:7881321/1‑149 (MQ=255)
                                                                         cccACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCg                                                                        <  2:2726576/149‑1 (MQ=255)
                                                                         cccACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCg                                                                        >  2:8167159/1‑149 (MQ=255)
                                                                          ccACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGc                                                                       >  1:5818536/1‑149 (MQ=255)
                                                                           cACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCt                                                                      >  1:3131892/1‑149 (MQ=255)
                                                                            aCTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCtg                                                                     <  1:8632746/149‑1 (MQ=255)
                                                                                      aGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGc                                                           >  1:8900385/1‑149 (MQ=255)
                                                                                        gCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTCACCGTATAGTTATTGTTTGCTTGGGTTTCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCa                                                         <  2:6143527/149‑1 (MQ=255)
                                                                                              agGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCg                                                   <  1:6533725/149‑1 (MQ=255)
                                                                                                 tataGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTg                                                >  2:5412608/1‑149 (MQ=255)
                                                                                                  ataGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGc                                               >  2:9159872/1‑149 (MQ=255)
                                                                                                  ataGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGc                                               <  1:4582783/149‑1 (MQ=255)
                                                                                                    aGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAAGTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAg                                             <  2:58190/149‑1 (MQ=255)
                                                                                                    aGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAg                                             <  2:7959302/149‑1 (MQ=255)
                                                                                                    aGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAg                                             >  1:4013214/1‑149 (MQ=255)
                                                                                                       aGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTg                                          >  2:8881002/1‑149 (MQ=255)
                                                                                                       aGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTg                                          <  1:670663/149‑1 (MQ=255)
                                                                                                       aGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTg                                          <  1:8042825/149‑1 (MQ=255)
                                                                                                        ggCGTGTTGAGTGCCCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGt                                         <  1:6395296/149‑1 (MQ=255)
                                                                                                         gCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCAAGTGTCGCTGCAGTTGTc                                        >  2:1656346/1‑149 (MQ=255)
                                                                                                             gtTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCAc                                    <  1:8164855/149‑1 (MQ=255)
                                                                                                              ttGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAANTTGCCAGTGTCGCTGCAGTTGTCGCAcc                                   <  1:5099291/149‑1 (MQ=255)
                                                                                                              ttGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCAcc                                   <  2:6568452/149‑1 (MQ=255)
                                                                                                               tGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCa                                  <  2:4901212/149‑1 (MQ=255)
                                                                                                                gAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAg                                 >  1:2419695/1‑149 (MQ=255)
                                                                                                                gAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAg                                 >  1:7395618/1‑149 (MQ=255)
                                                                                                                 aGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGt                                >  1:5954466/1‑149 (MQ=255)
                                                                                                                 aGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGt                                >  1:1519997/1‑149 (MQ=255)
                                                                                                                   tGCTCCTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGtat                              <  2:3845286/149‑1 (MQ=255)
                                                                                                                       ccTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGt                          <  2:8779993/149‑1 (MQ=255)
                                                                                                                        cTGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTa                         <  2:2211581/149‑1 (MQ=255)
                                                                                                                         tGAGACGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTAc                        <  2:8589830/149‑1 (MQ=255)
                                                                                                                          tataCGTAGCGATAGTTGAAGGTTACCTTTTGTGCTGCTGCTTTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTTCGCTAACTGTTCGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTTCATTTGTCGCACCAGTATAAGTACt                       <  2:1242725/146‑1 (MQ=255)
                                                                                                                          gagaCGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACt                       >  1:6113090/1‑149 (MQ=255)
                                                                                                                          gagaCGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACt                       <  1:5658547/149‑1 (MQ=255)
                                                                                                                          gagaCGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACt                       >  1:3435528/1‑149 (MQ=255)
                                                                                                                           ggaCGTAGCCATAGTTGAAGGTTACCGTTTATGGTGCTGCTTGCAGCGTGTCCGTCTTGGTATTGTTTGCCTGGGTTGCGGTAACTGTTGGGTTGGCCGCTGTGAATGGTGCCAGTGTCGCTGCAGTTGTCGCTCCAGTATAAGTtctc                      <  2:4663205/148‑1 (MQ=255)
                                                                                                                            gaCGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCg                     >  2:841173/1‑149 (MQ=255)
                                                                                                                             aCGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGt                    >  2:4115593/1‑149 (MQ=255)
                                                                                                                             aCGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGt                    <  2:8346640/149‑1 (MQ=255)
                                                                                                                              cGTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTg                   <  2:8497677/149‑1 (MQ=255)
                                                                                                                               gTAGCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGc                  >  1:5546262/1‑149 (MQ=255)
                                                                                                                                  gCCATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCCCTTTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCAt               <  2:6708684/149‑1 (MQ=255)
                                                                                                                                   ccATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCAtt              <  1:1735465/149‑1 (MQ=255)
                                                                                                                                   ccATAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCAtt              <  1:1314008/149‑1 (MQ=255)
                                                                                                                                     aTAGTTGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCATTAg            <  1:5140205/149‑1 (MQ=255)
                                                                                                                                         ttGAAGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCATTAGTGAt        <  1:5234204/149‑1 (MQ=255)
                                                                                                                                            aaGGTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCATTAGTGATAGa     <  2:707467/149‑1 (MQ=255)
                                                                                                                                               gTTACCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCATTAGTGATAGACGt  >  1:2049708/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGATTTTATTAGCCACTGTCGCTGTGGATGCTGTATTCCAGCCAGTTGGTGGTGTGACATTTGCTCCTGTCGTCCCACTCGCTGACAGGCTAGTAGGTATAGCAGGCGTGTTGAGTGCTCCTGAGACGTAGCCATAGTTGAAGGTTGCCGTTTGTGCTGCTGCTGTCAGCGTGACCGTATAGTTATTGTTTGCTTGGGTTGCGCTAACTGTTGGGTTGGCCGCTGTGAAAGTTGCCAGTGTCGCTGCAGTTGTCGCACCAGTATAAGTACTCGTGCCATTAGTGATAGACGT  >  AM406671/1782996‑1783287

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