Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,456,500 A→G Q76Q (CAA→CAG)  unknown → ParA‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,456,5000AG98.5% 218.3 / ‑4.2 67Q76Q (CAA→CAG) unknownParA‑like protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (33/33);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TTTCAAAGGAGACCCTTTATTATGAAAGTCATATCCGTTTTAAGTCAAAAAGGTGGTAGCGGTAAAAGTACGCTCTCAATTAATATCGCTCGTTGTTTGCAGCTCAAAGATTTTGAAGTGGCATTAATTGACACTGACCCGCAAGCGTCCGCGCGTGAATGGAATGCGCTTGCCGGTGATGATTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTCTCCGAAAAAGAAATCCGTGCGCTTGGTCGCCAAGCGGATTTTTTAATCATTGATGGTGCGCCACGTATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAA  >  1/1456252‑1456742
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
tctttaaGGAGACCCTTTATTATGAAAGTCATATCCGTATTAAGTCAGAAAGGCGGTAGCGGTAAAAGCACGCTTTCTATCAATATTGCCCGCAGCCTGCAACTTAAAGGTTTTGATGTCGCATTAATTGATACCGACCCGCAAGCGTCCGCGCGTGAATGGAATGCGCTTGCCGGTGATGATTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTCTCCGAAAAAGAAATCCGTGCGCTTAGTCGCCAAg                                                                                                                                                                                                                                                   >  1:2884/6‑250 (MQ=255)
 cttaaaGGAGGCCCTTTATCATGAAAGTTATATCTGTATTAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCg                                                                                                                                                                                                                                                 >  1:585774/4‑251 (MQ=255)
 cttaaaGGAGGCCCTTTATCATGAAAGTTATATCTGTATTAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCg                                                                                                                                                                                                                                                 >  1:658698/4‑251 (MQ=255)
            cccTTTATCATGAAAGTTATATCTGTATTAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAAt                                                                                                                                                                                                                                      >  1:544201/1‑251 (MQ=255)
                         aaGTTATATCTGTATTAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATggag                                                                                                                                                                                                                           >  1:310955/1‑247 (MQ=255)
                          agttatatCTGTATTAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATggagc                                                                                                                                                                                                                          <  1:2741/245‑4 (MQ=255)
                           gttatatCTGTATTAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATggagctc                                                                                                                                                                                                                        >  1:583782/4‑245 (MQ=255)
                               tatctgtattAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATggagctccccgc                                                                                                                                                                                                                   <  1:533598/243‑11 (MQ=255)
                               tatctgtattAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATggagctccccg                                                                                                                                                                                                                    >  1:554259/9‑241 (MQ=255)
                                 tctgtattAAGCCAGAAAGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATggagctccccgcat                                                                                                                                                                                                                 <  1:278310/245‑13 (MQ=255)
                                            ccagaaaGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCa                                                                                                                                                                                                      >  1:663622/5‑251 (MQ=255)
                                                aaaGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGAc                                                                                                                                                                                                  <  1:432086/251‑1 (MQ=255)
                                                 aaGGTGGCAGTGGTAAAAGCACGCTTTCTATCAATATTGCCCGCTGCCTGCAACTTAAAGGCTTTGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGAcc                                                                                                                                                                                                 <  1:708465/251‑1 (MQ=255)
                                                                                                ctgcaacttaaaggctttGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACt                                                                                                                                                                                   >  1:791679/16‑218 (MQ=255)
                                                                                                   caacttaaaggctttGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATcac                                                                                                                                                <  1:602088/238‑1 (MQ=255)
                                                                                                      cttaaaggctttGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGAtc                                                                                                                                                                >  1:19193/10‑230 (MQ=255)
                                                                                                            ggctttGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGAc                                                                                                                                                                                                  <  1:293382/188‑1 (MQ=255)
                                                                                                            ggctttGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGatat                                                                                                                                      >  1:163389/4‑251 (MQ=255)
                                                                                                            ggctttGATGTCGCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGat                                                                                                                                        <  1:269000/246‑1 (MQ=255)
                                                                                                                   atgtcgCATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGc                                                                                                                                <  1:241666/245‑1 (MQ=255)
                                                                                                                         cATTAATTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCaaa                                                                                                                         >  1:470649/1‑251 (MQ=255)
                                                                                                                             aaTTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTaa                                                                                                                                                                                                                                       <  1:766070/137‑1 (MQ=255)
                                                                                                                             aaTTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAAc                                                                                                                                                                                    >  1:160304/1‑188 (MQ=255)
                                                                                                                             aaTTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGAtc                                                                                                                                                                <  1:347801/208‑2 (MQ=255)
                                                                                                                             aaTTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGatat                                                                                                                                      >  1:213797/1‑234 (MQ=255)
                                                                                                                             aaTTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGAtt                                                                                                                     >  1:534074/1‑251 (MQ=255)
                                                                                                                             aaTTGATACTGACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGAtt                                                                                                                     >  1:127422/1‑251 (MQ=255)
                                                                                                                                       gACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCg                                                                                                                                                                              >  1:710526/1‑184 (MQ=255)
                                                                                                                                       gACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTcgatagag                                                                                                             >  1:687922/1‑242 (MQ=255)
                                                                                                                                       gACCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTcgatagag                                                                                                             >  1:243854/1‑242 (MQ=255)
                                                                                                                                        aCCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGa                                                                                                                       >  1:126054/1‑238 (MQ=255)
                                                                                                                                        aCCCACAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTcgatagag                                                                                                             >  1:615720/1‑241 (MQ=255)
                                                                                                                                           cacAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTg                                                                                                       <  1:249534/249‑1 (MQ=255)
                                                                                                                                           cacAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTg                                                                                                       <  1:214139/249‑1 (MQ=255)
                                                                                                                                            acAAGCGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTg                                                                                                       >  1:285738/2‑250 (MQ=255)
                                                                                                                                                 cGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGc                                                                                                  <  1:891707/250‑1 (MQ=255)
                                                                                                                                                 cGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCa                                                                                                 <  1:716751/251‑1 (MQ=255)
                                                                                                                                                 cGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCa                                                                                                 <  1:655757/251‑1 (MQ=255)
                                                                                                                                                 cGTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCa                                                                                                 <  1:832655/251‑1 (MQ=255)
                                                                                                                                                  gTCCGCGCGTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATgg                                                                                                                                                                                                                             <  1:497361/126‑1 (MQ=255)
                                                                                                                                                     cgcgcgTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTa                                                                                             <  1:651967/251‑1 (MQ=255)
                                                                                                                                                          gTGAGTGGAATGCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGaaa                                                                                                                                                       >  1:117399/1‑188 (MQ=255)
                                                                                                                                                           tgagtggaatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTg                                                                                                       <  1:807745/231‑1 (MQ=255)
                                                                                                                                                           tgagtggaatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCaaa                                                                                       >  1:87197/5‑251 (MQ=255)
                                                                                                                                                           tgagtggaatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCaaa                                                                                       >  1:380005/5‑251 (MQ=255)
                                                                                                                                                           tgagtggaatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGGGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCaa                                                                                        <  1:429625/246‑1 (MQ=255)
                                                                                                                                                             agtggaatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTc                                                                                                                                                                     >  1:874230/3‑171 (MQ=255)
                                                                                                                                                                 gaatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGCTCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATg                                                                                 >  1:425015/1‑251 (MQ=255)
                                                                                                                                                                  aatgCGCTCGCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGa                                                                                <  1:225076/251‑1 (MQ=255)
                                                                                                                                                                       gctcgCCGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGatat                                                                                                                                      <  1:476719/188‑1 (MQ=255)
                                                                                                                                                                             cGGTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAg                                                                     <  1:336443/251‑1 (MQ=255)
                                                                                                                                                                               gTGATGACTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGc                                                                    <  1:165671/250‑1 (MQ=255)
                                                                                                                                                                                        ttttCCCGGTGTTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGtc                                                          <  1:583945/251‑1 (MQ=255)
                                                                                                                                                                                                gtgtTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAAt                                                                                                                                             >  1:176844/1‑160 (MQ=255)
                                                                                                                                                                                                gtgtTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGc                                                                                                  <  1:289830/203‑1 (MQ=255)
                                                                                                                                                                                                gtgtTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGtc                                                          <  1:761184/243‑1 (MQ=255)
                                                                                                                                                                                                  gtTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATcac                                                                                                                                                >  1:386318/1‑155 (MQ=255)
                                                                                                                                                                                                  gtTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTg                                                                                                       <  1:329330/196‑1 (MQ=255)
                                                                                                                                                                                                  gtTTGCTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCaaa                                                >  1:890376/1‑251 (MQ=255)
                                                                                                                                                                                                       cTTGTGATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCaa                                                                                        >  1:823734/1‑206 (MQ=255)
                                                                                                                                                                                                          gtgATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGcc                                                                                                                                                                               >  1:215311/1‑116 (MQ=255)
                                                                                                                                                                                                          gtgATAAAGGGCTTTCCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAAc                                               <  1:320682/244‑1 (MQ=255)
                                                                                                                                                                                                                    gctttcCGAAAAAGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGTTTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATaaacaaaccaat                               <  1:259343/246‑11 (MQ=255)
                                                                                                                                                                                                                              aaaGAAATCAGGGTGCTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATaaacaaacc                                  >  1:813934/1‑230 (MQ=255)
                                                                                                                                                                                                                                         ggtgcTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttgag          <  1:661241/247‑4 (MQ=255)
                                                                                                                                                                                                                                          gtgcTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttgagtt        >  1:508998/3‑246 (MQ=255)
                                                                                                                                                                                                                                                 gTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTAttaa  <  1:461247/250‑4 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
TTTCAAAGGAGACCCTTTATTATGAAAGTCATATCCGTTTTAAGTCAAAAAGGTGGTAGCGGTAAAAGTACGCTCTCAATTAATATCGCTCGTTGTTTGCAGCTCAAAGATTTTGAAGTGGCATTAATTGACACTGACCCGCAAGCGTCCGCGCGTGAATGGAATGCGCTTGCCGGTGATGATTTTTTCCCGGTGTTTGCTTGTGATAAAGGGCTCTCCGAAAAAGAAATCCGTGCGCTTGGTCGCCAAGCGGATTTTTTAATCATTGATGGTGCGCCACGTATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAA  >  1/1456252‑1456742

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