Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 90,432 A→G E113E (GAA→GAG rsmH → 16S rRNA m(4)C1402 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091390,4320AG98.8% 258.9 / ‑4.4 82E113E (GAA→GAGrsmH16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (39/41);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

GTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAG  >  NC_000913/90191‑90676
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
gTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCttt                                                                                                                                                                                                                                              <  1:473746/251‑1 (MQ=255)
    cgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCt                                                                                                                                                                                                                                          >  1:124069/1‑251 (MQ=255)
    cgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCt                                                                                                                                                                                                                                          <  1:299005/251‑1 (MQ=255)
      gtggtCACTCACGTCTGATCCTCTCGCAGCTTGGCGGAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCCATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGAGCTGGCCGAATACGTTGCCGAGCGCGATATTATCGGCAAGGTCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCttt                                                                                                                                                                                                                                        >  1:624133/1‑251 (MQ=255)
      gtggtCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCttt                                                                                                                                                                                                                                        >  1:45042/1‑251 (MQ=255)
           cactcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGCTCGCTCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGTGCGCGATCTTATCGGCCAGCTCGACGGCATTCCACTCGATCTTGGCGTCTCTTCCACGCAAATTGATGATGCTGAGCGCGGCTTTTCCTTTATgcg                                                                                                                                                                                                                                   >  2:63337/1‑251 (MQ=255)
               cacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGAt                                                                                                                                                                                                                               >  1:648552/1‑251 (MQ=255)
                ccGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATg                                                                                                                                                                                                                              <  2:17358/250‑1 (MQ=255)
                acGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGGGCGCGGCTTTTCCTTTATGCGCGATg                                                                                                                                                                                                                              >  2:292778/1‑251 (MQ=255)
                  gTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGt                                                                                                                                                                                                                            >  1:386198/1‑251 (MQ=255)
                  gTCTGAGCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGt                                                                                                                                                                                                                            <  2:187351/251‑1 (MQ=255)
                             tcGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACAt                                                                                                                                                                                                                 >  1:483831/1‑251 (MQ=255)
                               gCAGCTTGGCGAAGGGGGGCGTTTGCGGGCGAGCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCAGCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGc                                                                                                                                                                                                               <  2:388518/251‑1 (MQ=255)
                                 aGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGGAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCtttt                                                                                                                                                                                                                                             >  2:25433/1‑219 (MQ=255)
                                 aGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCtttt                                                                                                                                                                                                                                             <  1:25433/219‑1 (MQ=255)
                                  gctgggCGAAGAGGGGCGTTTGCTGGGGATCGATCGCGACCCGCAGGCTATCGGCATTGCGAGGACTATTGAGGATCCGCGCTTCTCCATCAGCCACGGACCTTTTGCCGCGCTGGTCGAATACGTTGCCGAGCGCGATCTGATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTa                                                                                                                                                                                                            <  2:578879/247‑1 (MQ=255)
                                        cGAAGAGGTGCGTTTGCTAGCGATCGATCGCGACGCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCGTCTACAGCATCCACGGACCCTTCGCCGCGCTGGCAGAATACGTTGCCGAGCGCGATCTTATCGGCAAGAGCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCGTTTCCTTGAGGCGCGATGGTCCGCTGGACATGCGTATGGAcc                                                                                                                                                                                                      <  2:276977/251‑1 (MQ=255)
                                          aaGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGAccc                                                                                                                                                                                                     >  2:193198/1‑250 (MQ=255)
                                            gagGGGCGTTTGCGGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAAc                                                                                                                                                                                                  <  2:488431/251‑1 (MQ=255)
                                             atgggCGTGTGCTGGCGAGCGATCGCGACCCGCAGGCTATCGACGTTGCGAGGACTATTGATGACCCTCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTCGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTCTCCGTGATGCGCGATGGTCCGCTGGACATTCGTATGGTCCCAAcc                                                                                                                                                                                                 <  2:104558/249‑1 (MQ=255)
                                              gggggttttGATGGCTATCGATCGCGCCCCGCAGGCTATCGCCGTTGCGAAGACTAGTGAGGATCCGCGCTTCTCCATCATCCACGGACCTTGCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAAccc                                                                                                                                                                                                <  2:686845/245‑1 (MQ=255)
                                               gggCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCg                                                                                                                                                                                               <  1:372519/251‑1 (MQ=255)
                                                   gTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTggg                                                                                                                                                                                           >  1:629906/1‑251 (MQ=255)
                                                    gtttcTGGCGATCGAGAGCGACCCGTAGGCTTTCGTAGTTGTGAAGACTATTGATGATCCGCGCTTCTCCATAATCCACGGCCCGTTCTCCGCGCGGGGCGAAGACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCCCCGCATCTTGATGATGCTGAGCGCGGCTTTTCCTGTATGCACTATGGCCCGTTGGACATGCGTATGGACCCAACCCGTgg                                                                                                                                                                                            <  2:533071/245‑1 (MQ=255)
                                                     ttGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc                                                                                                                                                                                          <  1:519821/250‑1 (MQ=255)
                                                      tGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc                                                                                                                                                                                          >  1:288159/1‑249 (MQ=255)
                                                        gtGGCGATCGGTCGCGACCCGTAGGCTACCGCCGTTGCGAAGACTTATGATGACCCGCGCTTCTCCATCAGTCACGGACCTTTCGCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGAGCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCTCCGCAACTTGATTATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGGCCGCTTGACATGCGTATGGACCCACCCCGTGGTCAGTc                                                                                                                                                                                      <  2:124069/250‑1 (MQ=255)
                                                            cgatcgatCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGcc                                                                                                                                                                                  <  1:49648/251‑1 (MQ=255)
                                                             gatcgatcGCGACCCGCAGGCTATCGCCGTTGCTAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCATCCCGTGGGCAGTCAGCCg                                                                                                                                                                                 <  2:648552/251‑1 (MQ=255)
                                                             gatcgatcGCGAACAGCAGGCTATCGACGTTGAGAAGACTATTGAAGATACGCGCTTATCCATCATCCACGGACCTATCTCCGTGCTGGGAGAATAAGGTGCCGAGCGCGCTCGTATCGGAAAGATCGCCGGCAATCTCCTCGATCTTGGCGTCTCTTCACCGCACCTTGATTCTGCTGAGCGCGGCATTTCCTTTATGCGCGATGGTACGCTGGAAAGACGTATTGACCCAACCCGTGGGCAGtgaaccg                                                                                                                                                                                 >  2:120666/1‑245 (MQ=255)
                                                                cgatcgCGGCCCGCGGGCTATCGCCGTGGCGAATACTATTGATGGTCCGCGCTTCTCCATCATCCACGGACCTTGCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCGCCTCCATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTg                                                                                                                                                                              <  2:434977/251‑1 (MQ=255)
                                                                      cgACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGCTCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCACCCCGTGGGCAGTCAGCCGCTGAATGGc                                                                                                                                                                        >  2:47219/1‑251 (MQ=255)
                                                                      cgACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCAGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCACGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGGTGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACCTGCGTATGGACCCAACCCGTGGGCAGCCCTCCGCTGAGTGGc                                                                                                                                                                        >  2:83344/1‑251 (MQ=255)
                                                                       gACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCt                                                                                                                                                                       >  1:57625/1‑251 (MQ=255)
                                                                       gACCCGCAGGCTATCGCCGTTGCGAAGACTATTGAGGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCt                                                                                                                                                                       >  1:183089/1‑251 (MQ=255)
                                                                        aCCCCCAGGCTTTCGGCGGTGCGAAGACCATTTGAGATCCGCGCCTCCTCATCAGACCCGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGAGCTTATCGGCAAGATCGGAGGCCTTCCTCCCGATCTTGGCGTCTCTTCACCGTAACCTGATGATGCCGAGCGGCGTCTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTa                                                                                                                                                                      <  2:80628/251‑1 (MQ=255)
                                                                           cGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACa                                                                                                                                                                    >  2:133734/1‑250 (MQ=255)
                                                                            gCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCTCTTCTCCATCATCCACGGACCTTGCTCCGCGCTGGGCGAATACGGTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTGGGCGTCTCTTCACCGCAGCTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGTTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACaaa                                                                                                                                                                  <  2:57051/251‑1 (MQ=255)
                                                                            gCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACaaa                                                                                                                                                                  <  1:398117/251‑1 (MQ=255)
                                                                            ccAGGCTATCGCCGGTGCGAGGTATATTGATGGTACGCGCTTATCCATCATCCACGGACCGTTCTCCGCGCTGGGCGAATACGGTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAATTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCGGGACATGTGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACaaa                                                                                                                                                                  <  2:139852/250‑1 (MQ=255)
                                                                                gCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGc                                                                                                                                                              <  1:193198/251‑1 (MQ=255)
                                                                                             cGAAGACTATTGATGATCCGCGCTTCTCCATCAGCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGat                                                                                                                                                 <  1:47219/251‑1 (MQ=255)
                                                                                              gAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGata                                                                                                                                                <  1:393685/251‑1 (MQ=255)
                                                                                              gAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGat                                                                                                                                                 >  1:138467/1‑250 (MQ=255)
                                                                                               aaGACTATTGATGATCCGCGCTTCTCCATCATCCACGGCCCTTTCTCCGCGCTGGGCGCATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCAGCTGGACATGCGTATGTACCCAACCCGTGGGAAGTCAGCCGCGGAATGGCTACAAACCGCAGAAGAAGCCGctat                                                                                                                                               >  2:638997/1‑251 (MQ=255)
                                                                                               aaGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGatat                                                                                                                                               <  1:291725/251‑1 (MQ=255)
                                                                                                       tgatgaTCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGGAGCCGAGATCGCCTggg                                                                                                                                       >  1:133514/1‑251 (MQ=255)
                                                                                                          atgatcggtgcctctcCATCAGCCTCGTACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGTGCGATCTGATCGGCAAGATCGACGGCAGTCTCCTCGATCTTGGCGTCTCTTGACCGAAACTGGATGATGCTGATCGCGGCTTTTCTTTTAGGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGGGGTCAGGCATCCGCTGAATGGCTACAAATCGCCTAAGAAGCCGATACCGTCTGGGt                                                                                                                                      <  2:457632/238‑1 (MQ=255)
                                                                                                          aggttCCGTGCTCCTCCCTCATCCACTGTCCTTTCCCCGCGCTGGGTGAATACTTTTCCGAGCGCGATCTTACCGTCAAGATCGACGGCATGCTCCCCGAGCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCTCGGGACATTCGTATGGGCCCAACCCGTGTGCAGTCAGCCTCTGAATGGCTACAAACCGCAGAAGTAGCCGATATCGCCGGGGTa                                                                                                                                     <  2:638987/247‑1 (MQ=255)
                                                                                                             atCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTg                                                                                                                                  >  1:123751/1‑251 (MQ=255)
                                                                                                                cgcgCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCACCCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGaa                                                                                                                                >  2:183155/1‑250 (MQ=255)
                                                                                                                  cgcTTCTCCACCATCCTCGGACCGTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTAGCGGCAAGATCGACGGCATTCCCCTCGATCTTGGCGTCTCTTCACCGCAACGTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAAc                                                                                                                             <  2:553456/251‑1 (MQ=255)
                                                                                                                   ccTTCTCCATCATACCCGTACCTGTCTCCGTGCGGGGCGAATTCGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGAGGGTCCGCTGGACATGCGTATGGCCCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAAcc                                                                                                                            <  2:386229/250‑1 (MQ=255)
                                                                                                                     ttccccATCATCCCCGGACCTTTCTCCGCGCGGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGaagaag                                                                                                                                                      <  2:53044/219‑1 (MQ=255)
                                                                                                                     ttCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGaagaag                                                                                                                                                      >  1:53044/1‑223 (MQ=255)
                                                                                                                         ccATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGAGCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGTTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATgg                                                                                                                       <  1:83344/250‑1 (MQ=255)
                                                                                                                              atcCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGa                                                                                                                    >  1:453738/1‑248 (MQ=255)
                                                                                                                                 cACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACAGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCg                                                                                                              >  1:278464/1‑251 (MQ=255)
                                                                                                                                   cGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGCTATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCggt                                                                                                            >  2:19779/1‑249 (MQ=255)
                                                                                                                                    ggACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGGGAAGAGCGttt                                                                                                           >  2:207276/1‑251 (MQ=255)
                                                                                                                                      aCCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTg                                                                                                         >  1:632228/1‑251 (MQ=255)
                                                                                                                                                           aTACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGcc                                                                                    >  1:30079/1‑250 (MQ=255)
                                                                                                                                                                    cGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGc                                                                          <  1:68661/251‑1 (MQ=255)
                                                                                                                                                                                    ggCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGGGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACacac                                                          >  2:84375/1‑248 (MQ=255)
                                                                                                                                                                                     gCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCg                                                         >  2:173367/1‑251 (MQ=255)
                                                                                                                                                                                         gATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTCTGGCCCCAACCCGTGGGCAGTCAGCCGCTGCATGGCTACAAACCGACGTAGAAGCCGATATCGCCGGGGTATTGAAATCCTAGGGTGAAGCGCGTTTTGCCCAACCGCTTGTCCGCGCCATTTTCGCGCTCAACCGCGAACATCCggtg                                                      >  2:565082/1‑247 (MQ=255)
                                                                                                                                                                                         gATCGACGGAATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGCGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAACCCGAAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGACAACCGCATTGCCAGCGCAATTGTCGAGCGTAACCGCGAACAGCCGATg                                                      >  1:68159/1‑250 (MQ=255)
                                                                                                                                                                                          aTCGACGGCTTTCTCCCCGAGTTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGAc                                                    <  2:288159/251‑1 (MQ=255)
                                                                                                                                                                                          aTCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTg                                                                                                                     >  1:484223/1‑187 (MQ=255)
                                                                                                                                                                                          aTCGACGGCATTCTCCACGATCTTGGCGTCTCTTCCCCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATTGCCTGGGTATTGAAAACCTATGGTg                                                                                                                     <  2:484223/187‑1 (MQ=255)
                                                                                                                                                                                               cGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCa                                               <  1:500017/251‑1 (MQ=255)
                                                                                                                                                                                                   aTTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATAATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCa                                            >  1:200850/1‑250 (MQ=255)
                                                                                                                                                                                                        ccctcGATCTTGGTTTCTCTTCACCGCAATTTGAGGATTCTGAGCGCGGCTTTCCCTTTATGCGCGATGGCCCGCTGGACATGCGTATGGACCCAACCCGTGGGCCGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGaa                                       <  2:629906/250‑1 (MQ=255)
                                                                                                                                                                                                               tCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCgg                                >  1:498775/1‑251 (MQ=255)
                                                                                                                                                                                                                  tggggTCTCTGCACCGCAACTGGGTGATGCTGAGCGGGGTTTTGCCTTTATGCGCGTTGGTCCGCTGGACATGCGAATGGACCCAACCCGGGGGCAGTCAGACGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGGGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGa                               <  2:512010/245‑1 (MQ=255)
                                                                                                                                                                                                                  tGGCGTCTCTGCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTGATGCGCGATGGTCCGCTGGACAGGCGGATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGc                                  <  2:66074/246‑1 (MQ=255)
                                                                                                                                                                                                                        ctctTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTTGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTgg                       >  1:673671/1‑251 (MQ=255)
                                                                                                                                                                                                                         tctTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgc                      >  1:192337/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCAGAAGTCGTGGCTGCTGCa               <  1:34314/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tgatgCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGTGCGTAACCGCGAACAGCCGATGACCCGCACCACAGAACTGGCGGAAGGCGAGGCTGCTGCAACGACGGTgg     >  2:575939/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gatgCTGAGCGCGGCTTTTCCTTTATGCGCGATGGGCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGa     <  1:571671/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAg  <  1:19779/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAG  >  NC_000913/90191‑90676

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