Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 525,571 T→G 14.9% W771G (TGG→GGG)  rhsD → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913525,5710TG14.9% 178.8 / 5.4 85W771G (TGG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (44/23);  new base G (7/10);  total (52/33)
Fisher's exact test for biased strand distribution p-value = 9.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACA  >  NC_000913/525336‑525805
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
tGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGGGGGAAACCCGGGGACGGGGGAA‑ATGGTGGGgc                                                                                                                                                                                                                                           >  2:333491/1‑239 (MQ=255)
  aGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCCAGGACGATAAAGGCCGCCTG‑ACCGGCGAAGGCCAGACGGTGGGGAACCCGGAGAACGGGGAA‑CAGCAGGGgc‑a                                                                                                                                                                                                                                         >  2:236531/1‑238 (MQ=39)
       gggggACTTGTCCACCTCTGGTCATACGATGAATGGGGACGTATCACTCACCGCACGGGGAACCGCGAACCGGCAGAGCAGTGGCAGTAGGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATga                                                                                                                                                                                                                                   <  2:110432/237‑1 (MQ=37)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGCAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGCTAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGGGGGGAAACCGGGGACGGGGGAA‑CTGCTAGGGC‑GGCATgag                                                                                                                                                                                                                                  >  2:195974/1‑239 (MQ=255)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGCTAAAGGCCGCCTG‑ACAGGCGAAAGCCACACGGGGGGGAACCCGGAAACCGGGGAA‑CTGCGGAGgc‑agaaaggg                                                                                                                                                                                                                                  >  2:177718/1‑233 (MQ=255)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCCGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATggg                                                                                                                                                                                                                                  >  1:61288/1‑237 (MQ=255)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGGGACGGGGGAA‑CTGCtgtg                                                                                                                                                                                                                                             >  1:834749/1‑229 (MQ=255)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATgag                                                                                                                                                                                                                                  >  1:969864/1‑239 (MQ=255)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGGGGAGAACCCGGAGACGGGGGAA‑CTGCTGGGGA‑CGCATgag                                                                                                                                                                                                                                  >  1:919311/1‑239 (MQ=255)
        aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCAGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGCTAAAGGCCGCCTG‑ACCGGCGACTGCCAGACGGGGGCGACCACGACGACGGGGGAA‑TTGCTGTggggaactgag                                                                                                                                                                                                                                   >  1:90024/1‑230 (MQ=255)
          ggACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGGGGC‑GGCATGAGAc                                                                                                                                                                                                                                >  1:582857/1‑239 (MQ=255)
              ttggccaccTCTGGTACTACGGTGAATCGGACCGTACCCCTCACCGCACGGTGAACGGCGAACCGGTAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGACACGTGAGCGAAGGCAAACGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGaaa                                                                                                                                                                                                                            <  1:67664/232‑1 (MQ=37)
              ttGTCATCCTCTGGTACTACGATTAATCGGACCGTATCACTCACCGCACGGTGGACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCCGTCTGCCAGACATCAGCCACCTGAGCGAAGGCCACCGCGTTGCCGTCCACTATGGCATGGACCATAAAGGCCGCCTG‑ACCCGCGAAGGACAGCCGGTGGAGACCCCGGAAAAGGGGGGA‑AAGCGGGGGC‑AGCAAGGGAAGaaa                                                                                                                                                                                                                            >  2:552109/1‑239 (MQ=255)
              ttGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGaaa                                                                                                                                                                                                                            <  1:256583/239‑1 (MQ=255)
                   aTCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGc                                                                                                                                                                                                                       <  1:154962/239‑1 (MQ=255)
                           gTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCCGTGGC‑ACCATGAGACGAAAAACGCAAACAACg                                                                                                                                                                                                               >  1:375908/1‑239 (MQ=255)
                             tctacGGTGAACCGGACCGTTTTCCCACCCGCACGGGGAACGGCAAACCGGCAGAGCAGGGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGGGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAg                                                                                                                                                                                                             <  2:969864/238‑1 (MQ=25)
                               ttcGATGAATCGGACCGTATCCCTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCa                                                                                                                                                                                                           <  1:798200/237‑1 (MQ=255)
                                   aTGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAAAACGCATAAAACGAGAAggg                                                                                                                                                                                                        >  1:837693/1‑238 (MQ=255)
                                           gACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCaa                                                                                                                                                                                                >  1:314605/1‑238 (MQ=255)
                                                tATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGCAAAACGCATACAACGGGCAGGGGCTGGCAAACcccg                                                                                                                                                                                          >  1:906033/1‑239 (MQ=255)
                                                tATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGCGCAGGGGCTGGCAAACcgc                                                                                                                                                                                           >  1:1019142/1‑238 (MQ=255)
                                                tATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGGGTTACCGTCCACTATGGCTATGACGGTAAAAGCCGCCTG‑ACCGCGGAATGCCAGACGGGGGAAGACCCGGAGACGGGGGAA‑ATGATGGGGC‑AGCATCAAACCAAACACACAAACAACAGACAGGGGAGGGCaaaagcgg                                                                                                                                                                                          >  2:760453/1‑234 (MQ=255)
                                                gtgcactcacCGGACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgcg                                                                                                                                                                                          <  1:336852/236‑1 (MQ=255)
                                                        acCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGGGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCCGGTCACGcc                                                                                                                                                                                   >  1:800639/1‑238 (MQ=255)
                                                        acCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAAGACGCATACACCGAGCCGGGGCTGGCAAACCGCGTCACGCCg                                                                                                                                                                                  >  1:292427/1‑239 (MQ=255)
                                                          cgcgcgggaaaCGGAGAACCGGGAGAGCAGGGGGGGTGTGATGGCCACGGCGGGCTGACAGAAATCGGCCACCTGAGCGAAGGCCACCGTGTGGCCGGCCGCTATGGCTGTGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGa                                                                                                                                                                                <  2:997138/230‑1 (MQ=21)
                                                          cGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGGGCCCGGGGAA‑CTGCTGTGGG‑AGCATGGGACGAAACACCCCAACAACCAGAAGGGGGTGGCAAACCGCGCCAcccccag                                                                                                                                                                                >  2:815369/1‑236 (MQ=37)
                                                          cGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGa                                                                                                                                                                                <  1:333491/239‑1 (MQ=255)
                                                          cGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACAAGCATACAACGAGCAGGGGCTGGCAAACCGCGGCACGCCgg                                                                                                                                                                                 >  1:997138/1‑238 (MQ=255)
                                                            cACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAAAACGCATACAACGGGCACGGGCGGGCAAACCCCGTCCCGCCGGAc                                                                                                                                                                               >  1:1003192/1‑238 (MQ=255)
                                                                gTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                                                                                                                             >  1:871011/1‑236 (MQ=255)
                                                                              ggCAGAGCAGTGGCAGTATGGTGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTAGGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGt                                                                                                                                                            <  1:760453/239‑1 (MQ=255)
                                                                                  gagCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCCCCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGc                                                                                                                                                         <  1:177718/238‑1 (MQ=255)
                                                                                         ggCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGGGCAGGGGCTGGCAAACCCCGTCACCCCGGACAGCCGGCCCCCGGTGGAGGGGCTGACGCAt                                                                                                                                                 >  1:396924/1‑239 (MQ=255)
                                                                                            cgtgtggtgGCCCTGGGTGGCGGACGAACATCCGCCGCCTGGGCGAAGGGCACCGGTGTGCCGCCCTCTGTGGCTATGACCATAAAGGCCGCCGG‑CCCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGGGGC‑AGCATGAGACGAAACCCGCATCCAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATgg                                                                                                                                               <  2:890224/231‑1 (MQ=17)
                                                                                             gTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGc                                                                                                                                              <  1:171029/238‑1 (MQ=255)
                                                                                              tATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTgccgc                                                                                                                                                                     <  1:167525/214‑1 (MQ=255)
                                                                                                        aCGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTg                                                                                                                                  >  1:820351/1‑239 (MQ=255)
                                                                                                        aCGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGCCGTATGGCAGCGGTTACCTg                                                                                                                                  >  1:993206/1‑239 (MQ=255)
                                                                                                          ggctggctGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCCTGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGCACCGCCGGACAGCCTGCCGCCGCGGGAGGGGCTCACGGATGGCAGGCGTGACCTCGc                                                                                                                                >  1:749279/1‑239 (MQ=21)
                                                                                                          ggctggctGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACa                                                                                                                                                                                                                  >  1:954972/1‑157 (MQ=255)
                                                                                                          ggctggctGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACa                                                                                                                                                                                                                  <  2:954972/157‑1 (MQ=255)
                                                                                                                         aTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGCAACACGCATACAACGCGCAGGGGCTGGCCAACCGCGTCACGCCGGACAGCCCGCCGCCGGGGGAGCGGCTGGCGTATGGCCGCGGTTACCTGGCGGG‑GAAGAAACgggg                                                                                                                >  2:42884/1‑239 (MQ=18)
                                                                                                                         aTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACCAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACCCCGGACAGCCTGCCGCCGGTGTACTGGCTGACGTATAGAAGCGGTCAACTGGACGG‑AATGCGGCTggg                                                                                                                >  1:1018940/1‑239 (MQ=18)
                                                                                                                             gCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACTATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCa                                                                                                                                                                                                                      <  2:6715/134‑1 (MQ=255)
                                                                                                                             gCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCa                                                                                                                                                                                                                      >  1:6715/1‑134 (MQ=255)
                                                                                                                                   tGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGGcgcc                                                                                                       >  1:234837/1‑238 (MQ=255)
                                                                                                                                       cGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGGACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACAGCGTCCCGCCGGCCAGCCTGCCGCCGGTGGAGTg                                                                                                                                                           >  1:262052/1‑183 (MQ=255)
                                                                                                                                        gAAGGCCCCCGTGTTCGCGTCCCCTGGGGCTAGGACGATAAGGGCCGCCTG‑ACCGGCGAATGCCAGACGGGGGAAAACCCGGAGACGGGGGAA‑CTGCTGGGGC‑AGCAGGAGACGAAACACGCATAAAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGt                                                                                                 <  2:234837/239‑1 (MQ=17)
                                                                                                                                                  gtgtTGCCGCCCTCTATGGCTATGAGGAAAAAGGCCGCCTG‑CCCGGCGAATGCCAGGCGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATAcg                                                                                       <  2:820351/239‑1 (MQ=255)
                                                                                                                                                   tgtTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGCACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTgc                                                                                                                                                                        >  2:266941/1‑158 (MQ=255)
                                                                                                                                                   tgtTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTgc                                                                                                                                                                        <  1:266941/158‑1 (MQ=255)
                                                                                                                                                       gCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGACAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGGCGGGTATCCGCgggc                                                                                  >  1:288660/1‑238 (MQ=255)
                                                                                                                                                       ccccgccatcttGGCTTTGACTAAAAGGGCGCCCTG‑CCCGGGGACGCCCGGAGGGGGGAGCACCGGGGGCCGGGGGAA‑CGGTGGGGGC‑AGCATGAGACGAAACACGCAACCACCGAGCAGGGGCTGGCAAACCGCGTCACGCTGGACGGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGa                                                                                  <  2:871011/228‑1 (MQ=255)
                                                                                                                                                        ccGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGAc                                                                                 >  1:313185/1‑239 (MQ=255)
                                                                                                                                                           tCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGGGCAGGGGCTGGCAAACCGCGTCACGCCGGCCAGCCTGCCGCCGGGGGGTGGGCTGCCGCATGGCAGGGGGTCCCTTGCGGG‑AGTGAAGCCGGGGCGGACGCCGCGGGTGGAGGAAACGCggggaagg                                                                              >  2:1002971/1‑234 (MQ=17)
                                                                                                                                                           tCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAg                                                                               >  1:311813/1‑238 (MQ=255)
                                                                                                                                                             ccctgtGGCTATGACGATAAAGGCCTGCTG‑CCGGGTGAATTCCAGGGGGGGGAGACCCCGGAGAGGGGGGAA‑CTGCTGGGGC‑AGCGTGAGACGAAACGCGCATACAACGGGCAGGGGCGGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGGATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCt                                                                            <  2:396924/234‑1 (MQ=255)
                                                                                                                                                             cACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTgg                                                                                                                                                                                                                                            >  1:1006182/1‑81 (MQ=38)
                                                                                                                                                             cACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTgg                                                                                                                                                                                                                                            <  2:1006182/81‑1 (MQ=38)
                                                                                                                                                              tcttggcctttGACGATAAGGGCC‑CCTGCACCGGCGAATGCCAGGGGGGGGAGAACCCGGAGTCGGGGGAG‑CTGCGGTGGC‑AGTAAGAGACGAAACCCGCATAAATCGGGCAGGGGCGGGCAACCCGCGGCACGGCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGAAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTg                                                                           <  2:265302/229‑1 (MQ=255)
                                                                                                                                                                  tGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CCGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGGGTGGCTGACGTATGGCAGCGGTTACCTCGGGGG‑GATGAAGGTGGGAGGGGCGCCGCTGGGCGGGGATAACGGGGGCAGGcgcgacc                                                                       >  2:606414/1‑233 (MQ=17)
                                                                                                                                                                     cTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGGATACCCGGGACAGGCTGCACCggg                                                                    >  1:698634/1‑237 (MQ=255)
                                                                                                                                                                        tgaacaaaaaGGCGGCGTA‑ACCGGGGGTGGCCAGGCGGCGGAGACCCCGG‑GAGGGGGGAAGCTGTGGTGGC‑GAAAGGAGACAAACCACGCATACAAGGGGGAGGGGCGGGCAATCCGCGCCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTgaga                                                                 <  2:298903/232‑1 (MQ=255)
                                                                                                                                                                        tGACGATAAAGGCCGCCTG‑ACCGGGGAATGCCAGGCGGGGGAGAACCCGGAGACGGGGGAG‑CTGCTGGGGC‑AGCATGGGACGAAACCCGCATCCAACGAGCAGGGGCGGGCAAACCGGGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTgaga                                                                 <  2:163458/239‑1 (MQ=255)
                                                                                                                                                                        tGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTgaga                                                                 >  1:163458/1‑239 (MQ=255)
                                                                                                                                                                        tGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGAATGGCACCGGTTACCTGGCGGG‑AATGAAGCGGGGCGGGACCCCGCGGGCCGAGTATACGCGGCAACGGCTGAAACGGgagc                                                                 >  1:669154/1‑238 (MQ=18)
                                                                                                                                                                        ggACGAGAAAGGCCGCCTG‑ACCGGCAAAGGCCGGACGGGGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACCCGCATACAACGGGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTgaga                                                                 <  1:42884/238‑1 (MQ=255)
                                                                                                                                                                          aCGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGGTACCTGGCGGG‑AATGAAGCCGGGCGGGACGCCGCTGGTCGGGGATACGCGGGACAGGCTGCACCGGgaga                                                                 >  1:265302/1‑237 (MQ=255)
                                                                                                                                                                              tAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGGgc                                                           >  1:298903/1‑239 (MQ=255)
                                                                                                                                                                              tAAAGGCCGCCTG‑ACCGACGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGCAAGGCTGCACCGTGAGAGCGTgc                                                           >  1:835397/1‑239 (MQ=255)
                                                                                                                                                                                       ccTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTTGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGGATACGCGGGACAGGCTGGACCGTGAGACGGTGCGAAGCTTCg                                                  >  1:1000760/1‑239 (MQ=255)
                                                                                                                                                                                       ccTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGCGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTTGCAACCCGCGTCACGCCGGCCAGCCGGCCGCCGGAGCAGTGGCTGACGTATGGCAGCGGGTACCTGGCGGGTAAGGAAG‑TGGGCGGGACGCCGCCGGTCCCGGTCAACCGGGGCAGGCGGCACCGGGGGCCGGGgcgcgagcatg                                                  >  2:605109/1‑232 (MQ=255)
                                                                                                                                                                                         gttgccGGGGAATGGCCGGGGGGGGAGACCCCGGAGACGGGGGGA‑CGGCTGGGGC‑GGCATGAGACGAAACCCGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCgg                                                 <  2:993206/235‑1 (MQ=18)
                                                                                                                                                                                                    aTGCCAGACGGTGGAGAACCCGGAGCGGGGGGAA‑CTGCGGGGGC‑AGCATGAAACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                       <  1:195974/238‑1 (MQ=255)
                                                                                                                                                                                                     tGCCAGACGGTGGAGAACCCGGGGACGGGGGAA‑CTGCGTTGGC‑AGCATGAGCCGAAACACGCATACACCGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                     <  1:895712/239‑1 (MQ=255)
                                                                                                                                                                                                                        gggAGAGGGGGGAA‑CGGCTGGGGCGAGC‑GGAGACCAAACCCGCATACACCGAGCAGGGGCTGGCAAACCGCGTCCCGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGaa                  <  1:1002971/238‑1 (MQ=18)
                                                                                                                                                                                                                         ggAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGGGGGACGCCGCTGGCCGAGTATACGCGGGACAGGCGGACCCGGGGGAGGGGGGGCAGCTCCGGGGGCATGGGAGGGAGAAAAGCCggcaaacaac                 >  2:44366/1‑230 (MQ=18)
                                                                                                                                                                                                                         ggAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCa                                                                         <  2:777903/183‑1 (MQ=255)
                                                                                                                                                                                                                         ggAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCa                                                                         >  1:777903/1‑183 (MQ=255)
                                                                                                                                                                                                                         ggAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGAGGTATGGCAGCGGGTACCTGGCGGG‑GATGAAGGTGGGGGGGACGCCGCTGGGCGAGTTAACGGCGGACCGGCGGGACCGTGGGACGGGGGGCAGCTTCGGGAGACTGGGCGGAGGTAATGCCGcgacggagg                 >  2:862851/1‑231 (MQ=17)
                                                                                                                                                                                                                           ggaCGGGGAAA‑CGGCTTGGGC‑AGCATGGGACGAACCCCGCAACCAGCGAGCAGGGGCTGGCAAACCGCGTACCGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACt                <  1:17949/237‑1 (MQ=18)
                                                                                                                                                                                                                                ggggAA‑CTGCTGGGGC‑AGCATGAGACGAAACCCGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCa           <  1:44366/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tgctgTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATCAGCGGGACAGGCTGCACCGGGAGACGGGGCGCAGCTTCCGCAGCCTGGCCGGCGGGAATGCCGCCACCGAAAGGACCAGccaacaa   >  2:1016696/1‑233 (MQ=18)
                                                                                                                                                                                                                                          ctgTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAACGGAACTGACCAggcaaaaaa  >  2:1023020/1‑230 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
TGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTG‑ACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAA‑CTGCTGTGGC‑AGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGG‑AATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACA  >  NC_000913/525336‑525805

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