Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,569,293 T→C E580E (GAA→GAG)  unknown ← Phosphoglycerate transport system sensor protein PgtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,569,2930TC100.0% 255.4 / NA 84E580E (GAA→GAG) unknownPhosphoglycerate transport system sensor protein PgtB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/40);  total (44/40)

TAATTAAACATTATCTACCACCTTAAATTCCAAAATAACCAGCGCATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATTTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATATTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGTATTTGTGGTTTGGCGGGAAAAATGGCGCAGTCTTTTCAC  >  1/1569047‑1569535
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
tAATTAAACATTATCTACCACCTTAAATTCCAAAATAACCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCaa                                                                                                                                                                                                                                                >  1:437118/1‑251 (MQ=255)
        cATTATCTACCACCTTAAATTCCAAAATAACCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGt                                                                                                                                                                                                                                        >  1:308608/1‑251 (MQ=255)
            aTCTACCACCTTAAATTCCAAAATAACCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGttttt                                                                                                                                                                                                                                    >  1:470412/1‑251 (MQ=255)
                         aaTTCCAAAATAACCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGc                                                                                                                                                                                                                       <  1:811138/251‑1 (MQ=255)
                           ttCCAAAATAACCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTCATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCtt                                                                                                                                                                                                                     <  1:423417/251‑1 (MQ=255)
                                 aaTAACCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTg                                                                                                                                                                                                               <  1:240723/251‑1 (MQ=255)
                                    aaCCAGCGTATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCtt                                                                                                                                                                                                                     >  1:841573/1‑242 (MQ=255)
                                                        aTGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAat                                                                                                                                                                                        >  1:22927/1‑251 (MQ=255)
                                                                           gTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAAttt                                                                                                                                                                     >  1:790703/1‑251 (MQ=255)
                                                                                        tgttgCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        <  1:329450/251‑1 (MQ=255)
                                                                                         gttgCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAatactt                                                                                                                                                                                    >  1:168341/1‑219 (MQ=255)
                                                                                            gCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAGACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCtt                                                                                                                                                                                                                     >  1:218880/1‑186 (MQ=255)
                                                                                              aTAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAgg                                                                                                                                                  >  1:537109/1‑251 (MQ=255)
                                                                                                  tGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTaaa                                                                                                                                              >  1:597795/1‑251 (MQ=255)
                                                                                                   ggATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGaa                                                                                                                                                                                                     >  1:501984/1‑195 (MQ=255)
                                                                                                   ggATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGc                                                                                                                                                                                  <  1:640240/214‑1 (MQ=255)
                                                                                                   ggATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAAt                                                                                                                                             >  1:610717/1‑251 (MQ=255)
                                                                                                     aTAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATc                                                                                                                                            <  1:549327/250‑1 (MQ=255)
                                                                                                     aTAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATc                                                                                                                                            <  1:773502/250‑1 (MQ=255)
                                                                                                                gATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        >  1:167030/1‑227 (MQ=255)
                                                                                                                gATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAgaga                                                                                                                                >  1:467098/1‑251 (MQ=255)
                                                                                                                 aTAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        >  1:820915/1‑226 (MQ=255)
                                                                                                                    aTCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATaa                                                                                                                             >  1:240768/1‑250 (MQ=255)
                                                                                                                    aTCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATaa                                                                                                                             >  1:761726/1‑250 (MQ=255)
                                                                                                                     tCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAAc                                                                                                                            >  1:228256/1‑250 (MQ=255)
                                                                                                                           gCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        <  1:517741/216‑1 (MQ=255)
                                                                                                                           gCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATc                                                                                                                                            >  1:493284/1‑228 (MQ=255)
                                                                                                                                  ttAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCg                                                                                                                          <  1:183987/239‑1 (MQ=255)
                                                                                                                                     aTGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATAttt                                                                                                           <  1:2283/251‑1 (MQ=255)
                                                                                                                                       ggATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        <  1:614342/204‑1 (MQ=255)
                                                                                                                                        gATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAAc                                                                                                                            <  1:831347/231‑1 (MQ=255)
                                                                                                                                        gATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTAt                                                                                                         <  1:301004/250‑1 (MQ=255)
                                                                                                                                        gATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATg                                                                                                        >  1:203167/1‑251 (MQ=255)
                                                                                                                                         aTTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTcc                                                                                                                                                                                                                   <  1:604628/143‑1 (MQ=255)
                                                                                                                                         aTTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGt                                                                                                       >  1:294489/1‑251 (MQ=255)
                                                                                                                                         aTTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTAGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTc                                                                                                                                                                                                                    >  1:485255/1‑142 (MQ=255)
                                                                                                                                              ctactGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTgg                                                                                                  >  1:871769/1‑251 (MQ=255)
                                                                                                                                                 ctGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCt                                                                                                >  1:569344/1‑250 (MQ=255)
                                                                                                                                                   ggAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGaa                                                                                                                                                                                                     <  1:497784/147‑1 (MQ=255)
                                                                                                                                                   ggAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAATGAGTTCAGGaa                                                                                                                                                                                                     <  1:654688/147‑1 (MQ=255)
                                                                                                                                                    gagaAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        <  1:267217/191‑1 (MQ=255)
                                                                                                                                                         aGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAgg                                                                                                                                      >  1:683399/1‑204 (MQ=255)
                                                                                                                                                            ttGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATcc                                                                                     <  1:294710/250‑1 (MQ=255)
                                                                                                                                                              gtagtaATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGAt                                                                                                                               <  1:656177/206‑1 (MQ=255)
                                                                                                                                                              gtagtaATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCGCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGt                                                                                                       <  1:433029/230‑1 (MQ=255)
                                                                                                                                                                    aTTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCatta                                                                            <  1:171198/251‑1 (MQ=255)
                                                                                                                                                                     tttATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGt                                                                                                       >  1:670505/1‑223 (MQ=255)
                                                                                                                                                                        aTCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCAttattta                                                                        <  1:361926/251‑1 (MQ=255)
                                                                                                                                                                           tGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATTTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATac                                                                     <  1:449364/251‑1 (MQ=255)
                                                                                                                                                                           tGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATac                                                                     <  1:347197/251‑1 (MQ=255)
                                                                                                                                                                            gACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATac                                                                     <  1:575646/250‑1 (MQ=255)
                                                                                                                                                                             aCAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATaca                                                                    <  1:472408/250‑1 (MQ=255)
                                                                                                                                                                                 cGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGctc                                                                                                                                                                                >  1:23656/1‑138 (MQ=255)
                                                                                                                                                                                  ggCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCg                                                                                                                          <  1:410449/191‑1 (MQ=255)
                                                                                                                                                                                        ccGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATcc                                                                                     >  1:439908/1‑222 (MQ=255)
                                                                                                                                                                                              ccATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATac                                                                     >  1:673854/1‑232 (MQ=255)
                                                                                                                                                                                               cATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAg                                                                                                                                       <  1:8904/165‑1 (MQ=255)
                                                                                                                                                                                                        aaaCCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGt                                         <  1:874735/250‑1 (MQ=255)
                                                                                                                                                                                                          aCCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGat                                                                                                               >  1:516597/1‑178 (MQ=255)
                                                                                                                                                                                                           cccATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTTTCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAAt                                                                                                                                                                       >  1:41809/1‑121 (MQ=255)
                                                                                                                                                                                                           cccATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        <  1:6275/136‑1 (MQ=255)
                                                                                                                                                                                                           cccATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATcc                                                                                     >  1:827456/1‑203 (MQ=255)
                                                                                                                                                                                                           cccATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGc                                                                                <  1:68853/208‑1 (MQ=255)
                                                                                                                                                                                                              aTAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATGCATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGTATTtgtg                                  <  1:143368/251‑1 (MQ=255)
                                                                                                                                                                                                               tAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATcc                                                                                     >  1:662852/1‑199 (MQ=255)
                                                                                                                                                                                                                 aTTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATc                                                                                                                                            >  1:883436/1‑142 (MQ=255)
                                                                                                                                                                                                                 aTTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGa                                                                                                                                     >  1:602039/1‑149 (MQ=255)
                                                                                                                                                                                                                 aTTCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCt                                                                                                <  1:639797/186‑1 (MQ=255)
                                                                                                                                                                                                                  ttCTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACa                                                                                                                                                       <  1:718846/130‑1 (MQ=255)
                                                                                                                                                                                                                    cTACTCTCTGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATcc                                                                                     >  1:339478/1‑194 (MQ=255)
                                                                                                                                                                                                                        tctctGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        >  1:754744/1‑123 (MQ=255)
                                                                                                                                                                                                                         ctctGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCg                                                                                                                          >  1:709355/1‑152 (MQ=255)
                                                                                                                                                                                                                           ctGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAaattaa                                                      <  1:829073/218‑1 (MQ=255)
                                                                                                                                                                                                                           ctGATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGTATTTGTGGTTTGGCGGGaaa                     <  1:705316/251‑1 (MQ=255)
                                                                                                                                                                                                                             gATCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATAc                                                                                                                                                        <  1:561642/118‑1 (MQ=255)
                                                                                                                                                                                                                              aTCCGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATac                                                                     >  1:212646/1‑200 (MQ=255)
                                                                                                                                                                                                                                ccggtatatTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCtt                                                                                                                                                                                                                     <  1:830932/50‑1 (MQ=255)
                                                                                                                                                                                                                                ccGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGc                                                                                                                                                                                                                       >  1:864765/1‑52 (MQ=255)
                                                                                                                                                                                                                                ccGATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATc                                                                                                                                                     >  1:359372/1‑118 (MQ=255)
                                                                                                                                                                                                                                  gATATATTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGGGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATc                                                                                                                                            >  1:749836/1‑125 (MQ=255)
                                                                                                                                                                                                                                   atatatTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCg                                                                                                                          >  1:3723/1‑142 (MQ=255)
                                                                                                                                                                                                                                   atatatTTTGTATTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGTATTTGTGGTTTGGCGGGAAAAATGgcg              <  1:754381/250‑1 (MQ=255)
                                                                                                                                                                                                                                              aTTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCAttattta                                                                        >  1:154322/1‑181 (MQ=255)
                                                                                                                                                                                                                                              aTTTTTATCTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATACTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGTATTTGTGGTTTGGCGGGAAAAATGGCGCAGTCTTTTcac  <  1:765206/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
TAATTAAACATTATCTACCACCTTAAATTCCAAAATAACCAGCGCATTTTTTGTCAATGTAGAAGCAATTAATAAGTTTCCTTGACATTGTTGCATAATGGATAGACTTATTGATAATCCTAAGCCTAAATTAATGGATTTACTACTGGAGAAAGGTTGTAGTAATTTATCTGACAACGGCCAACCGGAACCATTATCAGAAACCCATAATTCTACTCTCTGATCCGATATATTTTGTATTTTTATTTCAATCTCCGGTTTTTCATGCTGCATAGCTTCCAGTGAGTTCAGGAAAATATTCACAAATATTTGCTCTAATAAAATTTCATCACCTAATACATCAGGTAAATCTGTAGGAAGAGATAACCGAGCATTCAGATATTTATGTTTTGGCTCAAGTAAATCCCACGCATTATTTATACACATCAATAAATTAATTTTCTGTAATGTATTTGTGGTTTGGCGGGAAAAATGGCGCAGTCTTTTCAC  >  1/1569047‑1569535

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