Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 653,674 C→T V355M (GTG→ATG)  unknown ← Fe2+ ABC transporter, permease protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1653,6740CT98.2% 124.5 / NA 54V355M (GTG→ATG) unknownFe2+ ABC transporter, permease protein 2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/32);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCAAGCCGGTGAGAATGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCACTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCGCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGCCCACGATGGCGATGAGAGT  >  1/653431‑653917
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gtCCAAGCCGGTGAGAATGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAATGCGCTGCGGGGCGGGACGCAAACAGCGATAAACGCCACttataaa                                                                                                                                                                                                                                               >  1:509357/1‑250 (MQ=255)
gtCCAAGCCGGTGAGAATGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGGACGATAACTGCGGTAAAAGCCCTtaataaa                                                                                                                                                                                                                                               >  1:544800/1‑250 (MQ=255)
        cGGTGAGAATGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATcc                                                                                                                                                                                                                                       <  1:463408/251‑1 (MQ=255)
               aaTGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATgg                                                                                                                                                                                                                                  <  1:494001/249‑1 (MQ=255)
               aaTGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATgg                                                                                                                                                                                                                                  <  1:298874/249‑1 (MQ=255)
                     cATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGc                                                                                                                                                                                                                          <  1:56528/251‑1 (MQ=255)
                     cATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGc                                                                                                                                                                                                                          <  1:216370/251‑1 (MQ=255)
                                ttttCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCa                                                                                                                                                                                                                <  1:128602/250‑1 (MQ=255)
                                   tCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGgcg                                                                                                                                                                                                            <  1:856280/251‑1 (MQ=255)
                                         aGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGGACGATAAATGCGATAAAAGCCATTAATAAAGATAAACCGATgg                                                                                                                                                                                                                                  >  1:599477/1‑223 (MQ=255)
                                                         gatgGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTTTATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGTACGATCAATGCGATAAACGCCATTAATAAAGATAAACCGATGGGGATCGGGAACACCGGGAGGGGCAAGTCGATGTAGGCTTTGaaa                                                                                                                                                                                      >  1:423796/1‑250 (MQ=255)
                                                               aacaaATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGGAAGATAAACGCGGTAAAAGCCATTTAAAAAGATAATCCGATGGGGATAGGCAAAAAAGGCAGGGGCATggccttgctctttttcacatttcc                                                                                                                                                                                 >  1:835973/1‑227 (MQ=255)
                                                                       ccccTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGGACGATAACTGCGATAAAAGCCCTTAATAAGGATAAACCGATGGGGATCGGCAACACCCGCAGGGGCAAGTCGATGTAGGGTTTGAAAAC‑GGGTAAACCggg                                                                                                                                                                       >  1:166866/1‑249 (MQ=255)
                                                                                      gAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGAt                                                                                                                                                                                                                                            <  1:787930/168‑1 (MQ=255)
                                                                                                  cTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGa                                                                                                                                                                                                                                             <  1:12131/155‑1 (MQ=255)
                                                                                                         tGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGc                                                                                                                                       <  1:717890/249‑1 (MQ=255)
                                                                                                                 cTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTgatgat                                                                                                                             <  1:712312/251‑1 (MQ=255)
                                                                                                                    gCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGGACGATAACTGCGATAAAAGCCCTTAATAAGGATAATCCGATGGTGATCGGCAacttcagcaacaatatcc                                                                                                                                                                                                       >  1:38620/1‑159 (MQ=255)
                                                                                                                     cccGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTa                                                                                                                                          <  1:389064/234‑1 (MQ=255)
                                                                                                                     cccGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGTACGATAAATGCGATAAAAGCCATTAATAAAGATAATCCGATGGGGATCGGCAACACCCGCAGGCGCATGTCGATGTAGGCTTTTAAGAA‑GGGTAAACCGAGTACTTGCGCCCGTAAAAAACCGATAAATAAACCGCTGGTGATcgcg                                                                                                                         >  1:684768/1‑251 (MQ=255)
                                                                                                                     cccGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGGACGATAACTGCGATAAAAGCCCTTAATAAGGATAATCCGATGGGGGTCGGCAACACCGGCAGGCGCATGTCGATGGAGGCCTTGAAGAC‑GGGTAAGCCGGGGACTTGCGCCAGGAAATAACCCATAAATAAACCGCTGATGAAcgcg                                                                                                                         >  1:584518/1‑251 (MQ=255)
                                                                                                                            gttttATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACGCATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGAt                                                                                                                                                                                                                                            <  1:18140/129‑1 (MQ=255)
                                                                                                                                  tataAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGAt                                                                                                                                                                                                                                            <  1:185992/124‑1 (MQ=255)
                                                                                                                                    taAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAAtt                                                                                                          <  1:206108/251‑1 (MQ=255)
                                                                                                                                            gtgtTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAg                                                                                                                                                                                      <  1:441652/168‑1 (MQ=255)
                                                                                                                                                        cGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGGACGATAACTGCGATAAACGCCATTAATAAGGATAAACCGATGGGGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGGTAAACCGAGGACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACgg                                                                                                                  >  1:761608/1‑223 (MQ=255)
                                                                                                                                                          gTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAAAGCCATTAATAAGGATAATCCGATGGGGATCGGCAACAACGGCAGGGGCATGTCGATGTAGGCTTTTAAGAC‑GGTTAAGCCGAGTACTTGCGCCCGGAAATAACCGATAATTAAACCGCTGATGAACGCGCATACGGCGATAATTAAGGTTTCCTCGCCGAtttgtt                                                                                    >  1:159033/1‑251 (MQ=255)
                                                                                                                                                          gTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGGACGATAACTGCGATAAAAGCCCTTAATAAGGATAAACCGATGGGGGTCGGCAACACCCGCAGGCGCATGTCGGTGTAGGGTTTGGAGACGGGTTAA‑CCCAGGACCTGCGCCCGTAAATAAACGATAATTAAACCGCTGGTGATCGCGCATACGGCGGTAAATTAGGGTTCCTCGCCGAtttgtt                                                                                    >  1:11922/1‑251 (MQ=255)
                                                                                                                                                          gTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGGACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGGGATCGGCAACACCGGCAGGCGCCTGTCGATGTAGGCTTTGAAGa                                                                                                                                                                                     >  1:291218/1‑155 (MQ=255)
                                                                                                                                                            tttAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGAtttgtttg                                                                                  <  1:588043/251‑1 (MQ=255)
                                                                                                                                                                 ccATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTtgatga                                                                              <  1:482525/250‑1 (MQ=255)
                                                                                                                                                                 ccATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTtgatga                                                                              <  1:472284/250‑1 (MQ=255)
                                                                                                                                                                   aTGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGtt                                                                                                    <  1:58298/226‑1 (MQ=255)
                                                                                                                                                                   aTGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGTACGATAAATGCGATAAAAGCCATTAATAAGGATAATCCGGTGGGGGTCGGCAACACCGGCAGGCGCCTGTCGAAGTAGGCTTTTAAGAC‑GGTTTAACCGAGTACTTGCGCCCGGAAATAACCGATAATTAAGCCGCTGATGATCGCGCCTACGGCGATAATTAAGGGTTCCTCGCCGATTTGTTTGATGATGt                                                                           >  1:689238/1‑251 (MQ=255)
                                                                                                                                                                     gTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGGCGGGACGATAACTTCGGTAAACGCCCTTTATAAGGATAAACCGGTGGGGATCGGCAACACCGGCAGGCGCCTGTCGGTGGAGGCTTTGAAGAC‑GGGTAAGCCCAGGACTTGCGCCCGGAAATAAACGATAATTAAGCCGCCGATGATCGCGCATACGGCGATAAt                                                                                                           >  1:776482/1‑217 (MQ=255)
                                                                                                                                                                      ttATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGa                                                                                                                                                                                                                                             <  1:463792/87‑1 (MQ=255)
                                                                                                                                                                             ccTTTTAACACATTGGCGGTTCGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTgatgat                                                                                                                             <  1:425267/191‑1 (MQ=255)
                                                                                                                                                                                     cacaTTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAAAGCCCTTAATAAGGATAATCCGATGGGGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGGAAATAAACGATAATTAAACCGCTGATGATCGCGCCTACGGCGATAATTAAGGGTTCCTCGCCGATTTGTTTGATGATGTCCCGGGTTTTTGCGCCgg                                                         >  1:596508/1‑250 (MQ=255)
                                                                                                                                                                                               gTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGGACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTg                                               >  1:456566/1‑251 (MQ=255)
                                                                                                                                                                                                       gTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATgg                                                                                                                                                                                                                                  >  1:403555/1‑65 (MQ=255)
                                                                                                                                                                                                       gTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACgg                                                                                                                  <  1:560225/176‑1 (MQ=255)
                                                                                                                                                                                                       gTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACgg                                                                                                                  <  1:629364/176‑1 (MQ=255)
                                                                                                                                                                                                               cgcgGCGGGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCtttt                               <  1:25242/251‑1 (MQ=255)
                                                                                                                                                                                                                   gcggGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGcc                            <  1:296809/250‑1 (MQ=255)
                                                                                                                                                                                                                    cggGTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCg                          <  1:392131/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggTCGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGtt                        >  1:577338/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGGTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGc                                                                                                                                                                           >  1:449869/1‑101 (MQ=255)
                                                                                                                                                                                                                           gTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGcc                    >  1:524989/1‑250 (MQ=255)
                                                                                                                                                                                                                           gTACGATAACTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTAAAGAC‑GGTTAAGCCg                                                                                                                                                                         <  1:607358/101‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGCGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAgg                                                                                                      >  1:850144/1‑159 (MQ=255)
                                                                                                                                                                                                                                       cGATAAACGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGCCCACGATGGCGATg       <  1:899857/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGCCCACGATGGCGATGagag   <  1:815354/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGCCCACGATGGCGATGAGAgt  <  1:664460/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCGCCATTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCTCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGCCCACGATGGCGATGAGAgt  <  1:391954/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GTCCAAGCCGGTGAGAATGTTCATTAAGGTGGTTTTCCCCGAGCCGGAAGCGCCCATGATGGCAACAAATTCCCCTTCGGCAATTTGAATATTAATGTCTTCAAGTGCGGTCACTTGCCCGAATCTTTTATATAAATGCTGTGTTTCAATCACGGTTTTAGCCATGTTATTCCCCTTTTAACACATTGGCGGTTTGAATGTCTAAGGCGCGGCGGGTCGGTACGATAACTGCGATAAACGCCACTAATAAGGATAATCCGATGGTGATCGGCAACACCGGCAGGCGCATGTCGATGTAGGCTTTGAAGAC‑GGTTAAGCCGAGTACTTGCGCCAGTAAATAACCGATAATTAAGCCGCTGATGATCGCGCATACGGCGATAATTAAGGTTTCCGCGCCGATTTGTTTGATGATGTCCCGGCTTTTCGCGCCGAGGGCTTTTTGCAACGCGAATTCTTTTGCCCGTTCGCCCACGATGGCGATGAGAGT  >  1/653431‑653917

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