Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,466,733 A→G pseudogene (1342/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,7330AG97.9% 434.2 / ‑4.0 145pseudogene (1342/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (79/62);  total (80/65)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 145 total aligned reads displayed.

TCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCAT  >  NC_000913/1466441‑1467022
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                  
    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACgg                                                                                                                                                                                                                                                                                          <  1:1095355/299‑1 (MQ=255)
        gCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAGCTCGGAAGAGCGCCGCGTCGGGAAAGAGTGTAGATGCCCGCAGATGACGGATCATGACACCACACAAAACAACAACGCCAAAACCCCCCTTGACGGCAAGGGAACCCCCGCACATGAAATTGCAGGCAAAAACGGGCAGGTGGTTCAGGACGGCGATc                                                                                                                                                                                                                                                                                     >  2:755779/1‑300 (MQ=0)
         cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGCCGGCAAGGATTCCACCGGTACGGAAATTGCTGGCATCAAAGGGAAGGTGGTTCAGGCCGGCGaa                                                                                                                                                                                                                                                                                      >  2:801083/1‑297 (MQ=255)
             ctggcagTCAGGGGCGTCGGTGAGGGGATTGACAGCACCGGTGACAGGTAAACGGTAGATAACAAGGGCAGCATGACCGTCACCGATCCGGAGTCCAACGGTATCCAGGTAGACGGTGACCAGGCGGTTGTTAAGAACGAAGGCGAGAGCGCCATCGCCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATg                                                                                                                                                                                                                                                                               <  2:925104/295‑1 (MQ=17)
              tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGACCTgg                                                                                                                                                                                                                                                                                  >  1:959385/1‑297 (MQ=255)
                  tgtaggcggcggctGGAACGGTATAGATAGCACCGGCGGAAGGGCAACGGCGGATAACAAAGGCAGCATGACCGTCACCGATCCGGAGTACATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACCGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTCGATGTCAgc                                                                                                                                                                                                                                                                         <  2:1244562/296‑1 (MQ=17)
                  tGTCAGCGGCGGCGGTCACGGTATTGATATCACCAGCGACAGCGCACCGGTGGATAGCAAGGGCACTATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCCAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGAGAGCTGGagttcacg                                                                                                                                                                                                                                                                          >  2:348570/1‑294 (MQ=255)
                      aGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAgcggcg                                                                                                                                                                                                                                                                     >  1:868280/1‑301 (MQ=255)
                                gTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAAGCCAGGGGTATCGACGTGGAGGGCGAGCAGGCGGTTGTTAACAACGAACGCGATAGCGCCAACACCAACAGTGGCAACGGCACGCAGATTAACGGTGATGACGCCACGGCACACAACAACGGCAAAACCACCGTTGACGGTAAGGATTCCCCCGGTACGGAAATTGCTGGAACTCACGGGCAGGCGGTTCATGAATGCGATCATGCTGTCCgcggccgccgatcacg                                                                                                                                                                                                                                                           >  2:527201/1‑293 (MQ=17)
                                      ttcagacgtgtgctcttccgatctcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                                                                                                                                                                      <  1:1130377/276‑1 (MQ=255)
                                        tttGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                                                                                                                                                                   <  1:1100028/300‑1 (MQ=255)
                                         ttGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGGCACGGTATTGgtaa                                                                                                                                                                                                                                                  >  2:1098701/1‑300 (MQ=255)
                                          tGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAAAAACGGCAAAACCGCCGTTGACGGCAAGGATTCCGCCGGTACGGAAATTGCCGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                                                                                                                                                                                 <  1:226849/301‑1 (MQ=255)
                                           gATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGATCGTCAATAACGAAGGCGAGCGGGCGGCCACGTGCGGTGACACCGGGACGCACATTAACGGCGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGAACGGAAATTGCTGTCAAGAACGGGAAGGTGGATCAGGCCGGCGAGCTGGAGGTCGGGGGCGGGGGACACGGTAtgtaaaac                                                                                                                                                                                                                                                 >  2:632231/1‑293 (MQ=14)
                                            tatccaccgttgcgctgtcgccGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGACGACGCCCCGGCAAACAACAACGGCAACACCACCGTTGACGGCAAGGCTTCCACGGGTACGGAAATTGCGGGCAAAAACGGGATGGTGGTTCAGGCCGGCGCACTGATGGTCACAGCCGGCGCTTCCGGTTTTGTTATCa                                                                                                                                                                                                                                                >  2:1082484/22‑300 (MQ=17)
                                             tatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                                                                                                                                                                              <  1:349004/301‑1 (MQ=255)
                                             tatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCACCGGCGGAGGTCAAGGTATTGATATCAc                                                                                                                                                                                                                                               >  1:1276738/1‑300 (MQ=255)
                                               ccACCGGCGACAGCGCAACGGTGGATAACAAGGGAACCAGGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGAACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCCCCGTTGACGGCAAGGATTCCCCCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                                                                                                                                                                            <  2:1266423/300‑1 (MQ=255)
                                                 aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCg                                                                                                                                                                                                                                          <  1:795636/301‑1 (MQ=255)
                                                  ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGcta                                                                                                                                                                                                                                         >  1:556910/1‑299 (MQ=255)
                                                  ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGa                                                                                                                                                                                                                                         >  1:173992/1‑301 (MQ=255)
                                                    ggCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGATACGGAAATTGATGGCAATAACGGGAAGGTGGTGAAGGACGGCGATCTGGATGGTAGGGGCGGCGGTCACGGTAGTGGTAGCACCGATGACa                                                                                                                                                                                                                                       >  1:269618/1‑301 (MQ=17)
                                                     cacctcacctacccgacgcccttccgagctccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAg                                                                                                                                                                                                                                      <  2:1111983/271‑1 (MQ=17)
                                                        aCAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAGGGTATTGATATCACCGGTGACAgag                                                                                                                                                                                                                                    >  1:1046875/1‑298 (MQ=255)
                                                              cAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTGAATGGGTATCCAGGTTGACGGCGAAGAGGCGGGTGCTAACAACGAGGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTACACCGGTACGGAAATGGCTTGCAATAAAGGGAAGGAGGTTCAGGACGGCGAGCGGGCAGTCAGCGGGGGGGGTCACGGTATTGatacacgcggcacacacgcaacggg                                                                                                                                                                                                                             >  2:423935/1‑279 (MQ=17)
                                                               taCGTTGGATAACAAGGTCACCATGACCGTAACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCCAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTg                                                                                                                                                                                                                            <  1:948635/300‑1 (MQ=255)
                                                               aaCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTg                                                                                                                                                                                                                            <  1:712937/301‑1 (MQ=255)
                                                                aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTCACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATCACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTGACGGTATTGTACTCAACGGCGACCGCGCTACGGTga                                                                                                                                                                                                                           >  2:947483/1‑300 (MQ=17)
                                                                 attccgatcgcAAGGGCACCATGAACGTAACCGAGCAGGAGTCCAAAGGTATCCAGGTTGACGGCGACCAGGCGGTGGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGa                                                                                                                                                                                                                          <  2:1259494/291‑1 (MQ=17)
                                                                     ccaccaggaagccgcgtgtccggTCCCCGATCCGGAGGCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAACAACGAAGGCGAGCGCGCCATCCCCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                                                                                                                                                                      <  2:1112033/279‑1 (MQ=17)
                                                                       caacacgacgctcttccgatctaaCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                                                                                                                    <  2:1220743/278‑1 (MQ=255)
                                                                         aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAgg                                                                                                                                                                                                                  >  1:745338/1‑301 (MQ=255)
                                                                           cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCTTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGaagatcgg                                                                                                                                                                                                                   >  1:1261026/1‑291 (MQ=255)
                                                                                   tggagttcagacgtgtgctcttccgatctggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAc                                                                                                                                                                                                        <  1:390426/272‑1 (MQ=17)
                                                                                          cggcACCGACCCGCGGTCCCTACGTATCCAGGCTGACCGGCACCAGGCGGATGTCAACAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACAGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                                                                                                                                 <  2:594486/298‑1 (MQ=14)
                                                                                             cACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaaagatcggaagagcacacgt                                                                                                                                                                                                 >  1:1222237/1‑279 (MQ=255)
                                                                                               ccGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACagatcggaagagcacacgtctgaac                                                                                                                                                                                             >  1:1113360/1‑276 (MQ=255)
                                                                                                  ttCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                                                                                                                         <  1:527810/300‑1 (MQ=255)
                                                                                                  aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGAACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACAACCGTTGACGGCACGGATTCCACCGGTACGGAAATTGCGGGCCATCACGGGATGGTGGTTCAGGACGGCGACCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCTGCGACAGCGCAACGGTTGATAAAAAGGGCACCATGACCGTTACCGATACggg                                                                                                                                                                                         >  1:1017287/1‑300 (MQ=17)
                                                                                                   tgcGGAGTCCACCGGTATCGAGGTTGACGGGGACCAAGCGGTTGTTAACAACGAAGGCGAGAGCGCCATTACCAACGGGGCCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGATGGCAAGGATTCGAGCGATTCGGAAATTGCTGTGTATTACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                                                                                                                                                                        <  1:42825/299‑1 (MQ=17)
                                                                                                   tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                                                                                                                                                                        <  1:424425/301‑1 (MQ=255)
                                                                                                          tCCAGCGGTATCCAGGTGGACGGCGCCCAGGCGGTTCGTAACAACGAAGGCGAGAGCGCCATCACCACCGGCGGCCCCGGCACGCAGAGGAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGTAAGGTGGTTCAGGACGGCGATCTGTATGCTCGTGCCCTCCGCCACGCTCTTCCTATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCg                                                                                                                                                                                 <  2:173763/301‑1 (MQ=17)
                                                                                                           ccATCGGTATCCAGGTTTACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGAGCCTTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACAGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCTGTCACGGTAGTGATATCACCGGCGAAAGCGCTACGGTGGATAACAAGGGCGCGATGACCGTATCCGGTCCGGAgcagatcgg                                                                                                                                                                                >  1:541902/1‑293 (MQ=17)
                                                                                                           ccATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCgg                                                                                                                                                                                >  1:17166/1‑301 (MQ=255)
                                                                                                            cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                                                                                                                                                                               >  1:212757/1‑301 (MQ=255)
                                                                                                                ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCgatcggcagagcgtcgtgtagggaaagg                                                                                                                                                                           >  2:389969/1‑274 (MQ=255)
                                                                                                                      cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTg                                                                                                                                                                     >  1:584180/1‑301 (MQ=255)
                                                                                                                      cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGAAGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGCTAACAAGGGCACCATGACCGTCACCGATCCGGATTCCATCTGTATCCAGAttt                                                                                                                                                                     >  2:1246511/1‑300 (MQ=255)
                                                                                                                             aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGCCGGCGa                                                                                                                                                              >  1:1018854/1‑301 (MQ=255)
                                                                                                                                 cGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAAAAGGGCACCATGACCGTCACCGATCTGGAGTCCATCGGTATCCCGATTGACGGCGAcccg                                                                                                                                                          >  1:780015/1‑299 (MQ=255)
                                                                                                                                  tgctcttccgatcctAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAgg                                                                                                                                                         <  1:401298/287‑1 (MQ=255)
                                                                                                                                   aCCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAgg                                                                                                                                                         >  1:749859/1‑300 (MQ=255)
                                                                                                                                    ccAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGcc                                                                                                                                                       <  1:556183/301‑1 (MQ=255)
                                                                                                                                       ggCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTTCACCGGTACGGAAATTGCTGGCAATAACGGGAATGTGGTTCCGGACGGCGATCTGGCTGTCAGCGGCGGCGGTCACGGTCTTGATCTCACCGGCGACAGAGCCACGGTGGATAACAAGGGCTCCATGACCGTCCCCGCTCCTGTTTCCCTCGGTATCCAGATTCTCAGCGACACGGCCac                                                                                                                                                     >  2:333552/1‑299 (MQ=17)
                                                                                                                                             tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                                                                                                              <  1:756668/301‑1 (MQ=255)
                                                                                                                                             tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGCAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                                                                                                              >  1:926363/1‑301 (MQ=255)
                                                                                                                                             tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATATGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCTGCGACAGCGCAACGGTGGATAACAAGGGCACCATGCCCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGCCGGCTACCAGGCCACCGTCAAt                                                                                                                                              >  2:47152/1‑301 (MQ=255)
                                                                                                                                                tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCTGGAGATCgggagaacg                                                                                                                                           >  2:400841/1‑293 (MQ=255)
                                                                                                                                                  acacagAAGGCGAGAACTCCCACGACAACGGGGACAGCGCCACGCAGAATAACTATGACGCAGCCACGGCAGACAACAACGGCAACACCAGCGTTGACGGCAAGGTTTCTACCAGGACGGGGAGTGCTGCGAGAAACGCGCTGGTGCTCTAGACGACGCTCTTGCGATCTAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGa                                                                                                                                         <  2:867273/296‑1 (MQ=9)
                                                                                                                                                     acGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAgg                                                                                                                                      <  1:632950/301‑1 (MQ=255)
                                                                                                                                                        aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgg                                                                                                                                   >  1:463804/1‑300 (MQ=255)
                                                                                                                                                          ggCGAGAGCACCATCACCAACGGTGGCACCGGAACGCAGATTACCGGGAATGACGCGACGGCAAACAACAACGGGAACACTTTCGTTGTCAGCGATGCTTCCACCGCTACGGCAAATGCTCGTAACCTCGACGACGCGGTTCCGGTCTGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgag                                                                                                                                  <  2:428927/300‑1 (MQ=14)
                                                                                                                                                                 gcaccATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAAGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                                                                                                                           <  1:333940/297‑1 (MQ=255)
                                                                                                                                                                  cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                                                                          >  1:41694/1‑300 (MQ=255)
                                                                                                                                                                    ccATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                                                                                                                       <  1:47201/301‑1 (MQ=255)
                                                                                                                                                                      aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                                                                                     <  1:1248005/301‑1 (MQ=255)
                                                                                                                                                                         ccgaaCGGTGGCACCGGCACGCAGATTAACGGTGATGAAGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCGGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                                                                                                                  <  2:463273/298‑1 (MQ=255)
                                                                                                                                                                            aaCGGTGGCAGCGGCACGAAGATTAACGGTGATGACGCAACGGCAAACAACAACGGCAAAACCACCGGTGACGGCAAGGATTCCAACGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                                                               <  2:748971/301‑1 (MQ=255)
                                                                                                                                                                                  ggCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                                                                         >  1:1272393/1‑301 (MQ=255)
                                                                                                                                                                                      ccGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACTGTCACCGATCCGGATTCCAATGGTATCTAGATCGACGGCGACCAGGCCATCCTCAATAACGAAGGCGAGCGCAATATCACCAATGGccgcaccggca                                                                                                      >  2:747537/1‑300 (MQ=17)
                                                                                                                                                                                       cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcaga                                                                                                    >  1:811652/1‑299 (MQ=255)
                                                                                                                                                                                             gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACCCAGac                                                                                               >  2:1176745/1‑299 (MQ=255)
                                                                                                                                                                                                  atAACGGGGATGACGCCACGGGAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATACCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                                                                          <  2:582148/299‑1 (MQ=255)
                                                                                                                                                                                                                     ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGaaa                                                                      >  1:84698/1‑300 (MQ=255)
                                                                                                                                                                                                                      gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCCGGACGGCGATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCATCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAGCCGGAGTCCAACGGTAGCCAGAGGGACGCAGACCAGGCCACCGTCAATAACGAAGGCGATGACAATATCAACAATGGCGGCACCGGCACGCAGAGCAAAGGTAACGAACCCACcgcgtaca                                                                     >  2:939239/1‑297 (MQ=17)
                                                                                                                                                                                                                       cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCACCGGTAACGACGCCACCGCGAATaa                                                                    >  1:1011114/1‑301 (MQ=255)
                                                                                                                                                                                                                         aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                                                                  <  1:121156/301‑1 (MQ=255)
                                                                                                                                                                                                                           caccaacGCCAAAACAACCGCTGGCGGCAAAGCTGCCACGGGACCGGAAAGTACTGGCAAGACCGGGAAGGTAGTTAAGGACGGCGACCCGGAGGTCAGCTGCGCAGGCCACGGTATTGATATCACGGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGt                                                                <  2:1045636/298‑1 (MQ=17)
                                                                                                                                                                                                                            aacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAACCGGAGGCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGAGCCCATATCAACGGTCACGACGACACAGCGAATAACaagg                                                               >  2:545680/1‑298 (MQ=17)
                                                                                                                                                                                                                             aatctggCAAAACCACCGTTGAGGGAAAGGATTCCACCGGTACGGAAATAAACGGCAATAACGTGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTAACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                                                              <  2:52171/296‑1 (MQ=17)
                                                                                                                                                                                                                                  ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa                                                         >  1:353458/1‑301 (MQ=255)
                                                                                                                                                                                                                                  ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa                                                         >  1:585643/1‑301 (MQ=255)
                                                                                                                                                                                                                                   tgctcttccgatctGACGGCAAGGATTCCCCCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGATATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa                                                         <  1:33535/287‑1 (MQ=255)
                                                                                                                                                                                                                                   gCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa                                                         <  1:1178167/300‑1 (MQ=255)
                                                                                                                                                                                                                                    cAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGCGAGCACTATCACCAATGGCGGCACCGGCACTCCGCTCAACGGTAACGACGCCACCGCGAATAACCGTGGCAAAAc                                                        >  2:800089/1‑300 (MQ=255)
                                                                                                                                                                                                                                    cAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAc                                                        >  1:934323/1‑300 (MQ=255)
                                                                                                                                                                                                                                      aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCAc                                                     <  1:1246271/301‑1 (MQ=255)
                                                                                                                                                                                                                                         ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTACCGACGCCACCGCGAATAACAGTGGAAAAACCa                                                      >  1:819542/1‑297 (MQ=255)
                                                                                                                                                                                                                                         ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCa                                                      >  1:993273/1‑297 (MQ=255)
                                                                                                                                                                                                                                                     gCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGa                                      >  2:988935/1‑301 (MQ=255)
                                                                                                                                                                                                                                                      cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCCCTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGAt                                     >  2:572902/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  cgatACGGAAATTGCTGGCAAGAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTAcc                         <  1:1047620/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cGGTACGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCAAGGGTAc                          >  2:1157925/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                       cGGAAATTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGTCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTCACGACGCCTCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAACATTCCACGGGTACCaaaaa                    >  2:175821/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                           aaTTGCTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATcgcg                <  1:990134/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cTGGCAATAACGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCaa           >  1:1262592/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 tGGCAATAACGGGAAGGTGGTTCAGGAAGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATCACAAGGGCACCAGGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCaaa          >  1:883636/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                     aaaaaCGGAGAGGTGGTGCGGTACGGCGAGCTGGAGGTCGGCGGCGGAGGTCACAGTAGTGATATCAGCGGCGACAGAGCAACGGTGGATAACAAGGGCACCAGGACCTTCACCGATCCGGTGTCCGACGGTATCCAGATTGACGGCGACCAGGCGATCGTAAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCg      <  2:1017616/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aTAACGGGAAGGTGGTGCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGAAACGGTGGATAACAAGGGCACAATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCgg     <  2:992074/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         aCGGGAAGGTGGTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCCGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCAATCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGACAAAACACTGTTGAGGGAAAAGATTCCACCGGTACCAAAAACGCGAGAAATCTCGGCAt  >  2:10803/1‑301 (MQ=18)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                  
TCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCAT  >  NC_000913/1466441‑1467022

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