Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,151,540 G→A E473K (GAA→AAA)  unknown → Autoinducer 2 (AI‑2) kinase LsrK (EC 2.7.1.‑)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,151,5400GA90.3% 134.7 / ‑4.4 62E473K (GAA→AAA) unknownAutoinducer 2 (AI‑2) kinase LsrK (EC 2.7.1.‑)
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (29/27);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAAAATGCTGCTATCGTTTCTTCCTGCAATTTGGATCAGGTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATTCGACCTTAGCCGAGGTGGGGGAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATACTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTGGAAGGCACCAGGGTTGTAGCGAACGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACATCGGTAAAACCC  >  1/1151297‑1151777
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
gaagaaAATGCTGCTATCGTTTCTTCCTGCAATTTGGATCAGGTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAAc                                                                                                                                                                                                                                         <  1:751800/250‑1 (MQ=255)
                    ttcttcCTGCAATTTGGATCAGGTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAAtt                                                                                                                                                                                                                               <  1:482650/240‑1 (MQ=255)
                               aTTTGGATCAGGTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAAt                                                                                                                                                                                                         <  1:363283/251‑1 (MQ=255)
                               aTTTGGATCAGGTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAAATGGGTAAATTTGAACGGCAACATCAACCGAAt                                                                                                                                                                                                         >  1:598396/1‑251 (MQ=255)
                                          gTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGGAAAAACTGGTTAAATTTTAACGGCAACATCAACCGAATCCCGATAACTa                                                                                                                                                                                              >  1:872627/1‑251 (MQ=255)
                                          gTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTa                                                                                                                                                                                              >  1:439114/1‑251 (MQ=255)
                                             cAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGc                                                                                                                                                                                           <  1:870574/251‑1 (MQ=255)
                                                  aTTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTtat                                                                                                                                                                                      >  1:584872/1‑251 (MQ=255)
                                                    tctcAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGGCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATAc                                                                                                                                                                                    <  1:442103/251‑1 (MQ=255)
                                                               aaaTTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACaa                                                                                                                                                                         <  1:859423/251‑1 (MQ=255)
                                                                 aTTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAg                                                                                                                                                                       >  1:774717/1‑251 (MQ=255)
                                                                    aaCCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGcac                                                                                                                                                                    <  1:348556/251‑1 (MQ=255)
                                                                         ccATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGcac                                                                                                                                                                    <  1:423088/246‑1 (MQ=255)
                                                                         ccATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATgg                                                                                                                                                               <  1:195019/251‑1 (MQ=255)
                                                                           aTTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGc                                                                                                                                                             <  1:824343/251‑1 (MQ=255)
                                                                              gTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAAtt                                                                                                                                                                                                                               <  1:158142/182‑1 (MQ=255)
                                                                              gTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAAcc                                                                                                                                                                                                             <  1:407676/200‑1 (MQ=255)
                                                                                cTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATa                                                                                                                                                        <  1:864446/251‑1 (MQ=255)
                                                                                               gtttAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGGAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAAAAAt                                                                                                                                         >  1:871066/4‑251 (MQ=255)
                                                                                               gTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCTCAATGGGCAGATATTTATAAAAAACAAt                                                                                                                                         >  1:216211/1‑251 (MQ=255)
                                                                                                tCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGGAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCa                                                                                                                                                                                  >  1:400175/1‑209 (MQ=255)
                                                                                                      ggCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATAc                                                                                                                                                                                    <  1:728920/201‑1 (MQ=255)
                                                                                                       gCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAAc                                                                                                                                   <  1:245507/249‑1 (MQ=255)
                                                                                                        cAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGGAAAAACTGGGTAAATTTGAACGGCAACATCAAACGAATTCCGATAACTAAGGTTTTTACCAAGCTCAAAAAGCAAAATGGGCAGAAATTTATAAAAAACAATTACAACTgg                                                                                                                                >  1:78209/1‑251 (MQ=255)
                                                                                                                 gAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGGAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAAAAATTACAACTGGTTGATAGCg                                                                                                                       >  1:201116/1‑251 (MQ=255)
                                                                                                                  aGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATgg                                                                                                                                                               <  1:198405/210‑1 (MQ=255)
                                                                                                                       gATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATGACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTg                                                                                                                                 >  1:429467/1‑235 (MQ=255)
                                                                                                                       gATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTaa                                                                                                                 >  1:775665/1‑251 (MQ=255)
                                                                                                                        aTCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGGAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAAc                                                                                                                                                                                                >  1:631501/1‑171 (MQ=255)
                                                                                                                        aTCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAATTAGTTAAATTTGAACGGCAACATGAACCCAATCCCGATAACAATGCTTTATACCCAGCTCACAAAGCACAAAGGGCAGATATTTATAAAAAACAATTAAAACTGGATGAGAGGGGATTaa                                                                                                                 >  1:784371/1‑250 (MQ=255)
                                                                                                                                 gATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGGTGATAGCGGATTAACTACTTCttt                                                                                                       >  1:99631/1‑251 (MQ=255)
                                                                                                                                           gATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAGTTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGc                                                                                              <  1:357319/250‑1 (MQ=255)
                                                                                                                                           gATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTTAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATTCCCAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCa                                                                                             >  1:23436/1‑251 (MQ=255)
                                                                                                                                           gATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATAc                                                                                                                                                                                    >  1:90446/1‑164 (MQ=255)
                                                                                                                                                aCGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGcac                                                                                                                                                                    <  1:346797/175‑1 (MQ=255)
                                                                                                                                                aCGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTtat                                                                                                                                                   >  1:677111/1‑192 (MQ=255)
                                                                                                                                                aCGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAAATGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAAAGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAgg                                                                                        >  1:599234/1‑251 (MQ=255)
                                                                                                                                                       aCATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAg                                                                                 <  1:712643/251‑1 (MQ=255)
                                                                                                                                                        cATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATAc                                                                                                                                                                                    >  1:91793/1‑151 (MQ=255)
                                                                                                                                                                 tGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATaa                                                                       >  1:678313/1‑251 (MQ=255)
                                                                                                                                                                                 gcgcTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCCACATCAACCGAATCCCGATAACTATGc                                                                                                                                                                                           >  1:823838/1‑119 (MQ=255)
                                                                                                                                                                                 gcgcTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAAc                                                                                                                                                                                                              <  1:652220/100‑1 (MQ=255)
                                                                                                                                                                                 gcgcTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGc                                                                                                                                                                                           >  1:820824/1‑119 (MQ=255)
                                                                                                                                                                                 gcgcTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCa                                                                                                                                                            <  1:135464/150‑1 (MQ=255)
                                                                                                                                                                                 gcgcTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTtgtag                                                                                                   >  1:125175/1‑205 (MQ=255)
                                                                                                                                                                                 gcgcTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAgggggg                                                       >  1:831482/1‑251 (MQ=255)
                                                                                                                                                                                      gggTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTtat                                                                                                                                                   >  1:280073/1‑154 (MQ=255)
                                                                                                                                                                                                  tGCCGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAg                                                                                                                         <  1:250204/168‑1 (MQ=255)
                                                                                                                                                                                                    ccGCAGGCGTCGGAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAAACCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACAGGTTGATAGCGGATTAACTACTTCTTTGTAGAAAGCACCAGGGGCGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGAttttt                                    >  1:336902/1‑251 (MQ=255)
                                                                                                                                                                                                                gAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGTGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATAc                         <  1:253458/250‑1 (MQ=255)
                                                                                                                                                                                                                gAGTCGGCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATAc                         <  1:292046/250‑1 (MQ=255)
                                                                                                                                                                                                                      gCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTaaaa                  >  1:194914/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCATTTATCCGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAACTCACAAAGCACAAAGTGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTATGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTaaaa                  >  1:551592/1‑251 (MQ=255)
                                                                                                                                                                                                                           tatccGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAacta                                                                                                              >  1:125405/5‑154 (MQ=255)
                                                                                                                                                                                                                           tatccGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACAt             >  1:33208/5‑251 (MQ=255)
                                                                                                                                                                                                                            atccGACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACAt             >  1:113542/4‑250 (MQ=255)
                                                                                                                                                                                                                                gACCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTtgtag                                                                                                   >  1:410807/1‑158 (MQ=255)
                                                                                                                                                                                                                                 aCCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGc                                                                                                                                                                               >  1:609791/1‑83 (MQ=255)
                                                                                                                                                                                                                                 aCCTTAGCCGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACATCGGTaa       >  1:538333/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACATCGGTAAAAccc  <  1:303988/249‑1 (MQ=255)
                                                                                                                                                                                                                                        ccGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACATCGGTAAAAccc  <  1:822620/249‑1 (MQ=255)
                                                                                                                                                                                                                                        ccGAGGCGGGGAAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATATTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTAGAAGGCACCAGGGTTGTAGCGAATGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACATCGGTAAAAccc  <  1:548765/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
GAAGAAAATGCTGCTATCGTTTCTTCCTGCAATTTGGATCAGGTGCAAAAATTCTCAGGAGTTAAATTAACCTCCATTGTCTTTGCCGGCGGCGGGTCTAAAGGCAAATTATGGAGCCAGATCTTAGCCGATGTCACCGGATTAACGGTAAACATTCCGATTGTCAAAGAAGCTACAGCGCTGGGTTGTGCCATTGCCGCAGGCGTCGGAGTCGGCATTTATTCGACCTTAGCCGAGGTGGGGGAAAAACTGGTTAAATTTGAACGGCAACATCAACCGAATCCCGATAACTATGCTTTATACCAAGCTCACAAAGCACAATGGGCAGATACTTATAAAAAACAATTACAACTGGTTGATAGCGGATTAACTACTTCTTTGTGGAAGGCACCAGGGTTGTAGCGAACGATAAAAAAATCGTAGGGGGGACAACTCTCTACGATTTTTTGTTGATATACCGTAAAAAACATCGGTAAAACCC  >  1/1151297‑1151777

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