Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,007 A→T S181S (TCT→TCA)  rph ← truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,0070AT98.7% 271.7 / NA 78S181S (TCT→TCA) rphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (46/31);  minor base C (1/0);  total (47/31)
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

GTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCG  >  NC_000913/3815762‑3816253
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGt                                                                                                                                                                                                                                                    >  1:322144/1‑251 (MQ=255)
   tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCCCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCtt                                                                                                                                                                                                                                                 >  1:260646/1‑251 (MQ=255)
         cTTCGCTACTCATCTTACTTTTCTACAGACAACAAACAGGCGACTCATCAGTCGCCTTAACAATCAGTTTGCCAGAGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGACAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGGCTTCAACgg                                                                                                                                                                                                                                           >  2:206813/1‑250 (MQ=255)
            cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTAtt                                                                                                                                                                                                                                        >  1:414893/1‑251 (MQ=255)
                   cATCTTACTTTTCTACAGACAAAAAAAAGNCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCGTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACTTATTCCAGa                                                                                                                                                                                                                                   <  2:249094/249‑1 (MQ=255)
                            tttCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaatt                                                                                                                                                                                                                          <  1:117613/249‑3 (MQ=255)
                             ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaatttct                                                                                                                                                                                                                       >  2:322001/1‑246 (MQ=255)
                               cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttc                                                                                                                                                                                                                     >  2:511074/1‑244 (MQ=255)
                                 aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttcac                                                                                                                                                                                                                   <  1:389571/251‑10 (MQ=255)
                                 aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttc                                                                                                                                                                                                                     <  1:471834/249‑8 (MQ=255)
                                    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:238185/1‑250 (MQ=255)
                                                                                                                                                                                                             gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                          >  1:30153/1‑250 (MQ=255)
                                                                                                                                                                                                             gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                          >  1:419242/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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCG  >  NC_000913/3815762‑3816253

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