Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 1,826,407 A→G Y12Y (TAT→TAC llmg_1845 ← metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066711,826,4070AG99.7% 4252.0 / NA 1044Y12Y (TAT→TACllmg_1845metal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  major base G (465/571);  minor base C (1/2);  total (468/576)
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 1046 total aligned reads displayed.

TAATGGCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGATACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAATTAT  >  AM406671/1826261‑1826551
                                                                                                                                                  |                                                                                                                                                 
 aaTGGCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTAc                                                                                                                                                >  1:3349908/1‑149 (MQ=255)
   tGGCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTAccc                                                                                                                                              >  2:2261461/1‑149 (MQ=255)
    ggCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCt                                                                                                                                             <  2:2540173/149‑1 (MQ=255)
     gCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTc                                                                                                                                            <  2:3162977/149‑1 (MQ=255)
           aaTAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCaaa                                                                                                                                      <  1:2868991/149‑1 (MQ=255)
             tAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAaca                                                                                                                                    <  2:4138660/149‑1 (MQ=255)
              aaCCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAacac                                                                                                                                   >  1:680226/1‑149 (MQ=255)
                ccGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCt                                                                                                                                 >  2:3599606/1‑149 (MQ=255)
                ccGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCt                                                                                                                                 <  1:1711299/149‑1 (MQ=255)
                  gCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAg                                                                                                                               >  2:1941711/1‑149 (MQ=255)
                  gCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAg                                                                                                                               <  2:9438826/149‑1 (MQ=255)
                     tCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAgcgc                                                                                                                            <  2:4811106/149‑1 (MQ=255)
                     tCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAgcgc                                                                                                                            >  2:8843562/1‑149 (MQ=255)
                     tCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAgcgc                                                                                                                            <  2:4537948/149‑1 (MQ=255)
                         aCTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGtt                                                                                                                        <  2:5684413/149‑1 (MQ=255)
                         aCTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGtt                                                                                                                        <  1:5790789/149‑1 (MQ=255)
                          cTGAGTCGGTTGTAACTCATGTTCAAGCNCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTa                                                                                                                       >  1:5350628/1‑149 (MQ=255)
                            gAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAAt                                                                                                                     <  2:9319868/149‑1 (MQ=255)
                                cGGTTGTAACTCATGTTCAAGCTCAGTCACCTCCCGAGACCCGAAGTCTAAAAGCATGTTTAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCa                                                                                                                 <  2:9109596/149‑1 (MQ=255)
                                    tGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTc                                                                                                                  <  2:8125391/144‑1 (MQ=255)
                                    tGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAtttt                                                                                                             <  1:224998/149‑1 (MQ=255)
                                      tAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAtttttt                                                                                                           >  1:2344976/1‑149 (MQ=255)
                                       aaCTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAttttttt                                                                                                          <  1:5459810/149‑1 (MQ=255)
                                         cTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtct                                                                                                        <  1:633065/149‑1 (MQ=255)
                                         cTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtct                                                                                                        >  2:5884235/1‑149 (MQ=255)
                                          tCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtctc                                                                                                       >  2:3006306/1‑149 (MQ=255)
                                            aTGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtctctc                                                                                                     <  2:527529/149‑1 (MQ=255)
                                              gTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCtt                                                                                                   <  2:766031/149‑1 (MQ=255)
                                               ttCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCttt                                                                                                  <  1:991165/149‑1 (MQ=255)
                                               ttCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCttt                                                                                                  <  1:7263134/149‑1 (MQ=255)
                                               ttCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCttt                                                                                                  <  1:3960685/149‑1 (MQ=255)
                                               ttCAAGCTCAGGCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCttt                                                                                                  <  1:4116895/149‑1 (MQ=255)
                                                 cAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCt                                                                                                <  1:7686233/149‑1 (MQ=255)
                                                   aGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAt                                                                                             <  2:7960728/149‑1 (MQ=255)
                                                      tCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttt                                                                                          <  2:7441212/149‑1 (MQ=255)
                                                        aGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttt                                                                                        >  1:3786598/1‑149 (MQ=255)
                                                         gTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAttttttt                                                                                       >  1:8008156/1‑148 (MQ=255)
                                                         gTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAttttttt                                                                                       >  2:2120681/1‑148 (MQ=255)
                                                           cACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttttct                                                                                     >  1:3073321/1‑146 (MQ=255)
                                                           cACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttttct                                                                                     <  2:2876406/149‑4 (MQ=255)
                                                             ccGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttttccta                                                                                   >  2:70457/1‑144 (MQ=255)
                                                             ccGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTACCGCTGTTAATTTCATTTTTTTCTCTCTTTCCAtttttttctta                                                                                   >  1:1927982/1‑144 (MQ=255)
                                                                cACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACGCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAATACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAAtt                                                                                 <  2:5416116/149‑1 (MQ=255)
                                                                 aCGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTc                                                                                >  1:7557825/1‑149 (MQ=255)
                                                                    agaCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtat                                                                             <  1:1028485/149‑1 (MQ=255)
                                                                    agaCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtat                                                                             >  2:9541260/1‑149 (MQ=255)
                                                                     gaCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtatt                                                                            <  2:8296540/149‑1 (MQ=255)
                                                                     gaCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtatt                                                                            <  2:8733941/149‑1 (MQ=255)
                                                                      aCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtatta                                                                           <  2:64991/149‑1 (MQ=255)
                                                                         cGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAg                                                                        <  1:4876800/149‑1 (MQ=255)
                                                                          gCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtat                                                                             >  2:7298211/1‑143 (MQ=255)
                                                                             gTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCa                                                                    >  1:2444593/1‑149 (MQ=255)
                                                                              tCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCaa                                                                   >  2:3954038/1‑149 (MQ=255)
                                                                                tAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAg                                                                 >  1:4387477/1‑149 (MQ=255)
                                                                                tAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAg                                                                 <  2:4327902/149‑1 (MQ=255)
                                                                                 aaaaGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGc                                                                >  1:5989925/1‑149 (MQ=255)
                                                                                  aaaGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCa                                                               <  2:6861822/149‑1 (MQ=255)
                                                                                   aaGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCAc                                                              >  2:1440108/1‑149 (MQ=255)
                                                                                    aGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCAcc                                                             <  2:5071120/149‑1 (MQ=255)
                                                                                      cATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTGCCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCtt                                                           >  1:7072164/1‑149 (MQ=255)
                                                                                      cATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCtt                                                           >  2:5548217/1‑149 (MQ=255)
                                                                                         gTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAAt                                                        >  1:2706895/1‑149 (MQ=255)
                                                                                          ttGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATc                                                       >  1:8559709/1‑149 (MQ=255)
                                                                                           tggaGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTGCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATGCTATTATAGTTCAAAGCACCTTAATCt                                                      <  2:8626286/146‑1 (MQ=255)
                                                                                            gAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCtt                                                     <  2:6679847/149‑1 (MQ=255)
                                                                                             aaGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTa                                                    <  2:1685867/149‑1 (MQ=255)
                                                                                                ccAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTg                                                 <  2:2428665/149‑1 (MQ=255)
                                                                                                ccAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTg                                                 >  2:7774603/1‑149 (MQ=255)
                                                                                                    aCTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCt                                                                               <  1:5194862/115‑1 (MQ=255)
                                                                                                    aCTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAg                                             <  2:113666/149‑1 (MQ=255)
                                                                                                       ccGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAttttttt                                                                                       <  1:7413893/103‑2 (MQ=255)
                                                                                                        cGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGt                                         <  1:8877780/149‑1 (MQ=255)
                                                                                                           ttGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCt                                      <  1:457816/149‑1 (MQ=255)
                                                                                                            tGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCta                                     <  2:4930629/149‑1 (MQ=255)
                                                                                                                aaaGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTAtatt                                 >  1:1152804/1‑149 (MQ=255)
                                                                                                                aaaGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTAtatt                                 <  2:5998130/149‑1 (MQ=255)
                                                                                                                aaaGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTAtatt                                 >  2:3418504/1‑149 (MQ=255)
                                                                                                                   gATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATattatt                              >  1:5617115/1‑149 (MQ=255)
                                                                                                                     tAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCt                            <  1:2143104/149‑1 (MQ=255)
                                                                                                                     tAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCt                            <  2:6391270/149‑1 (MQ=255)
                                                                                                                     tAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCt                            >  1:8274478/1‑149 (MQ=255)
                                                                                                                      aagaagTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTa                           >  2:3777860/1‑149 (MQ=255)
                                                                                                                      aagaagTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTa                           >  1:8084396/1‑149 (MQ=255)
                                                                                                                      aagaagTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTa                           <  2:2678629/149‑1 (MQ=255)
                                                                                                                        gaagTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTaaa                         <  1:4758835/149‑1 (MQ=255)
                                                                                                                            ttgttCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGc                     <  1:5062039/149‑1 (MQ=255)
                                                                                                                                  cTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTAt               >  1:6551356/1‑149 (MQ=255)
                                                                                                                                  cTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTAt               >  1:953458/1‑149 (MQ=255)
                                                                                                                                   tGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATc              <  1:7832621/149‑1 (MQ=255)
                                                                                                                                    gCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATcc             <  1:3178908/149‑1 (MQ=255)
                                                                                                                                    gCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATcc             <  1:1471719/149‑1 (MQ=255)
                                                                                                                                     cATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTa                                              <  1:9182460/115‑1 (MQ=255)
                                                                                                                                      aTCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCAc           >  1:3515479/1‑149 (MQ=255)
                                                                                                                                        cACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATTTTAGTGCTAGACGTCCTATCTTATTCTAAATGGCCCTTATCCCCtt         >  2:183867/1‑149 (MQ=255)
                                                                                                                                        cACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACtt         >  1:8954382/1‑149 (MQ=255)
                                                                                                                                         aCGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTc        >  1:9338251/1‑149 (MQ=255)
                                                                                                                                         aCGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTGTTTTCTCTCTTTCTATTTTTTTCTTTATTCTATTATAGTTCAAAGCACCTTAACCTTATTGCTGGACCTCCTATATTATTCTAAATGGCATTGATTCACTTc        >  2:6263337/1‑149 (MQ=255)
                                                                                                                                           gAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCaa      <  2:847327/149‑1 (MQ=255)
                                                                                                                                            aGTGGGGTACCCTCCCCAAACACCTAGCTCTGTTAATTTTATTCTTTTCTCTCTTTCTCTTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAAt     <  1:6921996/149‑1 (MQ=255)
                                                                                                                                              tGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAATta   <  1:4294963/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAATGGCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGATACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAATTAT  >  AM406671/1826261‑1826551

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: