Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,793 C→T R89R (CGC→CGT insA7 → IS1 family repressor protein InsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,7930CT96.2% 236.4 / NA 81R89R (CGC→CGTinsA7IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (48/29);  minor base A (3/0);  total (52/29)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGT  >  NC_000913/4518547‑4519036
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
aTGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTc                                                                                                                                                                                                                                                 >  2:257713/1‑251 (MQ=255)
       tcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCACAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAACTCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTATCGTCCCTTAAAAAACTCAGGACGCAGTCGGTAAcc                                                                                                                                                                                                                                          >  2:22265/1‑251 (MQ=18)
              ccgtaacgggtgggggggggCGTAACGGCAAAAGCTCCGCCGTACATCAGCGCTATCGCTGTTCTCACTGCTGTAAAACATGGCAACCGCAGTTCACTTACACCGGGTTTCAACCCGGGACGCACCAGAAAATCATTGATAGGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGTCTCAACACGGTTTGACGGCACTTAAAAAACGCAGGCCGTTGGCGGTATCCTcgcga                                                                                                                                                                                                                                    <  2:402373/231‑2 (MQ=255)
                 cTACTGCAGGCGTGGGGCGCAACGACAAAAGCAATGCCGGACATCAGCGCTATCTCTGCGCTCATTGCCGTAAAACATGGCAACTAAAGTTCACTTACACCGCCTCTCAGCCAGGTACGCACCCGCAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGACCTCAACACGGTTTAACGGCACTTCAAAAACTCCGACCGAAGTAGGCAACATCGCGCATAc                                                                                                                                                                                                                                >  2:318694/1‑251 (MQ=255)
                     tGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTCCGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGAATACAAcc                                                                                                                                                                                                                            >  2:416177/1‑251 (MQ=255)
                                tgCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGACTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatg                                                                                                                                                                                                                 >  1:53771/1‑251 (MQ=255)
                                    tAACGGCAAAAGCACTGCAGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGAGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCCATTACAAAACTCAGGCCGAATGCGGTAACCACACGCCTCCAACCGGGCAgagatctgc                                                                                                                                                                                                              >  2:527949/1‑246 (MQ=17)
                                                  cTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGaa                                                                                                                                                                                                >  1:20402/1‑250 (MQ=255)
                                                    gCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATg                                                                                                                                                                                             >  1:289915/1‑251 (MQ=255)
                                                        gACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGTTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACg                                                                                                                                                                                         >  1:361940/1‑251 (MQ=255)
                                                         aCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGGTTTACGTCCCTTAAAAAACTCAGGCCGTGGTCGATGACCTCGCGCATACAACCGGGCAGAGATGGGATTGACTGCCATGCAATGGACg                                                                                                                                                                                         >  2:494490/1‑250 (MQ=17)
                                                             cAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGAGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCAATGGACGAAc                                                                                                                                                                                      >  2:13274/1‑249 (MQ=255)
                                                               gcgcTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTg                                                                                                                                                                                  >  1:507770/1‑251 (MQ=255)
                                                                cgcttgctctGCCCTCATTGCCGTAAAACAGGGCAACTACAGTTACCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTgg                                                                                                                                                                                 <  1:324983/245‑1 (MQ=18)
                                                                      ccctcctctCATTGCCTTAAAACATGCCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTAc                                                                                                                                                                           <  2:261068/246‑1 (MQ=255)
                                                                        ctGCTCTCATTGCCGTAAAACATGGCCAATACAGTTAACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATAGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAACACCTCAGGCCTTCGTCGGAAACCTCGCGCATACAGCCGGACGGTGATGTGATTGTCTGCGATGAAACGGCCGACCGGTGGGGCTACgg                                                                                                                                                                         >  2:219263/1‑250 (MQ=17)
                                                                         tGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGAAGTCGGTAACCTCCCGCATACAACCGGGCAGTGTTGTGCTTGTCTGCGCTGAAATGGAAGAACGATGGGGCTActcc                                                                                                                                                                        >  2:329160/1‑248 (MQ=18)
                                                                          gCTCTCATCGCCGAAAAACTGGGCAACTACAGTTCCATTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGACATGGCCATGAATGGCGTCGGATTTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACgcc                                                                                                                                                                        <  1:417646/250‑3 (MQ=255)
                                                                          cctctcatggCCGTAAAACATGGCAACTACAGTTCCCGTGACCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTc                                                                                                                                                                        <  2:169527/241‑1 (MQ=18)
                                                                           ctctCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCgg                                                                                                                                                                      >  1:596941/1‑251 (MQ=255)
                                                                               cattGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGt                                                                                                                                                                     <  1:430381/246‑1 (MQ=255)
                                                                               cattGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCt                                                                                                                                                                  >  1:450883/3‑251 (MQ=255)
                                                                                     cGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATc                                                                                                                                                             >  1:675454/1‑250 (MQ=255)
                                                                                     cGTAAAACATGGAAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATc                                                                                                                                                             >  1:486985/1‑250 (MQ=255)
                                                                                      gTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCAc                                                                                                                                                           >  1:417208/1‑251 (MQ=255)
                                                                                          aacaaagCAACTACAGTTCACTTCCACCGACTCTCAGCCCGGTACGCACCAGAAAATCATTGATAGGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGAATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGAGATTGTCTGCGCTGAAATGGAAGAAAAGTGGGGCTACGTCGGAGCTAAATAACGTc                                                                                                                                                        >  2:389165/7‑250 (MQ=255)
                                                                                            cATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAgc                                                                                                                                                     >  1:479800/1‑251 (MQ=255)
                                                                                             aTGGCAACTACAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAgcg                                                                                                                                                    <  1:136689/251‑1 (MQ=255)
                                                                                                 cAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCgctgg                                                                                                                                                >  1:275894/1‑251 (MQ=255)
                                                                                                   aCTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGAGGGGCTACGTCGGTGCTAAATCACGTAAGCgctggct                                                                                                                                              >  1:134061/1‑251 (MQ=255)
                                                                                                      aCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAACTCACGTCAGCGCTGGCTGtt                                                                                                                                           >  1:195316/1‑251 (MQ=255)
                                                                                                      aCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGtt                                                                                                                                           <  1:416177/251‑1 (MQ=255)
                                                                                                       cAGTTCACTTACACAGCCTCTCAGCCCGGTACGAACCAGAAAATCATTGCTATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCCCGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGTCGAACAGTGGGGCTACGACGCTGCTAAATAACGTCAGCGCTGGCTGttt                                                                                                                                          >  2:305507/1‑251 (MQ=18)
                                                                                                                 tACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTAt                                                                                                                                 >  1:665050/1‑251 (MQ=255)
                                                                                                                  acacCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTAGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGGCAGCCCTTGCTGGTTTAcgcgattt                                                                                                                                >  2:597295/1‑247 (MQ=18)
                                                                                                                     ccGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGAc                                                                                                                              >  1:436495/1‑250 (MQ=255)
                                                                                                                      cGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTACCCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGGGGGGCTCCGTCGGTGCTAAATCACGTCAGCGGTGGCTGTTTTACGTGTCTGAc                                                                                                                              >  2:214023/1‑249 (MQ=18)
                                                                                                                             cAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACgg                                                                                                                     >  1:8208/1‑251 (MQ=255)
                                                                                                                              aGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGCTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACgg                                                                                                                     >  1:579290/1‑250 (MQ=255)
                                                                                                                              aGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCgctg                                                                                                                                                 <  1:340145/222‑1 (MQ=255)
                                                                                                                                cccGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCTCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGTACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCTTATGACAGGATACggag                                                                                                                   <  1:245217/250‑1 (MQ=255)
                                                                                                                                cccGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGCACAGTGGGGCTACGTCGGTGCTAAATCAAGTCAGCGCTGGCTGATTTACGCGTATGACAGGATACggagg                                                                                                                  >  1:327029/1‑251 (MQ=255)
                                                                                                                                cccGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGAGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGAGACggagg                                                                                                                  >  1:405996/1‑251 (MQ=255)
                                                                                                                                      acgCACCAGAAAAGCATTGATATGGCCATGAATGGCGTTGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGCAGGACGGta                                                                                                            <  2:20402/249‑2 (MQ=255)
                                                                                                                                         cACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGATTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTACGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTAAGCGCTGGCTGTTTTACGAGTATGACAGGATACGGAGGACGGTtgtg                                                                                                         >  1:247069/1‑251 (MQ=255)
                                                                                                                                                 aaTCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGt                                                                                                 >  1:93710/1‑251 (MQ=255)
                                                                                                                                                        gATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGt                                                                                          <  1:608144/251‑1 (MQ=255)
                                                                                                                                                                cATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCAc                                                                                  <  1:13274/251‑1 (MQ=255)
                                                                                                                                                                cATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTATTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCAc                                                                                  <  1:257713/251‑1 (MQ=255)
                                                                                                                                                                 aTGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAAc                                                                                                                                                                                                                             >  1:316604/1‑111 (MQ=255)
                                                                                                                                                                 aTGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAAc                                                                                                                                                                                                                             <  2:316604/111‑1 (MQ=255)
                                                                                                                                                                   gAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTTCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCActc                                                                                <  1:370500/250‑1 (MQ=255)
                                                                                                                                                                   gAATGGCGGCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACGCAGGCCTTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCActct                                                                               <  2:596941/251‑1 (MQ=255)
                                                                                                                                                                    aaTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTg                                                                              >  1:460297/1‑251 (MQ=255)
                                                                                                                                                                      ggtcGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGTCCTTAACTCGGTTTTATGTCACCTAAAAAACTCAGGCCTTAGTCGGTAACCTCGCGCTCACATCCGGGCATTGATGTGATTGTCGGCGCTGAAATGGACGAACAGTGGGGATACGTCGTTGATAAATCACGTCAGGGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGAACGTCTTCGGTGAACGCGCTCTGGc                                                                            <  2:148038/248‑1 (MQ=255)
                                                                                                                                                                          gTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAACTGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGGTACGGAGGACGGTTGTGGCGCACGTCTTCggcgatcggaagagcggcggg                                                                        >  2:384326/1‑232 (MQ=255)
                                                                                                                                                                          gTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCgg                                                                                           <  1:384326/232‑1 (MQ=255)
                                                                                                                                                                             ggATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGACACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCGGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTg                                                                     <  2:195316/251‑1 (MQ=255)
                                                                                                                                                                             ggATGTCGCGCCAGTGCACGCATTATGGGCGTTCGCCTCAACACGGTTTTACGTCCCTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTATGCGTATGACAGGATACGGAGGAGGGTTGTGGCGACCGTCTTCGGTGAAAGCCCTCGGGCcacgatg                                                                     >  2:77213/1‑247 (MQ=255)
                                                                                                                                                                              gATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTgg                                                                    >  1:193502/1‑251 (MQ=255)
                                                                                                                                                                               gcccataatgaGTGCACGCCTTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACCCTgg                                                                    >  1:431349/11‑250 (MQ=255)
                                                                                                                                                                                tGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGACCGCACTCTGGCCACACTgg                                                                    >  1:132989/1‑249 (MQ=255)
                                                                                                                                                                                                gCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCctgc                                                  <  1:432513/251‑1 (MQ=255)
                                                                                                                                                                                                  aTTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCctgctg                                                <  1:305507/251‑1 (MQ=255)
                                                                                                                                                                                                       gggCGTGGGCCTCAACACGGTTTTTCGTCACTTAAAAAACTCAGGCCGTAGGCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTTTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGGCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGTCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCgga                                           <  2:436495/251‑2 (MQ=21)
                                                                                                                                                                                                          cGTTGGCCTGAACACGGTTTTACGTCTCTGAAAAAACTCAGGACGTATTCGGTAACCTCGTGCATACCACCGGGCAGTGATGTGATTGTCTGCGCTGAATTGGACGAACAGTGGGGCTACGTCGGTGCGAAATCACGTCAGCGCTTGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCtt                                        <  2:75880/251‑1 (MQ=255)
                                                                                                                                                                                                           ggtGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTGCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCttt                                       <  1:39101/249‑1 (MQ=255)
                                                                                                                                                                                                           gTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCt                                                                                                                                                                             >  1:409638/1‑117 (MQ=255)
                                                                                                                                                                                                              ggCCTCAACACGGTTTTACGTCACTTAAAAATCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCATGGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCGACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGGATGACAGGATACGGATGACTGTTGTGGCGCACGTCGTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGctttctggctgtctttgag                                    <  2:612830/251‑19 (MQ=255)
                                                                                                                                                                                                              ggCCTCAACACGGTTTTACGTCACTGAAAAAACTCAGGCCGTCTTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCTCGTCTTCGGTTAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGa                                     <  2:395023/250‑1 (MQ=255)
                                                                                                                                                                                                                ccTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAggg                                  >  1:288771/1‑250 (MQ=255)
                                                                                                                                                                                                                ccTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCACATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                  >  1:598360/1‑251 (MQ=255)
                                                                                                                                                                                                                 cTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTTTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                 >  1:375402/1‑251 (MQ=255)
                                                                                                                                                                                                                   cAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                               <  1:527949/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaCACGGCTTTACGTCACTTAAAAACCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGCCGAACAGTGGGGCACCGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGAAGGTTGTGGCGCACGTTTTCGGTGCACGCAATCTAGCCACCATGGAGCGTCTTCAGAGCTCGCTGTCGGCCATGGTGGTCGTg                              >  2:31085/1‑251 (MQ=18)
                                                                                                                                                                                                                     ccacGGTTTTACGTGACTTAAAAAACTCAGGGCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAACCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCATGCTGGCCACACTGGAGTGTCTTCTGAGCCGGCTGTCGGCCTTTGAGGTCGTgg                             <  2:145022/250‑1 (MQ=255)
                                                                                                                                                                                                                      cacGGTTTTACGTCACTTAAAAAACGCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCCCACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTgggttcggggt                            >  2:360953/1‑241 (MQ=255)
                                                                                                                                                                                                                                cGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACg                  >  1:260157/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGAtgg             >  1:254480/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTCTAGATGACGGATGGCTGGCCg      >  1:363876/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGt  >  1:158286/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
ATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGT  >  NC_000913/4518547‑4519036

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