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.6% 3828.8 / NA 953Y12Y (TAT→TACllmg_1845metal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  major base G (450/495);  minor base T (0/4);  total (453/500)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 954 total aligned reads displayed.

TAATGGCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGATACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAATTAT  >  AM406671/1826261‑1826551
                                                                                                                                                  |                                                                                                                                                 
   tGGCTCAAAATAACCGCATCCAAGTCAGTCGGTTGTAACTCATGTTCAAACTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTAccc                                                                                                                                              >  2:5011332/1‑149 (MQ=255)
           aaTAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCaaa                                                                                                                                      >  2:3187924/1‑149 (MQ=255)
             tAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAaca                                                                                                                                    <  1:8204827/149‑1 (MQ=255)
              aaCCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAacac                                                                                                                                   >  1:1096991/1‑149 (MQ=255)
               aCCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACAcc                                                                                                                                  >  1:49256/1‑149 (MQ=255)
                ccGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCt                                                                                                                                 <  1:1324372/149‑1 (MQ=255)
                            gAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAAt                                                                                                                     <  2:3483517/149‑1 (MQ=255)
                             aGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAAtt                                                                                                                    >  2:7706953/1‑149 (MQ=255)
                             aGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAAtt                                                                                                                    >  2:5942959/1‑149 (MQ=255)
                              gTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAAttt                                                                                                                   <  1:2287389/149‑1 (MQ=255)
                               tctggtGTAACTCTTGTTCAATCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATTTTGAAGCCAGACTCCGATTGTAAAATATAAGAAGTTGTTCCTGCATCACGAGTGGGGTGCCCTCCCCAAACACCTAGCGCTGTTAATTTc                                                                                                                  <  2:5603612/144‑1 (MQ=255)
                                  gTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAtt                                                                                                               <  2:7476918/149‑1 (MQ=255)
                                   ttGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAttt                                                                                                              >  1:2506244/1‑149 (MQ=255)
                                   ttGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAttt                                                                                                              >  1:4429088/1‑149 (MQ=255)
                                    tGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAtttt                                                                                                             <  1:647874/149‑1 (MQ=255)
                                      tAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAtttttt                                                                                                           <  2:2321658/149‑1 (MQ=255)
                                      tAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCAtttttt                                                                                                           <  2:1187292/149‑1 (MQ=255)
                                        aCTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtc                                                                                                         <  1:4359862/149‑1 (MQ=255)
                                        aCTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtc                                                                                                         <  1:4242912/149‑1 (MQ=255)
                                         cTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtct                                                                                                        <  2:2473944/149‑1 (MQ=255)
                                           cATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtctct                                                                                                      <  2:2701896/149‑1 (MQ=255)
                                           cATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtctct                                                                                                      <  2:4307123/149‑1 (MQ=255)
                                            aTGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAACATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTtctctc                                                                                                     >  2:6209896/1‑149 (MQ=255)
                                             tGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTNGCGCTGTTAATTTCATTTTTTTctctct                                                                                                    <  1:5086084/149‑1 (MQ=255)
                                             tGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTctctct                                                                                                    <  1:7762656/149‑1 (MQ=255)
                                             tGTTCAAGCTCAGGCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTctctct                                                                                                    <  2:6310154/149‑1 (MQ=255)
                                              gTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCtt                                                                                                   <  1:8000303/149‑1 (MQ=255)
                                                tCAAGCTCAGTAACCGCACGAGACCCGCAGTATAAAAGCATGTTGAAGCCAGACTCCGATTGTACAAGATAAGAAGTTGTTCCTGCATCTCGAGTTGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTc                                                                                                 <  2:8261397/149‑1 (MQ=255)
                                                 cAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCCTGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCt                                                                                                <  1:5653546/149‑1 (MQ=255)
                                                    gCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtt                                                                                            >  1:1006351/1‑149 (MQ=255)
                                                    gCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGAGAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTGTCTAtt                                                                                            <  2:4467014/149‑1 (MQ=255)
                                                      tCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttt                                                                                          <  1:2577962/149‑1 (MQ=255)
                                                        aGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttt                                                                                        <  1:1898327/149‑1 (MQ=255)
                                                        aGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttt                                                                                        <  1:3924167/149‑1 (MQ=255)
                                                         gTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAttttttt                                                                                       <  1:3149715/149‑2 (MQ=255)
                                                          tCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttttc                                                                                      >  1:2118483/1‑147 (MQ=255)
                                                          tCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTAtttttttc                                                                                      <  1:8329546/149‑3 (MQ=255)
                                                            aCCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCGAGCGCTGTTAACTTCATTTTTTTCTCTCTTTCTAtttttttctt                                                                                    <  1:6307354/149‑5 (MQ=255)
                                                              cGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTTGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTaa                                                                                   <  2:6540424/149‑1 (MQ=255)
                                                                cACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAAtt                                                                                 >  2:184842/1‑149 (MQ=255)
                                                                cACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAAtt                                                                                 <  1:453317/149‑1 (MQ=255)
                                                                 aCGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTc                                                                                <  2:1375579/149‑1 (MQ=255)
                                                                   gagaCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCta                                                                              <  1:3023989/149‑1 (MQ=255)
                                                                      aCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTTtttta                                                                           >  2:1595367/1‑149 (MQ=255)
                                                                      aCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtatta                                                                           <  1:879405/149‑1 (MQ=255)
                                                                      aCCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtatta                                                                           <  1:4106432/149‑1 (MQ=255)
                                                                       cccGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCtattat                                                                          >  2:7306938/1‑149 (MQ=255)
                                                                         cGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAg                                                                        <  2:2866019/149‑1 (MQ=255)
                                                                         cGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAg                                                                        >  1:3347611/1‑149 (MQ=255)
                                                                         cGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAg                                                                        <  1:477493/149‑1 (MQ=255)
                                                                            aGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGttt                                                                     >  2:816006/1‑148 (MQ=255)
                                                                              tCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCaa                                                                   >  1:6806419/1‑149 (MQ=255)
                                                                              tCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCaa                                                                   <  2:5427942/149‑1 (MQ=255)
                                                                                tAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAg                                                                 <  2:9104283/149‑1 (MQ=255)
                                                                                   aaGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCAc                                                              >  1:157798/1‑149 (MQ=255)
                                                                                   aaGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCAc                                                              <  1:5962036/149‑1 (MQ=255)
                                                                                    aGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTATCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCAcc                                                             <  2:1936359/149‑1 (MQ=255)
                                                                                    aGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCttt                                                                                                  <  1:4031953/112‑1 (MQ=255)
                                                                                    aGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCAcc                                                             <  2:7986292/149‑1 (MQ=255)
                                                                                     gCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCt                                                            >  1:6152153/1‑149 (MQ=255)
                                                                                     gCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCt                                                            >  2:6836455/1‑149 (MQ=255)
                                                                                       aTGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTa                                                          <  2:4654137/149‑1 (MQ=255)
                                                                                        tGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTaa                                                         >  2:5172641/1‑149 (MQ=255)
                                                                                         gTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAAt                                                        <  2:7610747/149‑1 (MQ=255)
                                                                                         gTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAAt                                                        >  2:4043238/1‑149 (MQ=255)
                                                                                           tGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCt                                                      <  2:3062645/149‑1 (MQ=255)
                                                                                            gAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCtt                                                     <  1:5280284/149‑1 (MQ=255)
                                                                                              aGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTATTTCTATTATAGTTCAAAGCACCTTAATCTTAg                                                   >  2:3933094/1‑149 (MQ=255)
                                                                                              aGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAg                                                   >  1:466219/1‑149 (MQ=255)
                                                                                                   gACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTa                                              <  2:4842584/149‑1 (MQ=255)
                                                                                                    aCTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAg                                             >  2:7061882/1‑149 (MQ=255)
                                                                                                      tCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGAc                                           >  1:5419648/1‑149 (MQ=255)
                                                                                                       ccGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACg                                          >  2:2160560/1‑149 (MQ=255)
                                                                                                       ccGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACg                                          <  1:84843/149‑1 (MQ=255)
                                                                                                        cGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGt                                         <  1:7268335/149‑1 (MQ=255)
                                                                                                        cGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGt                                         <  1:2897906/149‑1 (MQ=255)
                                                                                                          aTTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTcc                                       <  2:207750/149‑1 (MQ=255)
                                                                                                           ttGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCt                                      >  2:5281538/1‑149 (MQ=255)
                                                                                                            tGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCta                                     >  1:7480384/1‑149 (MQ=255)
                                                                                                             gTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCtat                                    >  1:6580828/1‑149 (MQ=255)
                                                                                                             gTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCtat                                    >  2:3743355/1‑149 (MQ=255)
                                                                                                               aaaaGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCCCCTTAATCTTAGTGCTAGACGTCCTatat                                  >  2:8411732/1‑149 (MQ=255)
                                                                                                                 aaGATAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTAtatta                                >  1:2737788/1‑149 (MQ=255)
                                                                                                                     tAAGAAGTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCCTTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTATTCTTAGTGCTAGACGTCCTATATTATTCt                            <  2:806806/149‑1 (MQ=255)
                                                                                                                       agaagTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTaa                          <  1:7060723/149‑1 (MQ=255)
                                                                                                                         aagTTGTTCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAAt                        >  1:1757211/1‑149 (MQ=255)
                                                                                                                          agTTGTTCCTGCATCACGAGTGGTGTACCCTCCCCAAACTCCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATg                       <  2:872231/149‑1 (MQ=255)
                                                                                                                           gttgttCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATgg                      <  1:3639011/149‑1 (MQ=255)
                                                                                                                                tCCTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACtt                 <  1:9082694/149‑1 (MQ=255)
                                                                                                                                 ccTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTTTTTTTTTTTTAGTTCAAATCACCTTAATCTTAGTCCTAGACGTCCTGTATTATTCTAACTGGCACttt                >  2:4122470/1‑148 (MQ=255)
                                                                                                                                  cTGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTAt               >  1:3531824/1‑149 (MQ=255)
                                                                                                                                   tGCATCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATc              >  1:8933279/1‑149 (MQ=255)
                                                                                                                                       tCACGAGTGGGTTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAATCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACt          <  2:8880587/149‑1 (MQ=255)
                                                                                                                                       tCACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACt          >  1:5796771/1‑149 (MQ=255)
                                                                                                                                        cACGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTATTTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACtt         >  2:8545575/1‑149 (MQ=255)
                                                                                                                                         aCGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTc        >  1:934779/1‑149 (MQ=255)
                                                                                                                                          cGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTTTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCATAGCACCTTAATCTTAGTGCTTGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCa       >  2:1529337/1‑149 (MQ=255)
                                                                                                                                          cGAGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCa       >  1:7658663/1‑149 (MQ=255)
                                                                                                                                            aGTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAAt     >  2:5805324/1‑149 (MQ=255)
                                                                                                                                             gTGGGGTACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAAtt    <  1:8717887/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAATGGCTCAAAATAACCGCATCCAACTGAGTCGGTTGTAACTCATGTTCAAGCTCAGTCACCGCACGAGACCCGCAGTCTAAAAGCATGTTGAAGCCAGACTCCGATTGTAAAAGATAAGAAGTTGTTCCTGCATCACGAGTGGGATACCCTCCCCAAACACCTAGCGCTGTTAATTTCATTTTTTTCTCTCTTTCTATTTTTTCTTAATTCTATTATAGTTCAAAGCACCTTAATCTTAGTGCTAGACGTCCTATATTATTCTAAATGGCACTTATCCACTTCAATTAT  >  AM406671/1826261‑1826551

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: