Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,043 G→A G169G (GGC→GGT)  rph ← truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,0430GA96.2% 235.9 / NA 80G169G (GGC→GGT) rphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (49/27);  minor base C (2/1);  total (51/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTC  >  NC_000913/3815798‑3816285
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                                                                                                                                                                                                >  2:182970/1‑251 (MQ=255)
gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                                                                                                                                                                                                >  1:155674/1‑251 (MQ=255)
   aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGGCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGttt                                                                                                                                                                                                                                              >  2:422116/1‑249 (MQ=255)
    aaacaaaGGCGCCGCATCAGCCGCCTTAAAAATCAGTTGGCCAGCGCCGTCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCGC‑GGGCCAGCTCCAACAAGGTGAGTAGCTCTTCATGGGTTCACGGCTCGCCGTCTGCCATCCCCTGCACTGCAATTATGCGCCCGGCTTCGGTCATCACTCCGTGCATGTCGGTCTCTGCGGCTTATTCTTCAACGGATTCCAGATCGCAAAGTGCTTCACCGtttac                                                                                                                                                                                                                                            <  2:210915/247‑4 (MQ=255)
        aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTCCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCACCC‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCAATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAt                                                                                                                                                                                                                                         >  1:183388/1‑250 (MQ=255)
         aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtac                                                                                                                                                                                                                                       >  1:262752/1‑249 (MQ=255)
           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                                                                                                                                                                     <  1:440499/251‑1 (MQ=255)
           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                                                                                                                                                                     <  1:344055/251‑1 (MQ=255)
             cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                                                                                                                                                                   >  1:197973/1‑251 (MQ=255)
              gaatcataaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCATC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                                                                                                                                                                   >  1:293667/4‑250 (MQ=255)
                       gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCg                                                                                                                                                                                                                         >  1:510088/1‑251 (MQ=255)
                         cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGc                                                                                                                                                                                                                       >  2:384908/1‑251 (MQ=255)
                             ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGAATCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTaac                                                                                                                                                                                                                   >  2:303507/1‑249 (MQ=255)
                               aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGCTTACAATACCGACAGACACTGCGGCTAccct                                                                                                                                                                                                                 >  2:142792/1‑249 (MQ=255)
                                 aaaTCAGTTTGCAAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCACAAAGGCGAGTAGCTCTTGATGGGTGAACGGCTCGCCTTCTTCCGTCCCCTGCCCTTTAATTGTGCCCCCATCGATGGTGTTCACTACCTTCATGTCGGTCTATTCTGCTGACTCTTCAATGTATCGCAGGTCCCAAATTGCTTCAGCGTTNACAATCCCGCCAGAAACTGCGGccaccaaac                                                                                                                                                                                                               >  2:255360/1‑248 (MQ=255)
                                          gtcccgcgcCGCATTCTGCGTCGCTACAATGTATTCGATTC‑CCCCTC‑GGGCCAGCGCAAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCTAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCCTTCATCg                                                                                                                                                                                                      <  2:488228/246‑1 (MQ=255)
                                           gCCAGCGCCGCCTTCTTCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGGCTTCACCGTATTCCAGATCGCAACTTGCTTCCCCGTTTACACTACCTACAGAAACTGCGGCTACCACCCCTTTCCTCgg                                                                                                                                                                                                     >  2:66177/1‑251 (MQ=255)
                                              aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAAAGAAACTGCCGCCACCATCCCTTTCATCggggt                                                                                                                                                                                                  >  2:285705/1‑248 (MQ=255)
                                              aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCgggtt                                                                                                                                                                                                  <  1:322001/251‑4 (MQ=255)
                                                  ccgccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTACCCTCCCTTTTATCCggttggtt                                                                                                                                                                                              >  2:37546/1‑251 (MQ=255)
                                                  ccgccTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCGggttggtt                                                                                                                                                                                              >  1:308149/1‑251 (MQ=255)
                                                         cTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                                                                                                                                                                       >  2:378524/1‑251 (MQ=255)
                                                                  acaagagctgcgatgcc‑ccc‑TC‑GGGCCAGCGCCAACATGGGTAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTTTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCTTCCCTTTCAACGGGTTGGTTTTCAGCTTGCCGGTTTc                                                                                                                                                                             <  2:366404/235‑1 (MQ=255)
                                                                   aCAAGGGATTCGATTC‑CCC‑CCCGGGCCAGCGCCAACCAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCGCTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATATCGACAGAAACTGCGGCTACAATCCCTTTCATCGGGTTGGTGTTCAGCTTGCCGtttt                                                                                                                                                                              <  2:193/250‑1 (MQ=255)
                                                                    cAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTTCCGTTTTcc                                                                                                                                                                            >  2:19336/1‑251 (MQ=255)
                                                                            tcgattc‑CCCCTC‑GGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGACACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCtt                                                                                                                                                                    >  2:64700/1‑251 (MQ=255)
                                                                                         c‑GGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCGTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGATCATCTCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGTCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATc                                                                                                                                                        <  1:10630/251‑1 (MQ=255)
                                                                                            ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                                                                                                                                                                       <  1:53690/218‑1 (MQ=255)
                                                                                            ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                                                                                                                                                                       >  2:53690/1‑218 (MQ=255)
                                                                                            ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtg                                                                                                                                                      >  1:45310/1‑250 (MQ=255)
                                                                                             gCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtgc                                                                                                                                                     >  1:314931/1‑249 (MQ=255)
                                                                                             gCCAGCGCCAAAAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCGGACTTCCCCTGCACTTCAATCATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCGGCGGCTGAGTCTTCAACGGATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTGCCGCCAGCTTCTGTAGCGCATCtgc                                                                                                                                                     <  2:440436/251‑3 (MQ=255)
                                                                                                agcgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtgcc                                                                                                                                                    >  1:126448/4‑246 (MQ=255)
                                                                                                agcgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAg                                                                                                                                                  <  1:127587/248‑1 (MQ=255)
                                                                                                  cgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAgcg                                                                                                                                                >  1:209484/2‑251 (MQ=255)
                                                                                                      aacaaggcgaatAGCTTTTCATGGGTGAACGGCTCCCCTTCTGCCGTCCCCTGCACCTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGTTCTCTGCGGCTGAGGCTTCAACGTATTCCAGATCGCAAATTGCTTCTCCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCGTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCGTCTGTAGCGCATCTGCCAGCGCCAc                                                                                                                                            <  1:255360/240‑1 (MQ=255)
                                                                                                      aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCAc                                                                                                                                            <  1:511074/251‑1 (MQ=255)
                                                                                                              gAGTAGCTCTTAATGGGTGACCGGCTCGCTTTCTGCCGTCCCCTGCACTTCAATGATGAGCCCGTCTTCGGTCATCACTACGTTCCTGGCGGTTTCTGCGGCTGAGTCTTAAACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTTGCGCATCTGCCAGCGCCACCCAGGCa                                                                                                                                     <  2:103748/250‑1 (MQ=255)
                                                                                                                gTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCAccc                                                                                                                                  <  1:119354/251‑1 (MQ=255)
                                                                                                                  aGCTCTTCATGGGTGAACGGCTCGCCTTCTGNCGTCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCAGCTGCCAGCGCCACGCAGGCACCCGt                                                                                                                                <  2:237768/251‑1 (MQ=255)
                                                                                                                  aGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTNGCGCATCTGCCAGCGCCACGCAGGCACCCGt                                                                                                                                >  2:272513/1‑251 (MQ=255)
                                                                                                                        tCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCATCGCCACGCAGGCACCCGTAATCg                                                                                                                           >  2:17170/1‑250 (MQ=255)
                                                                                                                        tCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAAt                                                                                                                             <  2:510088/248‑1 (MQ=255)
                                                                                                                         cATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCCGGCACCCGTAATCGAc                                                                                                                         >  2:385577/1‑251 (MQ=255)
                                                                                                                              gtgagctgCTCGCCTTCTGCCGTCCCCTGCGCTTCAATTATGCGCCCGTCGTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTGCAACGTATTCCAGATCGCAAATTGCTTCCCCGCTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGt                                                                                                                    <  2:405514/244‑1 (MQ=255)
                                                                                                                                gAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTAc                                                                                                                  <  1:385577/251‑1 (MQ=255)
                                                                                                                                 aaCGGCTCGCCTTCTGCCCTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCCAATTTCTTCACCGTTTACAATAACGACAGAAACTGCGGCTACCATCCCCTTCCTTGGGTTGGTTTTCCGCTTGCCGTTTTCCACCCGCTTCTGTAGCGCCTCTTCCAGCGCCCCGCAGGCACCCGTACTCGACCCGGTAcg                                                                                                                 >  2:163497/1‑251 (MQ=255)
                                                                                                                                     gCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTg                                                                                                             >  1:225500/1‑251 (MQ=255)
                                                                                                                                                gCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGc                                                                                                  <  2:94070/251‑1 (MQ=255)
                                                                                                                                                            aCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTc                                                                                      >  2:429319/1‑251 (MQ=255)
                                                                                                                                                                    gatgCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGCCGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTccc                                                                              >  2:403602/1‑250 (MQ=255)
                                                                                                                                                                         gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                          >  1:419242/1‑250 (MQ=255)
                                                                                                                                                                         gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                          >  1:30153/1‑250 (MQ=255)
                                                                                                                                                                         gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                          >  1:238185/1‑250 (MQ=255)
                                                                                                                                                                             gTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                                                     <  1:422116/251‑1 (MQ=255)
                                                                                                                                                                               cTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTg                                                                   >  2:182807/1‑251 (MQ=255)
                                                                                                                                                                                 tCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGa                                                                  <  1:178271/250‑1 (MQ=255)
                                                                                                                                                                                   ggTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACt                                                               >  1:186754/1‑251 (MQ=255)
                                                                                                                                                                                    gttatcaCTACGTTCATGGCGGTCTCTGCGGCTGAGTCTTCAACGGATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTGCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTc                                                              <  2:480168/248‑1 (MQ=255)
                                                                                                                                                                                    gTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAGTTGCTTCACCGTTTACAATACCGCCAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCCGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGCAGCACGTCGCAGTCCAGCGTAATGGTGAACTc                                                              >  2:51956/1‑251 (MQ=255)
                                                                                                                                                                                    gTCATCACTACGTTCATGGCGGTCTCTGCGGCTGAGTCTGCAACGGATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTGCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTc                                                              <  1:51956/251‑1 (MQ=255)
                                                                                                                                                                                      catcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCCCCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAc                                                            >  2:223776/1‑251 (MQ=255)
                                                                                                                                                                                             aCGTTCATGTCGGTCTCTGCGGCTGATTCGTCTACGTATTCCAGATCCCAAATTTCATCACCGTTTACAATGCCGACAGAACCTGCGGCTCCCATCCCTTTCATCGGTTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAgc                                                       <  2:356335/249‑1 (MQ=255)
                                                                                                                                                                                             aCGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAgcgc                                                     >  2:212336/1‑251 (MQ=255)
                                                                                                                                                                                               gTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCCCCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATTGTGAACTCCCCCAGCGCtt                                                   >  2:277001/1‑251 (MQ=255)
                                                                                                                                                                                                ttCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGGTGGTTTTTAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCCTCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCttt                                                  >  2:37986/1‑251 (MQ=255)
                                                                                                                                                                                                               gCGGCTGAGTCTTCAACGTATTCCAGATCGNAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCcg                                    <  2:262752/250‑1 (MQ=255)
                                                                                                                                                                                                                cGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCTTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCcgcg                                  <  1:160727/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggctgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGAGGCTACCATCCCTTTCATAGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGAGAACTCACCCAGCGCTTTCAAATCTACCGCcgcgc                                 >  1:305858/5‑251 (MQ=255)
                                                                                                                                                                                                                  gctgagTCTTCAACGTATTCCAGATCGCAACTTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCCCCCGTAATCGACGCGGTACGCGTGCCACCATCCGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCCCCCACCGCTTTCAAATCTACTGCcgcgcg                                >  2:229944/4‑251 (MQ=255)
                                                                                                                                                                                                                   ctgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAACTCTACTGCcgcgcg                                >  1:276738/3‑250 (MQ=255)
                                                                                                                                                                                                                   ctgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCcgcgcgc                               >  2:468175/3‑250 (MQ=255)
                                                                                                                                                                                                                       gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAagag                           <  1:226028/251‑1 (MQ=255)
                                                                                                                                                                                                                         cTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGCACTCCCCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGc                          >  2:143417/1‑250 (MQ=255)
                                                                                                                                                                                                                              actgaTTCCAGATCGCAAATGTCGTCACCGTTTACAATTCCGACAGAAACTGCCGCTACCATCCCTTTCAGCGGGTTGGTTGTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCGGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTGCAAATCTACTGCCGCGCGAAGAGCACGGGc                    <  2:314931/247‑1 (MQ=255)
                                                                                                                                                                                                                              aCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGc                    >  1:61217/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGCAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACg           >  2:433527/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGACGCACGTCGCAGTCCAGCGTAATGGGGAACTCACCCAGCGCTTTCAAATCCNCTGCCGCGCGAAGAGCACGGGCGCTCAGACg           >  2:265941/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGc          >  1:226726/1‑251 (MQ=255)
                                                                                                                                                                                                                                                attgCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTc  >  1:465675/4‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
GACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTC  >  NC_000913/3815798‑3816285

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: ‑