Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,614,170 C→T F230F (TTC→TTT yjjV → putative metal‑dependent hydrolase YjjV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,1700CT97.7% 260.3 / NA 89F230F (TTC→TTTyjjVputative metal‑dependent hydrolase YjjV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (49/37);  minor base G (0/2);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTTACCGCGCACTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCG  >  NC_000913/4613924‑4614415
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cGATTTACCGCGCACTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCg                                                                                                                                                                                                                                                    >  1:424425/1‑251 (MQ=255)
 gATTTACCGGGCACTGGCGGGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGt                                                                                                                                                                                                                                                   <  2:365197/251‑1 (MQ=255)
 gATTTAACGAGCAGTGGCGCGGTTCACGGTTTTTCCGGCTGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCCTAGGCGGTACTATAACCTATCCACGCGCCCGTACAACACGCGATGTCATAGCCAAATTATAGCTGGCACAGTTATTGCGGGAAAACGACGCGCCGGATATGCCGCTCAACGGTCTTGCGGGACAGACTACACGCCCGGAGAAGCATGCACGTGTGTTTGCCgg                                                                                                                                                                                                                                                   >  2:374459/1‑250 (MQ=255)
     taccgagaaCTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGAACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCCGCCTTACCGCCCGGAGCCAGCTGCCCGGGGGTTTGCCGTGCtt                                                                                                                                                                                                                                               >  2:407098/9‑251 (MQ=255)
      aCCGCGCACTGGCGTGGTTCACGGTTTTTCTGGCAGCCTGCCACAGGCCGAACTGTTTGGACAGCGGGGCTACAAAATTGGAGTACGCGGTATTATCCCATATCCACGCGGCAGTAAACCCCGCGATGTCATCGCAGAATTTCCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATACGTCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCTCGTGTGTTTGCCGTGCtt                                                                                                                                                                                                                                               <  1:605776/250‑1 (MQ=255)
       cctcgcACTGGCGTGGTTCACTGTTTTTCCGGCAGCCTGCAACAGGCCGAACTGTTTGTACAGCTGGGCTACAAAATTGGCGGAGTCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAATCCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTg                                                                                                                                                                                                                                             <  2:464212/248‑1 (MQ=255)
       ccGCGCACTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTg                                                                                                                                                                                                                                             >  1:339030/1‑251 (MQ=255)
        cgcgCACTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCGGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGc                                                                                                                                                                                                                                            <  1:584802/251‑1 (MQ=255)
              cTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCCCAACGGTTTTCAGGGGCTGCCTAACCGCCCGGAGCAGGCTGCCCGTGGGTTTGCCGTGCTTTGCGAGTTg                                                                                                                                                                                                                                      >  2:140829/1‑251 (MQ=255)
                  cGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGcc                                                                                                                                                                                                                                  <  1:135993/251‑1 (MQ=255)
                    tGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCgg                                                                                                                                                                                                                                >  1:158288/1‑251 (MQ=255)
                         cACGGTTTTTCCGGCAGCCTTCAACAGGCCGAACGGTTTTTACATCTGGGCTACAAAATTGGCGTAGGCGGTACTATCCCCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAAcc                                                                                                                                                                                                                           <  1:462514/251‑1 (MQ=255)
                               tgttCCGGCATCATGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAATCCCGCGAGGTCAGCGCAAAATTACCGCTGGAATCGTTATTGCTGGAAACCGACGCGCCGGATATTCCGCTCAACGGTTTTCATGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACcggcgg                                                                                                                                                                                                                      <  1:572585/248‑1 (MQ=255)
                                   ccGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTGGCGCCGGGAACCGGCGGATga                                                                                                                                                                                                                  <  1:491576/250‑1 (MQ=255)
                                         gCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCCACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTgcg                                                                                                                                                                                                           >  2:418525/1‑251 (MQ=255)
                                          ccTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCAGCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTgcgc                                                                                                                                                                                                          <  1:266659/251‑1 (MQ=255)
                                             gCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCaa                                                                                                                                                                                                        <  2:151028/250‑1 (MQ=255)
                                             gCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCa                                                                                                                                                                                                         >  1:236947/1‑249 (MQ=255)
                                              cAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGAAGCGCCGGCTATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGACCGGAGCCGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAAACGGCGGGTGTGATTGCGCCAGc                                                                                                                                                                                                      >  1:363740/1‑251 (MQ=255)
                                                aCAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCCAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGGGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAgggg                                                                                                                                                                                                    >  2:278243/1‑248 (MQ=255)
                                                  aGGCCGAACGGTTGGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCCCCTATCCACGCGCCAGTAAAACCCGCGATGTCGTCGCAAAAGTACCGCTGGCATCGGTATTGCTGGAAACCGGCGCGCCGGATATTCCTGTCAACGGTTTTCAGGGGCAGCCGAACCGCCCGGAGCAGGCTGCGCGTGTGTGTGCCGGGCTTTGCGAGTTGCGCCGGGATCCGGCGGATGAGAGGGCGCAAGCGTTg                                                                                                                                                                                                  <  2:49952/251‑1 (MQ=255)
                                                   ggCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGc                                                                                                                                                                                                 <  1:577201/251‑1 (MQ=255)
                                                    gCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGCCATCGCAAAAGTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCGCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGACCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCt                                                                                                                                                                                                <  1:590413/251‑1 (MQ=255)
                                                     ccGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCCTCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGCGATTGCGCAAGCGTTGCtt                                                                                                                                                                                               >  2:204380/1‑251 (MQ=255)
                                                      cGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAACTTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGGGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCtt                                                                                                                                                                                               >  1:407924/1‑250 (MQ=255)
                                                       gAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTta                                                                                                                                                                                              >  1:539645/1‑250 (MQ=255)
                                                        aaCGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTtaat                                                                                                                                                                                            >  2:475227/1‑251 (MQ=255)
                                                            gTTTGTACAGCGGGGCAACAAAACTGGTGTAGGAGGTACTATCACCTCTCCACGCGCCAGGAAAACCCGCGATGTCATCGACAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCCGCGGATACGCCGCTCAACGGTTTTCAGGGGCCGCCGAACCGCCCGGAGCAGGCTGCCCGTGTGTGTGCCGTGCTTTGCGAGTTGCGCCCGGAACCGGCATAGGAGAGTGCGCAAGCGtgtctgaatacca                                                                                                                                                                                        <  2:86520/251‑13 (MQ=255)
                                                               tGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCCAGCGTTGCTTAATAACACg                                                                                                                                                                                      >  2:176056/1‑250 (MQ=255)
                                                                      cTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGGGATTGCGCAAGCGTTGCTTAATAACACGTATACgtt                                                                                                                                                                              >  1:381076/1‑251 (MQ=255)
                                                                        gggCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACgttgt                                                                                                                                                                            >  1:33200/1‑251 (MQ=255)
                                                                             ataaatatGGCGTAGGCGGTACTATCCCCTATCCACGCGCCATTAAAGCCCGCGATGTCATCGCAAAATTTCCGCTGGCAGCGTTATTGCTGGAATCCGACTCGCCGGAGTTGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTTCCGTTCTTTTCGAGTTGCGCCGTTAACCGGCTGATGAGATTGCGCAAGCGTTGCTTAATATCACGTATACGTTGTTTAAc                                                                                                                                                                       <  2:613259/244‑1 (MQ=255)
                                                                                aaaTTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGGGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTg                                                                                                                                                                    >  1:489713/1‑251 (MQ=255)
                                                                                  aTTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGcc                                                                                                                                                                  >  1:286117/1‑251 (MQ=255)
                                                                                  aTTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGGGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGCGATTGCGCAAGCGTTGCTTAATTACACGTATACGTTTTTTAAAGTGcc                                                                                                                                                                  >  1:630292/1‑251 (MQ=255)
                                                                                   ttGGCGTAGGCGGTCCTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCg                                                                                                                                                                 >  1:199794/1‑251 (MQ=255)
                                                                                        gTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAGCCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGGTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTGAATAACACGTATACGGTGGTTAACGTGCCGTAGGc                                                                                                                                                            <  1:228195/251‑1 (MQ=255)
                                                                                               gTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATCTGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCCAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTCGGCCgggta                                                                                                                                                      >  2:121091/1‑247 (MQ=255)
                                                                                                   tATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCg                                                                                                                                                 >  2:455694/1‑251 (MQ=255)
                                                                                                     tCACCTATCCACGCGCCAGTAAAAGCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATCGCTGGCAACCGCCGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTATCCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTCTGCGTGTTCCGCCGGGAACCGGCGGATGAGATTGCGCACGCGTTGCTTAATAACACGTATAACTTGTTTAAAGTGCCGTTGGACGGATAAAGCggt                                                                                                                                               >  2:222455/1‑249 (MQ=255)
                                                                                                           aTCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGTAAGCCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGTCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                        <  1:328697/251‑1 (MQ=255)
                                                                                                             ccACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                      >  1:580712/1‑251 (MQ=255)
                                                                                                                       gTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                                            >  2:179173/1‑251 (MQ=255)
                                                                                                                        tAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGGTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCAGtt                                                                                                                           >  1:148577/1‑251 (MQ=255)
                                                                                                                         aaaaCCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTg                                                                                                                          <  1:115853/251‑1 (MQ=255)
                                                                                                                                gcgATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGCTATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGCGTTGCGCCGGGAACCGGCGGATGAGCTTGCGCAAGCGTTGCTTAATCACACGTATCCGGTGTTTAACGTGCCGTAGGCCGGACAAGGCGTTCACCCCCCCTCCGGCCGTTGgcgccc                                                                                                                    >  2:136248/1‑248 (MQ=255)
                                                                                                                                       cATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGCGCAGGCTGCCCGTGTGTTAGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAAAACGTCTACGTTGTTTAACGTGCCGTTGGCCGGATCATGCGTTCACGCCGCATCCGGCAGTTCGCGCACCATGCCTg                                                                                                            >  1:622974/1‑251 (MQ=255)
                                                                                                                                           gtaaaaTTACCGCTGGCAGCGTTATTGCGGGAAACCGTCGCGCCGGATATGCCTCGCATCGGTTTTCAGGGGCAGCCTACCCGCTCGTAGCAGGCGGCACGTGTGTTTGCCTTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGTCCTGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCTCACAATGCCTGATGc                                                                                                        <  1:404672/249‑1 (MQ=255)
                                                                                                                                           gCAAAATTACCTCTGGCATCTTGATTCCGGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTGTCAGGGGCAGTCTAACCGCCCGGAGCAGGCTGCCCGTGTGTGTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGGGGAGATTGCGCAATCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGTTAAGGCTTTCACTCCGCATCCGGCAGTTGGCGCACAATGCCTTATGc                                                                                                        <  1:383245/251‑1 (MQ=255)
                                                                                                                                            cAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATACCACGTATACGTTGTTTAACGTGCCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCg                                                                                                       >  2:123338/1‑251 (MQ=255)
                                                                                                                                            cAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCg                                                                                                       <  1:290817/251‑1 (MQ=255)
                                                                                                                                                aTTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTATCCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGc                                                                                                   <  1:326557/251‑1 (MQ=255)
                                                                                                                                                          gcagcGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTc                                                                                         <  1:468292/247‑1 (MQ=255)
                                                                                                                                                             tCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGCTTTCAGCGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGACCCGGCGGACGAGGTTGCGCAAGAGTTGCTTAATAACACGTATACGTTGTTTACCGTGCCGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCtt                                                                                       >  1:41612/1‑250 (MQ=37)
                                                                                                                                                                  aTTGCTGGAACCCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                  <  1:478377/250‑1 (MQ=255)
                                                                                                                                                                    tggtgtaaaCCGACGCGCCGGATTTGCCGCTCATCGTTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGACGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACTTTGGTGAACGTGCCGTAGGCCGGATAAGGGTCTCACGCCGCATCCGGCAGTTGGCGCACAATGCCAGATGCGCCGCTTGGCTCGTCTTATCATGc                                                                                <  2:366512/244‑1 (MQ=37)
                                                                                                                                                                      cTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCt                                                                              >  1:160270/1‑250 (MQ=255)
                                                                                                                                                                                  gcgcCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGGGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGcc                                                                 <  2:385827/251‑1 (MQ=255)
                                                                                                                                                                                   cgcCGGATATGCCTCTCAACGGTTTACAGGGGCAGCCTACCCGCCCGGAGCAGGCGGCCTGTGTGTTTGCCGTGCTTTGCGATTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACTTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCg                                                                <  1:420888/251‑1 (MQ=25)
                                                                                                                                                                                      cGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGGGGGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                             <  2:580712/251‑1 (MQ=37)
                                                                                                                                                                                       ggATATGCCGCTCAACGGTTTACAGGGACAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACCATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaagactatgccgcacc                                                            >  2:445129/1‑236 (MQ=37)
                                                                                                                                                                                        gatagcccGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                           <  1:66773/245‑1 (MQ=37)
                                                                                                                                                                                        gATATGCCGCTAAACGGTTGTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGGGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                           <  1:6496/251‑1 (MQ=25)
                                                                                                                                                                                                   caccGGTTTTCAGGGGCATCCTAACCGCCCGTGGCAGGCTGCTCTTGTGTTTGCCGTGCTTTGCGTGTTGCGCGGGGAACCGGCGGATGAGATTGTGTAAGTGTTGCTTAAGATCACGTATACGTTGTTTACCGTGCCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGACCCGTAGGCCGGATa                                                <  2:630292/248‑1 (MQ=18)
                                                                                                                                                                                                   cAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGa                                                  <  1:222535/249‑1 (MQ=255)
                                                                                                                                                                                                    aaCGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTGGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATaa                                               <  1:341533/251‑1 (MQ=255)
                                                                                                                                                                                                     aCGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATa                                                >  1:586917/1‑249 (MQ=255)
                                                                                                                                                                                                        gTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATAAGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCg                                           >  2:440193/1‑251 (MQ=255)
                                                                                                                                                                                                         ttttCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTTTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGt                                          >  1:74620/1‑251 (MQ=255)
                                                                                                                                                                                                           ttCAGGGGCCGCCTACCCGCCCGGACCAGGCTGCACGTGCGTTTGCCGTGCTTTGCGAGCTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATACCACGTATACGTTGTTTAACGTGCCGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGCTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                        >  2:95600/1‑251 (MQ=21)
                                                                                                                                                                                                             cAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAc                                      >  1:397691/1‑251 (MQ=255)
                                                                                                                                                                                                              aGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACCCGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                                     >  2:523080/1‑251 (MQ=255)
                                                                                                                                                                                                               ggggGAGCCTAACCGCCCGGATGAGGGCGCCCGGGTGTGTGCCGGGCGTTGGGAGGTGCGCCGGGAACCGGAGGATGAGAGTGCGCAAGCGTTGCTTAATAGCTCGGATACGTTTTTTAACGTGCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTTTCGCGTCTTAGCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAc                                      <  2:381076/249‑1 (MQ=17)
                                                                                                                                                                                                               ggggCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                    >  2:564291/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGTGACGCTTGTCGCGTCTTATCATGCCTACAAATCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccgc                                 >  2:335265/1‑251 (MQ=25)
                                                                                                                                                                                                                   cAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGCTAAGGCGTTCAcgccgc                                 >  2:136247/1‑250 (MQ=255)
                                                                                                                                                                                                                      ccTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTACCGTGCCGTAGACCGGCTAAGGCGTTCACGCAGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCCTGCCTACAAGTCTGTGCCGAACCGTACGCCGGCTAAGGCGTTCACGCCGCCTcc                             >  1:623797/1‑251 (MQ=21)
                                                                                                                                                                                                                       cTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                            >  1:239582/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaCCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                           <  1:170531/250‑1 (MQ=255)
                                                                                                                                                                                                                           ccGCCCGGAGCAGGCTGCTCGTGTGTTGGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAg                        <  1:125193/251‑1 (MQ=21)
                                                                                                                                                                                                                           ccGCCCGGAGCAGGCTGCCCGGGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGGAGGCCGGATCAGGCGTTCACGCCGCATCCGGccg                        >  2:126126/1‑249 (MQ=21)
                                                                                                                                                                                                                               ccGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGtcgg                    >  1:15330/1‑248 (MQ=255)
                                                                                                                                                                                                                                    gCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATGGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGcat               <  1:335265/251‑2 (MQ=18)
                                                                                                                                                                                                                                     cAGGCTGCCGGTGTGGTGGCCGGGCTTTGCGAGGTGCCCGGGGAACCGGCGGAGGGTATGGCGCACGCGTGGCTTAATAGCCCGTGTACGTTGTTGAACGGGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTTCCATCCCTACAATTCTGTGCCGAACCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGtctgcgcata              <  2:33200/251‑10 (MQ=2)
                                                                                                                                                                                                                                     cAGGCTGCCCGTGTGTTTGCCGTGTTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGcata              <  1:278243/251‑3 (MQ=18)
                                                                                                                                                                                                                                     cAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGCTGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGcat               >  1:216096/1‑249 (MQ=18)
                                                                                                                                                                                                                                       ggCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAAt            >  1:344176/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGc          >  1:604882/1‑251 (MQ=18)
                                                                                                                                                                                                                                                 gtgtTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCg  >  1:427869/1‑251 (MQ=18)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CGATTTACCGCGCACTGGCGTGGTTCACGGTTTTTCCGGCAGCCTGCAACAGGCCGAACGGTTTGTACAGCTGGGCTACAAAATTGGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCG  >  NC_000913/4613924‑4614415

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