Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,495 C→T 7.5% pseudogene (1104/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,4950CT 7.5% 236.2 / 6.4 85pseudogene (1104/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (41/36);  new base T (4/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT  >  NC_000913/1466207‑1466780
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                               
cgcgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGc                                                                                                                                                                                                                                                                                                 >  1:589988/1‑289 (MQ=255)
           aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAAc                                                                                                                                                                                                                                                                                      <  2:589988/289‑1 (MQ=255)
             aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGCCAGCGCAACgg                                                                                                                                                                                                                                                                                    >  1:565792/1‑289 (MQ=255)
                  gggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAACGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑At                                                                                                                                                                                                                                                                              >  1:774304/1‑289 (MQ=255)
                    gCACAATGACCATCACCGATCCGGAGTCTGTGGGTATCCAGATCGATGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAGCATCACCAACGGTAGCACCGGCACCCAGATTAATGGTGATAACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGTGACAGCGCGACGGTGG‑At                                                                                                                                                                                                                                                                              <  1:982430/287‑1 (MQ=17)
                     cACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAAc                                                                                                                                                                                                                                                                           >  1:893688/1‑289 (MQ=255)
                          tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCCAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑CTCACCGGCGACAGCGCAACGGTGG‑ATAACAAgg                                                                                                                                                                                                                                                                       >  1:310565/1‑288 (MQ=255)
                          tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGTCGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAgg                                                                                                                                                                                                                                                                       >  2:623905/1‑288 (MQ=255)
                              agTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAATGGAGGCACCGGCACTCAGATCAACGGTGACGACGCCACCGCCAATAACCACGGACAAACAAACGTTGATGGCAAAGATTCCAAGGGGACCAAAATAGCGGGAAATAACGGAATCGTAAATAGGGAAGGCAACCAGGCAGTCAGAGGAGGTGCGCACGGTGTTGAG‑AGCAGTGGCGACacgggcacgggcacgcacctgagaag                                                                                                                                                                                                                                                                    >  2:7186/2‑263 (MQ=255)
                                  aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAggg                                                                                                                                                                                                                                                                      >  1:59299/1‑281 (MQ=255)
                                  aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAggg                                                                                                                                                                                                                                                                      <  2:59299/281‑1 (MQ=255)
                                  aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATg                                                                                                                                                                                                                                                               >  1:62236/1‑288 (MQ=255)
                                   ccGACCCGGAGTCCATTGGTAGCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGAc                                                                                                                                                                                                                                                             <  2:893688/289‑1 (MQ=255)
                                      aTCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGCTAACTCATTTAGGCCTGTGACCTTGCCGACAGATGCGGCGGTCACGGTATTGAT‑AGCACCGGCTACAGCGCAAAGGTGG‑ATAACAAGCGACACATGACCg                                                                                                                                                                                                                                                           >  2:606442/1‑288 (MQ=17)
                                       tCCGGAGTCCGTCGGTATCCAGGTCCACGGTGACAAGGCCTGCGTCAACAACGAAGGCGAGTGCCCCATCACCGACGGTGGTACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGTAATTAATGGTAATAACTGTAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTTACTGTATTGAT‑ATAGCCGGCGACAAGGGAACGGTGG‑ATAACAAGGGCACCATGACCGt                                                                                                                                                                                                                                                          <  2:774304/288‑1 (MQ=11)
                                              tCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGACGGAGGTCACGGTCTTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGAACTATGACAGTCACCGag                                                                                                                                                                                                                                                   >  2:115132/1‑287 (MQ=255)
                                                cattgcgatctcggtCGACAGAGACAAAGCCATCGTTAATAACAGAGCGGCGAGCACCAGCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGAAAACAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGTTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACTGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                                                                <  2:1152965/275‑1 (MQ=11)
                                                cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                                                                >  1:229153/1‑289 (MQ=255)
                                                        tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTACCCGATCCTGAGTCCa                                                                                                                                                                                                                                         >  1:1152965/1‑288 (MQ=255)
                                                        tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGAGGCGGCGGTCCCGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑CTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCa                                                                                                                                                                                                                                         >  1:619256/1‑288 (MQ=255)
                                                             aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATTAACGGTGATGACGCCACGGCAAATAACAACGGAAAAACCACCGTTGATGGAAAGGATTCCACCGGTACCGAAATCGACGGTAATAACGGAAATGTAAATCTGGATGGTGACCTGGATGTCACCGGCGGTGCGCAGGGTATTGAT‑AACCCGGCCGTCAGAGCAACGTGAG‑ATAACAGAGGAACCATGCACGTTACCGATCAGGAGTCCCTCGGTa                                                                                                                                                                                                                                   >  2:1016365/1‑289 (MQ=255)
                                                             aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                                                                                                                                                                                                                                   >  1:368180/1‑289 (MQ=255)
                                                               cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                                                                                                                                                                           <  1:261809/279‑1 (MQ=255)
                                                               cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                                                                                                                                                                           >  2:261809/1‑279 (MQ=255)
                                                               cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                                                                                                                                                                                                                                  >  1:296119/1‑288 (MQ=255)
                                                                  cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGCTGTCAGCTGCGGCGGTCCCGGGATTGAG‑ATCACCGGCGACAGCGCAACGGTGG‑CTAACCAGGGCACCATGACCGTCACCGATCCGTAGTCCATCGGGATCAAg                                                                                                                                                                                                                              >  2:1203253/1‑289 (MQ=17)
                                                                        cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGa                                                                                                                                                                                                                         <  1:528243/288‑1 (MQ=255)
                                                                          aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCCCGGCAAACAACAACGGCAAAACCACCCGACAAGGCAAGGAGTCCACCGCTACGGAAATCAACGGTCATAACGGAAAAGTTATTAAGGACGGCGATCAGGATGTCAGCGGCGGCGGTGACGGGATCGAT‑ATCACCGGCGACAGCGCAAACGGGG‑ATAACAAGGGCACCATGACCGGCACCGCACGGGAGTACATGGGTACCAAgggggaagg                                                                                                                                                                                                                      >  2:523234/1‑282 (MQ=255)
                                                                                  gTCAATACCGCAGGAGAGAGGACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACTCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCAACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAgg                                                                                                                                                                                                              <  1:1152332/289‑1 (MQ=255)
                                                                                   tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGTCGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAgg                                                                                                                                                                                                              <  1:623905/288‑1 (MQ=255)
                                                                                      aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgt                                                                                                                                                                                                          <  2:296119/289‑1 (MQ=255)
                                                                                                 gagaGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                                                                                                <  1:115132/288‑1 (MQ=255)
                                                                                                    agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGATGGCGACCAGGCGGTTGTTAACAACGAAg                                                                                                                                                                                            >  1:887661/1‑289 (MQ=255)
                                                                                                    agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAg                                                                                                                                                                                            >  1:671766/1‑289 (MQ=255)
                                                                                                             aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                                                                                                                                                                   <  2:368180/289‑1 (MQ=255)
                                                                                                                  cGGTGGCACCTGCCCCCATATAAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACCATc                                                                                                                                                                              <  1:952506/289‑1 (MQ=255)
                                                                                                                         accggcaccCAGAGTAATGGGGATAAAGGAACGGCAAACAACAACGCCAAAACCACCGTTTACGGAAAGGATTCCACCGGTACGTAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAG‑ATCAACGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                                                                                                                                                       <  1:586186/289‑1 (MQ=17)
                                                                                                                                   aGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                                                                                                                                                       >  2:228649/1‑279 (MQ=255)
                                                                                                                                   aGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                                                                                                                                                       <  1:228649/279‑1 (MQ=255)
                                                                                                                                       gaaTGGTGAGAACCACACGGCAAACAACAACGGAAAACCCACGGTTGACGGCAAGTGTTCCACCGGTACGGAAATCAACGGTAATAACGTAAAATTTTTTCATGACGGCGATCTGGATGTCAGAGGCCGCGGCCACGCTATTGGTCGGCACC‑GCGAGATGGCTACAAGAG‑ACAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAcc                                                                                                                                                         <  2:310565/288‑2 (MQ=17)
                                                                                                                                               atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAAc                                                                                                                                                                                                                                                                                      <  1:632318/157‑1 (MQ=255)
                                                                                                                                               atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAAc                                                                                                                                                                                                                                                                                      >  2:632318/1‑157 (MQ=255)
                                                                                                                                                 taCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATACCACCGGTACGGAAATCAACGATAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                               <  1:523234/288‑1 (MQ=255)
                                                                                                                                                 gaCGCCACGGCAAACACGACAGACAAAAACACCGTAGCGGTCAGTGTTTCCAGCGGTACGGAAATCAACGGTAATATCGGAAAGGTATTTCAGTTCTGTGATGTGGATGTGAGGGCCAGCGGTAACTGTACCGCT‑CTATCCGGCGGCAGCGTCAGGGTGGTATAA‑GAGAGAGCCATGACCGTAACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                               <  2:62236/289‑1 (MQ=255)
                                                                                                                                                      cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGCAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGCTGTCGCCTATACACACCTGACTGTT‑CCGACGAATAGCGAGGCATCGATGG‑CTAACGACGCCACCACGACCATCACCGAACCGGAGTACAGCGGTATCGAGGTTGACGGCGACAAGTCGGTTGTTAACAACAAAGGCCAGAGCAACAGCACCAACGGTGGCCCAGGCAAGcgagtgaaccggcgc                                                                                                                                          >  2:150822/1‑275 (MQ=9)
                                                                                                                                                      cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                                                                                                                                         <  1:335916/274‑1 (MQ=255)
                                                                                                                                                      cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                                                                                                                                         >  2:335916/1‑274 (MQ=255)
                                                                                                                                                         ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATaa                                                                                                                                                                                                                                                                            <  1:870766/156‑1 (MQ=255)
                                                                                                                                                         ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATaa                                                                                                                                                                                                                                                                            >  2:870766/1‑156 (MQ=255)
                                                                                                                                                            aaaCAACAAAGGCAAAACCACCGTTGACGGAAAGGAGTCCACCGGTACGTAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGc                                                                                                                                     <  1:1073102/288‑1 (MQ=255)
                                                                                                                                                             aacaacaatGGCAAAACCACCGTGGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGTTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   <  2:531758/289‑1 (MQ=255)
                                                                                                                                                             aacaacaacAGCAAAACCACCGTTGGCAGCAAGGATTCCAGGGGTACGGAAATCAACGGGAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGGGACAGCGCAACGATGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   <  1:271827/289‑1 (MQ=255)
                                                                                                                                                             aacaaaaacGGGAAAGCCACCGTTGACGGCAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   <  1:964474/289‑1 (MQ=255)
                                                                                                                                                                     tggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGTGACAGCGCGACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGTGACCAGGCCATCGTTAATAACGAAGGCGAGAGCGCCATCACAAACGGAGGCACCGGCACGTAGAGTagaggactcgtgacttgccgacgc                                                                                                                           >  2:586186/2‑266 (MQ=14)
                                                                                                                                                                     cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGGCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGCGGCACCGGCACGCAGATTAACGGTAACGACGCCACGGCAaac                                                                                                                           >  1:129714/1‑289 (MQ=255)
                                                                                                                                                                                     gACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGCGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGAGCCATAACCAGCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGAAAACAACAAAGGCAAAacc                                                                                                            >  1:483209/1‑288 (MQ=255)
                                                                                                                                                                                     gACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTCGTCAACAACGAAGGCGAACGCGCCATCACCAACGGCGGCACCGGAACGCCGATTAACGGGAACGACGCACCCGCAAATAACAACGGCAAAaaca                                                                                                           >  2:1248611/1‑289 (MQ=14)
                                                                                                                                                                                           aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaaca                                                                                                                       >  2:603172/1‑271 (MQ=255)
                                                                                                                                                                                           aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaaca                                                                                                                       <  1:603172/271‑1 (MQ=255)
                                                                                                                                                                                                 tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   >  1:189526/1‑253 (MQ=255)
                                                                                                                                                                                                 tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   <  2:189526/253‑1 (MQ=255)
                                                                                                                                                                                                        gTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTCAACAACGAAGGCGAGAGGGCCATCACCAACGGCGGCACCGGCACGCAGCTTAACGGTAACGCGGCCACGGCAAAAAACAACGGCAAAACCACG‑GTTGATGGCAATGAtta                                                                                       >  2:593345/1‑288 (MQ=14)
                                                                                                                                                                                                         tACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCCCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCCGGTTGACGGGGACCAGGCGGTTGTTAACAACGACGGCGAGAGCGATATCATCAAAGGCCGCCCCGGCCCGCAGCTCCATGGCAACGCCGCCACGGGAAACAACAATGGCAAAACCACT‑GTTGCAGGCAAGGAtcca                                                                                      >  2:1186053/1‑286 (MQ=14)
                                                                                                                                                                                                                           taaCGGAAAAGTTGTCCAGTACGGGGATCGGGTTGTCGGCGTCGGCGGTTACGGTTTAGAG‑AGAACCGGCGAAAGGGCAACGGTTG‑ATAACAAGGGCACCATGACCGGCACCGATCAGGAGTCCACTGGTATCCACGTCGGTGTCGAACGGGCAGTAGAGAAAAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                    <  2:887661/289‑1 (MQ=14)
                                                                                                                                                                                                                           taaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGGGGTCACGGTATTGTT‑CTCACCGGCGACAGGGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                    <  1:606442/289‑1 (MQ=255)
                                                                                                                                                                                                                               ggAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGAACCATGACGGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAAGCAGATGAACGGTGAGGACCACACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACCGAAATTGCTGGCa                                                                >  1:874925/1‑289 (MQ=255)
                                                                                                                                                                                                                                  aaaGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATc                                                                                                                                                                                                                                       <  1:979353/120‑1 (MQ=35)
                                                                                                                                                                                                                                                   cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCACCGGTGG‑ATAACCAGGACACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGCGCGCCATCACCAACGGTGGCCCCGGCCCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCTACCGGTACGGAACTTGCTGGCAACAACGGGAAGGTGATTCAg                                            >  1:588622/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      tccgggtGTCAGGGGCGCGGGTAACGGTCAGGAT‑ATCACCGGCGACAGCCCAACGGTGG‑ATAACAAGGGCACCACGACCGTAACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGTCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCGCC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                          <  1:449336/282‑1 (MQ=18)
                                                                                                                                                                                                                                                      tCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                          <  1:1203253/288‑1 (MQ=255)
                                                                                                                                                                                                                                                        tGGATGGCAGCGGCGGCGGTCACGGTATTGAT‑ATAACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAGCAACGGTGGAACCGGCTCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACgg                                       <  1:150822/289‑1 (MQ=255)
                                                                                                                                                                                                                                                             gTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   >  2:366831/1‑193 (MQ=255)
                                                                                                                                                                                                                                                             gTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                   <  1:366831/193‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGTGACAGCGCGACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                                                                                                       >  2:867062/1‑106 (MQ=12)
                                                                                                                                                                                                                                                                aGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                                                                                                       <  1:867062/106‑1 (MQ=32)
                                                                                                                                                                                                                                                                aGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTgg                               >  1:188482/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                         atCGCCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa       <  1:1163430/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          tCACCGGGGCAAGCGAAACGAGGG‑AGAAGAAGGGCACCACGCCCGTCACCTATCCGGAGGCCATCGGTATCCAGGTGGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGAACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt        <  2:1245808/287‑1 (MQ=17)
                                                                                                                                                                                                                                                                                            aCCGGTGACAGCGCGACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                                                                                                                                                                   >  1:845729/1‑115 (MQ=255)
                                                                                                                                                                                                                                                                                            aCCGGTGACAGCGCGACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                                                                                                                                                                   <  2:845729/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            aCCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATACCGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCCGCGGTCACGGTATTg    >  1:1245808/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                            aCCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                         >  2:1031207/1‑236 (MQ=255)
                                                                                                                                                                                                                                                                                            aCCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                         <  1:1031207/236‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            aCCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATCAACGGTGATGACGCCACGGCAAACCACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCAACCGGTACGGAAGTTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGCATTCCGCGGCGGCGGTCACGGTATTg    >  2:1163430/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                            aCCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTCGTTAACAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACTCAGATCAACGGTAATGACGCCACCGCGAACAACCAGTGCCTAATAACCATTTGAC‑GCAGGGGCCGAACAGAGACGGATAGTTCTCGCAATAACGGGTAGGTGCATCAGGAAGGCGATCTTGATGCAAGCGGCGGCGGTCACGGTATTg    >  2:323451/1‑289 (MQ=2)
                                                                                                                                                                                                                                                                                              cGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAACTGGATGTCAGCGGCGGCGGTCCGGGTATTGat  >  1:153872/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                               
CGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT‑ATCACCGGCGACAGCGCAACGGTGG‑ATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAT  >  NC_000913/1466207‑1466780

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