Predicted mutation
evidence seq id position mutation annotation gene description
RA 4 8,616 C→T S3S (AGC→AGT)  unknown → Twin‑arginine translocation protein TatC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*48,6160CT98.8% 223.0 / NA 83S3S (AGC→AGT) unknownTwin‑arginine translocation protein TatC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (41/41);  minor base G (1/0);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCACAAAAGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGCGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGTTGCCGTTTTTATTTC  >  4/8374‑8858
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
aGCCACAAACGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACg                                                                                                                                                                                                                                             <  1:398918/250‑1 (MQ=255)
aGCCACAAACGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACg                                                                                                                                                                                                                                             <  1:646530/250‑1 (MQ=255)
aGCCACAAACGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACg                                                                                                                                                                                                                                             <  1:14059/250‑1 (MQ=255)
aGCCACAAACGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACg                                                                                                                                                                                                                                             <  1:687691/250‑1 (MQ=255)
  ccACAAACGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGAt                                                                                                                                                                                                                                           <  1:15992/250‑1 (MQ=255)
        acggAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAAc                                                                                                                                                                                                                                    <  1:482905/249‑1 (MQ=255)
         cggAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAAcc                                                                                                                                                                                                                                   <  1:487815/250‑1 (MQ=255)
            aaGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAAcccc                                                                                                                                                                                                                                 <  1:899877/250‑1 (MQ=255)
                         cacaAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCAtcact                                                                                                                                                                                                                         >  1:649953/1‑245 (MQ=255)
                         cacaAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCAtcac                                                                                                                                                                                                                          >  1:471501/1‑244 (MQ=255)
                              aaaCAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGa                                                                                                                                                                                                              <  1:319678/251‑1 (MQ=255)
                                acaacGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTACGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAAc                                                                                                                                                                                                            <  1:659294/251‑1 (MQ=255)
                                   acGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTa                                                                                                                                                                                                          >  1:643909/1‑250 (MQ=255)
                                      gAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGc                                                                                                                                                                                                       <  1:702054/250‑1 (MQ=255)
                                      gAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCa                                                                                                                                                                                                      >  1:267232/1‑251 (MQ=255)
                                        aaaaaCGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCa                                                                                                                                                                                                      >  1:231475/1‑249 (MQ=255)
                                        aaaaaCGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCa                                                                                                                                                                                                      >  1:563118/1‑249 (MQ=255)
                                        aaaaaCGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACt                                                                                                                                                                                                    <  1:149323/251‑1 (MQ=255)
                                                 cGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCttatta                                                                                                                                                                                           <  1:568083/251‑1 (MQ=255)
                                                   cACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATTCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTAc                                                                                                                                                                                          <  1:455169/250‑1 (MQ=255)
                                                       ttCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAATTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGt                                                                                                                                                                                                                 <  1:711998/223‑1 (MQ=255)
                                                       ttCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCAtcactcac                                                                                                                                                                                                                      >  1:76600/1‑218 (MQ=255)
                                                       ttCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCttatta                                                                                                                                                                                           >  1:716152/1‑245 (MQ=255)
                                                            cggcggAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCattat                                                                                                                                                                                >  1:47120/1‑251 (MQ=255)
                                                              gcggAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCattatt                                                                                                                                                                               <  1:548515/250‑1 (MQ=255)
                                                               cggAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTAtttt                                                                                                                                                                             <  1:274163/251‑1 (MQ=255)
                                                                      ggCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTAttt                                                                                                                                                                              >  1:93457/1‑243 (MQ=255)
                                                                       gCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCg                                                                                                                                                                     >  1:108578/1‑251 (MQ=255)
                                                                                 aagaaATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGtgtttgt                                                                                                                                                           <  1:461149/251‑1 (MQ=255)
                                                                                   gaaATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTttgttg                                                                                                                                                         >  1:789088/1‑251 (MQ=255)
                                                                                    aaaTCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGtgtttgtt                                                                                                                                                          <  1:78639/249‑1 (MQ=255)
                                                                                           aCTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGtt                                                                                                                                                 <  1:896957/251‑1 (MQ=255)
                                                                                            cTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGttt                                                                                                                                                <  1:642256/251‑1 (MQ=255)
                                                                                                     aaaCTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTa                                                                                                                                       >  1:526224/1‑251 (MQ=255)
                                                                                                                ttACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACgc                                                                                                                                                                                        >  1:852908/1‑191 (MQ=255)
                                                                                                                   cATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAAt                                                                                                                         <  1:97701/251‑1 (MQ=255)
                                                                                                                       aatcaatACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAAc                                                                                                                                                                                                            <  1:647994/164‑1 (MQ=255)
                                                                                                                       aatcaatACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTAc                                                                                                                                                                                          >  1:20537/1‑182 (MQ=255)
                                                                                                                       aatcaatACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTtcgtc                                                                                                                                                                         <  1:378944/199‑1 (MQ=255)
                                                                                                                       aatcaatACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTtcgtc                                                                                                                                                                         <  1:640397/199‑1 (MQ=255)
                                                                                                                                gATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTAttt                                                                                                                                                                              >  1:206698/1‑185 (MQ=255)
                                                                                                                                gATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGa                                                                                                                                   >  1:686587/1‑228 (MQ=255)
                                                                                                                                gATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCa                                                                                                             >  1:13621/1‑250 (MQ=255)
                                                                                                                                 aTCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTa                                                                                                                                                                                                                >  1:250074/1‑150 (MQ=255)
                                                                                                                                     ggACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCg                                                                                                          >  1:411682/1‑248 (MQ=255)
                                                                                                                                                    ggCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGc                                                                                                  >  1:439439/1‑241 (MQ=255)
                                                                                                                                                              gCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTA‑TTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAg                                                                             >  1:505502/1‑251 (MQ=255)
                                                                                                                                                                 tCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           <  1:241837/251‑1 (MQ=255)
                                                                                                                                                                 tCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           <  1:707824/251‑1 (MQ=255)
                                                                                                                                                                 tCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           <  1:33348/251‑1 (MQ=255)
                                                                                                                                                                 tCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           <  1:657860/251‑1 (MQ=255)
                                                                                                                                                                    gCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCgtggtg                                                                                                                                                                >  1:488642/1‑163 (MQ=255)
                                                                                                                                                                       caTTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACtt                                                                                                                                            >  1:150481/1‑180 (MQ=255)
                                                                                                                                                                       caTTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTc                                                                                >  1:173894/1‑240 (MQ=255)
                                                                                                                                                                           tGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           >  1:340278/1‑241 (MQ=255)
                                                                                                                                                                                   gTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATc                                                          >  1:801748/1‑250 (MQ=255)
                                                                                                                                                                                                    aCCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           <  1:336985/216‑1 (MQ=255)
                                                                                                                                                                                                     ccGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTAccc                                       <  1:58668/251‑1 (MQ=255)
                                                                                                                                                                                                            ttGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAg                                                                                                                                                   <  1:88463/136‑1 (MQ=255)
                                                                                                                                                                                                            ttGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAAc                                                                                     <  1:410889/198‑1 (MQ=255)
                                                                                                                                                                                                            ttGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           >  1:510255/1‑208 (MQ=255)
                                                                                                                                                                                                            ttGCGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTaa                                >  1:433506/1‑251 (MQ=255)
                                                                                                                                                                                                               cGGGATAGCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCTACCAACATCGCCACGCCGTTTTTTACCCCGATTAAAc                              <  1:876819/250‑1 (MQ=255)
                                                                                                                                                                                                                      gCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCTTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGc                                                        <  1:453355/217‑1 (MQ=255)
                                                                                                                                                                                                                      gCCTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           >  1:84301/1‑198 (MQ=255)
                                                                                                                                                                                                                       ccTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACt                                                                                                                                                     <  1:107656/123‑1 (MQ=255)
                                                                                                                                                                                                                       ccTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGCCAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGgtg                     >  1:437485/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGCCAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGgtg                     >  1:432058/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccTCACATCAGGATAACTTATGACCAGGGTTGACGATACTCAACCCCTCATCACTCACTTGGAAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGGGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCCCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAAAACCGCCACGCCGTTTTTTACCCCGATTAAAcctgcctctt                     >  1:757649/1‑242 (MQ=255)
                                                                                                                                                                                                                        cTCACATCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGgtgt                    <  1:465351/251‑1 (MQ=255)
                                                                                                                                                                                                                          cacaTCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGAt                                                                       >  1:406873/1‑198 (MQ=255)
                                                                                                                                                                                                                            caTCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTg                                                                                                       <  1:48631/164‑1 (MQ=255)
                                                                                                                                                                                                                            caTCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGtt                >  1:833290/1‑251 (MQ=255)
                                                                                                                                                                                                                            caTCAGGATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGtt                >  1:152657/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggATAACTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           >  1:147152/1‑187 (MQ=255)
                                                                                                                                                                                                                                      aCTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGTTGCCGTTTTTa      >  1:521108/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCttatta                                                                                                                                                                                           <  1:758972/69‑1 (MQ=255)
                                                                                                                                                                                                                                       cTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCa                                                                           >  1:783418/1‑181 (MQ=255)
                                                                                                                                                                                                                                       cTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAAc                              >  1:17515/1‑226 (MQ=255)
                                                                                                                                                                                                                                       cTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGTTGCCGTTTTTAt     >  1:741883/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cTTATGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGTTGCCGTTTTTAt     >  1:858482/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGTTGCCGTTTTTATTTc  <  1:896042/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tGACCAGTGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGc                                                                                                                                                                                                 <  1:472446/59‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
AGCCACAAAAGGAAGCCCTCGCTGCCACAAAAACAACGGAAAAAACGACCGCACTTTCCCCGGCGGAGCTGGCGGAACAGGAAGAAATCGAACTGGACGAAAAACTCGCCCATTACATCAATCAATACGATCCGGACAATCCGCTGCCGGCTAATACGGCATCGGCACATTTGGCCAAAGTCTCAGAGAAAAGTTAACCGTGCTTTGCGGGATAGCCTCACATCAGGATAACTTATGACCAGCGTTGACGATACTCAACCCCTCATCACTCACTTGGTAGAACTACGCACTCGCTTATTACGCAGCATTATTTTCGTCGCCGTGGTGTTTGTTGCACTAGTTTACTTTTCTAACGAAATTTATAATTTTGTGGCGGCACCGCTGACAGCGAACCTGCCGAACGGTTCCAGCATGATTGCCACCAACATCGCCACGCCGTTTTTTACCCCGATTAAACTCACCGGTGTGGTTGCCGTTTTTATTTC  >  4/8374‑8858

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