Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 332,713 C→T S11F (TCC→TTC)  unknown → Mn2+ and Fe2+ transporters of the NRAMP family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1332,7130CT100.0% 113.5 / NA 43S11F (TCC→TTC) unknownMn2+ and Fe2+ transporters of the NRAMP family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/21);  total (22/21)

TGCCGGTTATTTATTTTTTATCTTAATATTATTTCCTGAAATTATTAAATAATCAGAGTTTTATTCTGCTATTTTCTCTGAGGTAGAATTTTATTTTCTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATCCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAGGGTATCGCCAAAAAGGGCGTTTCTATT  >  1/332510‑332952
                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
tGCCGGTTATTTATTTTTTATCTTAATATTATTTCCTGAAATTATTAAATAATCAGAGTTTTATTCTGCTATTTTCTCTGAGGTAGAATTTTATTTTCTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATgg                                                                                                                                                                                                  <  1:249026/251‑1 (MQ=255)
                                                                          tctcTGAGGTAGAATTTTATTTTCTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTg                                                                                                                        <  1:755800/251‑1 (MQ=255)
                                                                          tctcTGAGGTAGAATTTTATTTTCTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTg                                                                                                                        <  1:491915/251‑1 (MQ=255)
                                                                                       attttattttCTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATcc                                                                                                           <  1:366687/251‑1 (MQ=255)
                                                                                                 cTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAAttt                                                                                                 <  1:640517/251‑1 (MQ=255)
                                                                                                             gTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATcc                                                                                                           <  1:133418/229‑1 (MQ=255)
                                                                                                                     cTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCAtattat                                                                                                                                                                   <  1:558845/165‑1 (MQ=255)
                                                                                                                     cTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTc                                                                              >  1:115955/1‑250 (MQ=255)
                                                                                                                       attatAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGcc                                                                                                                                                   >  1:47903/1‑179 (MQ=255)
                                                                                                                       attatAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGa                                                                           >  1:413391/1‑251 (MQ=255)
                                                                                                                              aTGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGgtg                                                                    <  1:865850/251‑1 (MQ=255)
                                                                                                                                          gCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGAtt                                                                         >  1:875721/1‑234 (MQ=255)
                                                                                                                                           cATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCCGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGcc                                                                                                                                                   >  1:405097/1‑159 (MQ=255)
                                                                                                                                             tttttATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTgg                                                     >  1:794570/1‑251 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCAc                                                                                                                                                                                                                                              >  1:148949/1‑54 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATcc                                                                                                           >  1:488576/1‑185 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCaa                                                               <  1:543453/229‑1 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATa                                         >  1:367541/1‑251 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATa                                         >  1:838354/1‑251 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATa                                         >  1:819500/1‑251 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAGCGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGAGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGCTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACggg                                              >  1:739511/1‑246 (MQ=255)
                                                                                                                                                         aTTTATAAGAGGATATTTTATGAATGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCGTTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGGGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCtt                                                       <  1:491407/237‑1 (MQ=255)
                                                                                                                                                          tttATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGt                                                                                                                                    <  1:865814/159‑1 (MQ=255)
                                                                                                                                                          tttATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATcc                                                                                                           <  1:420325/184‑1 (MQ=255)
                                                                                                                                                                agagGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAAttt                                                                                                                                         >  1:445370/1‑148 (MQ=255)
                                                                                                                                                                 gagGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCttgttg                                 >  1:364507/1‑251 (MQ=255)
                                                                                                                                                                  agGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGt                                                                                                                                    <  1:309678/151‑1 (MQ=255)
                                                                                                                                                                  agGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGGAATACCttgttg                                 <  1:125272/250‑1 (MQ=255)
                                                                                                                                                                  agGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCACTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCttgttg                                 <  1:435145/250‑1 (MQ=255)
                                                                                                                                                                   ggATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGTGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAgctgc                                                                                                                                                                                        <  1:238716/98‑1 (MQ=255)
                                                                                                                                                                   ggATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGa                               >  1:290029/1‑251 (MQ=255)
                                                                                                                                                                   ggATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGa                               >  1:780639/1‑251 (MQ=255)
                                                                                                                                                                      tatTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAgg                            >  1:684654/1‑251 (MQ=255)
                                                                                                                                                                      tatTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTCAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATAc                                                          <  1:701958/221‑1 (MQ=255)
                                                                                                                                                                       atTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAggg                           <  1:782702/251‑1 (MQ=255)
                                                                                                                                                                       atTTTATGAGTGAAGTGTCAAATCTGGTAGAAAAATTCACTTGGGCGATTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAggg                           <  1:769419/251‑1 (MQ=255)
                                                                                                                                                                        ttttATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGt                                                                                                                                    >  1:538006/1‑145 (MQ=255)
                                                                                                                                                                        ttttATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATcc                                                                                                           >  1:307058/1‑170 (MQ=255)
                                                                                                                                                                        ttttATGAGTGAAGTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAGGGt                          >  1:75430/1‑251 (MQ=255)
                                                                                                                                                                                     gTTTCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTa                                                                                                              <  1:445383/154‑1 (MQ=255)
                                                                                                                                                                                      tttCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATgg                                                                                                                                                                                                  >  1:305312/1‑69 (MQ=255)
                                                                                                                                                                                      tttCAAATCTGGTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGc                                             >  1:720702/1‑218 (MQ=255)
                                                                                                                                                                                                 gTAGAAAAATTCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAGGGTATCGCCAAAAAGGGCGTTTCTAtt  <  1:287550/250‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
TGCCGGTTATTTATTTTTTATCTTAATATTATTTCCTGAAATTATTAAATAATCAGAGTTTTATTCTGCTATTTTCTCTGAGGTAGAATTTTATTTTCTTTACTCAATTGTTAAATTCTATTATAAATGGTTTCCATTGCATTTTTATGGTCAATTTATAAGAGGATATTTTATGAGTGAAGTTTCAAATCTGGTAGAAAAATCCACTTGGGCGAGTAAATTTTCTGCCTTAGGTCCGGGCATTGTCATGGCGTCAGCTGCTGTCGGTGGTTCGCATATTATCGCTTCAACCCAAGCCGGTGCAATTTATGGTTGGGAGTTGGTGAGCATTGTTATCCTGGTTAATTTATTCAAATATCCTTTTTTCCGATTTGGTGTGCAATATACCTTGGATACGGGCAATACCTTGTTGGAAGGGTATCGCCAAAAAGGGCGTTTCTATT  >  1/332510‑332952

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