Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 824,112 T→C P84P (CCA→CCG)  unknown ← 5'‑deoxynucleotidase YfbR (EC 3.1.3.89)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1824,1120TC98.5% 185.4 / ‑5.4 67P84P (CCA→CCG) unknown5'‑deoxynucleotidase YfbR (EC 3.1.3.89)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (44/22);  minor base G (1/0);  total (45/22)
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

CACTTTTGCCGATCTAAATTCCTGATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATTGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCGTGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGATCCAAACAGGCAAAAAAATGGCTGGTTTTTATT  >  1/823872‑824355
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
cACTTTTGCCGATCTAAATTCCTGATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGc                                                                                                                                                                                                                                            <  1:171305/250‑1 (MQ=255)
            tCTAAATTCCTGATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTaa                                                                                                                                                                                                                               >  1:409910/1‑251 (MQ=255)
                 aTTCCTGATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAAtt                                                                                                                                                                                                                          >  1:159268/1‑251 (MQ=255)
                    ccTGATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCCAAATATTTTTTTCATCATCCGAAAATTAGGGCCTATCCAGTTAGGGTAAAAAATCTCCCGGCAATTCTTCCGGCAACAACCTAATCAAAT‑CAACCTCGCCGCCGGCCTCAATTGGTTTGTAAGCATGGGAAATCTCTTCATTGAATTATTTAATGGGGGTCGGCAGAT‑GGCGGTaaaaaatttc                                                                                                                                                                                                                       >  1:482940/1‑245 (MQ=255)
                       gATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGAt                                                                                                                                                                                                                    >  1:857923/1‑251 (MQ=255)
                       gATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGAt                                                                                                                                                                                                                    >  1:547772/1‑251 (MQ=255)
                         ttACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCCCAGATTAAATTTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAGTATTTGATCGGTGTCGGCAgttc                                                                                                                                                                                                                                       >  1:706774/1‑227 (MQ=255)
                            ccATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATgtgt                                                                                                                                                                                                                >  1:120236/1‑250 (MQ=255)
                                 tCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGAt                                                                                                                                                                                                          >  1:514734/1‑251 (MQ=255)
                                                      gATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATg                                                                                                                                                                                      >  1:857858/1‑250 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATc                                                                                                                                                                                               <  1:650567/229‑1 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATgcgc                                                                                                                                                                                   >  1:151444/1‑241 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCa                                                                                                                                                                         >  1:728786/1‑251 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCa                                                                                                                                                                         >  1:411816/1‑251 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCa                                                                                                                                                                         >  1:36376/1‑251 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCa                                                                                                                                                                         >  1:276243/1‑251 (MQ=255)
                                                                  gATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCa                                                                                                                                                                         >  1:759361/1‑251 (MQ=255)
                                                                    ttAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGt                                                                                                                                                                       <  1:200770/251‑1 (MQ=255)
                                                                                gTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCAt                                                                                                                                                                            <  1:42650/234‑1 (MQ=255)
                                                                                 tttCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTaaa                                                                                                                                                          >  1:11194/1‑251 (MQ=255)
                                                                                    cacaATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTaaaaa                                                                                                                                                        >  1:558733/1‑250 (MQ=255)
                                                                                        atatGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGAtt                                                                                                                                                   >  1:770853/1‑251 (MQ=255)
                                                                                          atGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATa                                                                                                                                                                                                         >  1:260211/1‑195 (MQ=255)
                                                                                          atGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGAtttt                                                                                                                                                 >  1:457901/1‑251 (MQ=255)
                                                                                                                   tGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGc                                                                                                                                                                              >  1:890162/1‑197 (MQ=255)
                                                                                                                      cTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAg                                                                                                                                                                                                                                          >  1:799049/1‑134 (MQ=255)
                                                                                                                      cTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACaaa                                                                                                                     >  1:99399/1‑251 (MQ=255)
                                                                                                                       tATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAAt                                                                                                                    >  1:203941/1‑251 (MQ=255)
                                                                                                                        aTCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGgtgt                                                                                                                                                                                                                                                >  1:557730/1‑126 (MQ=255)
                                                                                                                        aTCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCAc                                                                                                                                                                >  1:140962/1‑206 (MQ=255)
                                                                                                                        aTCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATg                                                                                                                   >  1:290522/1‑251 (MQ=255)
                                                                                                                         tCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATg                                                                                                                   >  1:121933/1‑250 (MQ=255)
                                                                                                                                 gggTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATtaaggctaaggc                                                                                                                                 <  1:312771/228‑1 (MQ=255)
                                                                                                                                 gggTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATtaaggctaagg                                                                                                                                  >  1:382653/1‑227 (MQ=255)
                                                                                                                                  ggTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTaa                                                                                                         >  1:217756/1‑251 (MQ=255)
                                                                                                                                   gTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATa                                                                                                                                                                                                         >  1:453079/1‑154 (MQ=255)
                                                                                                                                   gTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGc                                                                                                                                                                                                   >  1:841507/1‑160 (MQ=255)
                                                                                                                                    tAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CCCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGc                                                                                                       <  1:208659/251‑1 (MQ=255)
                                                                                                                                         aaTCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGt                                                                                                                                                                       >  1:660098/1‑182 (MQ=255)
                                                                                                                                                        cttcCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTcc                                                                                   <  1:743713/251‑1 (MQ=255)
                                                                                                                                                              gCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATtaaggctaagg                                                                                                                                  <  1:389165/198‑1 (MQ=255)
                                                                                                                                                                aacaaACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAat                                                                           <  1:887786/251‑1 (MQ=255)
                                                                                                                                                                   aaaCTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATAtt                                                                        >  1:797934/1‑251 (MQ=255)
                                                                                                                                                                      cTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAg                                                                                                                                                                                                                                          >  1:647576/1‑86 (MQ=255)
                                                                                                                                                                       tGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGc                                                                    <  1:206415/251‑1 (MQ=255)
                                                                                                                                                                       tGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGc                                                                    <  1:119805/251‑1 (MQ=255)
                                                                                                                                                                       tGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGc                                                                    <  1:557678/251‑1 (MQ=255)
                                                                                                                                                                        gATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGc                                                                                                                                                                  >  1:256381/1‑156 (MQ=255)
                                                                                                                                                                        gATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTAc                                                                      >  1:203403/1‑248 (MQ=255)
                                                                                                                                                                        gATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCa                                                                   >  1:339632/1‑251 (MQ=255)
                                                                                                                                                                            aaaTGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAg                                                               <  1:596209/251‑1 (MQ=255)
                                                                                                                                                                               tGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGa                                                            <  1:736748/251‑1 (MQ=255)
                                                                                                                                                                                               gTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGc                                                                    <  1:672757/228‑1 (MQ=255)
                                                                                                                                                                                               gTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCCGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGaatcaa                                             >  1:347696/1‑251 (MQ=255)
                                                                                                                                                                                                 ctcAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGaatcaat                                            >  1:688010/1‑250 (MQ=255)
                                                                                                                                                                                                     aTTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATtaaggctaaggc                                                                                                                                 <  1:616694/161‑1 (MQ=255)
                                                                                                                                                                                                      tttgtttgTAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAg                                      <  1:491257/251‑1 (MQ=255)
                                                                                                                                                                                                         gtttgtAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCGATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGa                                   >  1:351251/1‑251 (MQ=255)
                                                                                                                                                                                                           ttgtAAGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCAt                                                                  >  1:797913/1‑218 (MQ=255)
                                                                                                                                                                                                                aGCATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGaatcaatc                                           >  1:709068/1‑236 (MQ=255)
                                                                                                                                                                                                                  cATGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATTGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGATCCAAAc                            <  1:436838/249‑1 (MQ=255)
                                                                                                                                                                                                                    tGGGTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTtg                                                                                                                                                                     <  1:187386/110‑1 (MQ=255)
                                                                                                                                                                                                                      ggTAATCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGATCCAAACAGGCaa                      >  1:405206/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaTCTCTTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGc                                                                    >  1:665358/1‑202 (MQ=255)
                                                                                                                                                                                                                              ctTCATTGAAATATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGcc                                                                                                                            <  1:715238/141‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaTATTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGATCCAAACAGGCaaaaaaa                 <  1:510093/240‑1 (MQ=255)
                                                                                                                                                                                                                                          atTTGATCGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCATGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGATCCAAACAGGCAAAAAAATGGCTGGTTTTTAtt  <  1:867849/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
CACTTTTGCCGATCTAAATTCCTGATTACCATTTCCCAGTTCAAATTTCACCTTGATATAAGCACAGATTAAATCTGCCTGTTTCACAATATGTTTTTCATCATCCGAAAAATAGTGACTATCCAGATAGGGTAAAAAATCTGCCTGCAATTCTTCCGGCAACAAACTGATCAAATGCAA‑CTCGGCAGCGGTCTCAATTTGTTTGTAAGCATGGGTAATCTCTTCATTGAAATATTTGATTGGTGTCGGCAGATCGCCGGTAAAAATTTCCGATGTGTCGTGATACATCGCCATCACCGCAATGCGCTCTGCATTCAGTTGTGCACCGTAAAAACGATTTTTAATTAAGGCTAAGGCGTGCCCGACAAATGCCACCTGTAAGCTATGTTCCGCCAGATTTTCCCGTTCAATATTACGCATAAGCGACCAGCGCTGAATCAATCTCAAGCGATCCAAACAGGCAAAAAAATGGCTGGTTTTTATT  >  1/823872‑824355

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