Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,250 G→C 5.4% G194G (GGG→GGC dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2500GC 5.4% 709.9 / 22.3 205G194G (GGG→GGCdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (98/96);  new base C (5/6);  total (103/102)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 208 total aligned reads displayed.

ATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGT  >  NC_000913/4261982‑4262521
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                               
gacaggagaTCAACCGGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                                                                                                                                                                                            <  2:19866/269‑1 (MQ=255)
                          ggcactcttGAGCGGGTGCAGCACGGCAGGTTTGGTGCGTGTGTGTTGGGTAAGGCGGAGCTGGTTGCCGACTGGGGGAAAGGGATGCGCTACGTGGTGTCGTTGCCGGTGACGGTGAAAACGCGTATTGGCGTGGACGACCAGAACAGCTATGAATTTCTCTGCGAGTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCGGAAAGAAAACCGTGAAATCccgcc                                                                                                                                                                                                                                                 <  4:194012/267‑1 (MQ=255)
                                        ggTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                                                                                                                                                                                                      >  2:419155/1‑240 (MQ=255)
                                             aGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCg                                                                                                                                                                                                                                                                        <  2:307284/233‑1 (MQ=255)
                                                    cATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                                                                                                                                                                                                       >  2:211063/1‑275 (MQ=255)
                                                                                    aGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                                                                                                                                                                                        >  2:280447/1‑274 (MQ=255)
                                                                                                                          gtgtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGt                                                                                                                                                                                                                                                                                >  6:658581/1‑148 (MQ=255)
                                                                                                                           tgtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTTATCAATACCGTTTCCGGCAAAGGCGTGTGTGAGATGTTCATCATCCGCGCACGTAAAGCCTGGCTTTCGGGGTTa                                                                                                                                                                                                                                                                              >  7:566168/1‑149 (MQ=255)
                                                                                                                            gtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTaa                                                                                                                                                                                                                                                                             >  7:514208/1‑149 (MQ=255)
                                                                                                                            gtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTa                                                                                                                                                                                                                                                                              >  5:14677/1‑148 (MQ=255)
                                                                                                                             tCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                                                                                                                                                                                            >  7:141814/1‑149 (MQ=255)
                                                                                                                                aTTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGGTTAACGGTGGTATCAAGTCGCTGGAAGAGGCAAAAGCAAACCt                                                                                                                                           >  1:151983/1‑275 (MQ=255)
                                                                                                                                   ccGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                                                         <  2:612824/274‑1 (MQ=255)
                                                                                                                                    cggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                                                                                                                                                                                                     <  7:659654/149‑1 (MQ=255)
                                                                                                                                    cggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCCCACCTGCaa                                                                                                                                       >  2:9851/1‑275 (MQ=255)
                                                                                                                                    cggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      <  6:233716/148‑1 (MQ=255)
                                                                                                                                       tgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcc                                                                                                                                                                                                                                                                  <  8:633357/149‑3 (MQ=255)
                                                                                                                                        gacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGc                                                                                                                                                                                                                                                                           <  6:429645/139‑1 (MQ=255)
                                                                                                                                         acggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgt                                                                                                                                                                                                                                 >  4:525898/1‑180 (MQ=255)
                                                                                                                                             tgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTTAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGt                                                                                                                                                                                                                                                            <  8:747941/149‑1 (MQ=255)
                                                                                                                                               aaaaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                                                                                                                                                                                                          <  8:112595/149‑1 (MQ=255)
                                                                                                                                                   cgcgTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGCTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGGGAAATCCCGCCGCCCGATTATCCGCGTGTGTAGCAACTGAAGCGTGACTTTCGGCACCTGACAATTTCGATTAACGGTGGTATCAAGTCGCTGGAAAAGGCCAAAGAACACATGCAAAATATGGATATCGTg                                                                                                                        >  4:511770/1‑275 (MQ=255)
                                                                                                                                                      gTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                                                                                                                                                                                                    <  6:479125/148‑1 (MQ=255)
                                                                                                                                                        aTTTGCATCGACGACCAGGACAGCTATGCATTTCGCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCCTAACCCACGACCGTAAAGCCTGGGTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAACCCCGCCGCTGGATTATCCGCGTTTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGt                                                                                                                   <  1:137922/275‑1 (MQ=255)
                                                                                                                                                                cgacgaCCAGAACAGCTATGAATTGCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTTGGGcgcg                                                                                                           <  2:382825/275‑1 (MQ=255)
                                                                                                                                                                   cgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                                                                                                                                                                                       >  5:6199/1‑148 (MQ=255)
                                                                                                                                                                      ccAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                                                                                                                                                                                                    >  4:659993/1‑132 (MQ=255)
                                                                                                                                                                          gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                                                                                                                                                                                                <  6:763941/148‑1 (MQ=255)
                                                                                                                                                                          gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                                                                                                                                                                                                <  6:400522/148‑1 (MQ=255)
                                                                                                                                                                                tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                                                                                                                                                                                                          >  5:575678/1‑148 (MQ=255)
                                                                                                                                                                                 ctGAATTTCTCTGCGTTTTCATCAATACCGTTTCCGGCAAAGGCGAGAGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATcc                                                                                           <  2:511772/273‑1 (MQ=255)
                                                                                                                                                                                 aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                                                                                                                                                                                                                         >  5:834188/1‑148 (MQ=255)
                                                                                                                                                                                 aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                                                                                                                                                                                                                         >  5:530162/1‑148 (MQ=255)
                                                                                                                                                                                    aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                                                                                                                                                                                                      >  5:534338/1‑148 (MQ=255)
                                                                                                                                                                                    aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTg                                                                                                                                                                                                                     >  7:706808/1‑149 (MQ=255)
                                                                                                                                                                                         ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCacac                                                                                                                                             >  2:105359/1‑216 (MQ=255)
                                                                                                                                                                                           ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGt                                                                                                                                                                                                               >  5:660354/1‑148 (MQ=255)
                                                                                                                                                                                            tgagATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGGGTGTATCAACTGAAGCGTGa                                                                                                                                                                                                             >  7:604142/4‑149 (MQ=255)
                                                                                                                                                                                            tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                                                                                                                                                                                             <  7:28141/149‑1 (MQ=255)
                                                                                                                                                                                            tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGAc                                                                                                                                                                                              >  2:660953/1‑164 (MQ=255)
                                                                                                                                                                                             gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      >  6:839310/1‑91 (MQ=255)
                                                                                                                                                                                             gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      <  5:839310/91‑1 (MQ=255)
                                                                                                                                                                                             gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGAACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAg                                                                                                                                                  <  2:344421/207‑1 (MQ=255)
                                                                                                                                                                                                 tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                                                                                                                                                                                                          >  8:243647/1‑131 (MQ=255)
                                                                                                                                                                                                 tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                                                                                                                                                                                                          <  7:243647/131‑1 (MQ=255)
                                                                                                                                                                                                 tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACttt                                                                                                                                                                                                         >  6:35436/1‑148 (MQ=255)
                                                                                                                                                                                                      tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTaa                                                                                                                                                                                                                                                                             <  7:374226/75‑1 (MQ=255)
                                                                                                                                                                                                      tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                                                                                                                                                                                                          <  5:701582/78‑1 (MQ=255)
                                                                                                                                                                                                      tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGCCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCAt                                                                                                                                                                                                   >  8:861892/1‑149 (MQ=255)
                                                                                                                                                                                                        aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt                                                                                                                                                                                                 >  8:565647/1‑149 (MQ=255)
                                                                                                                                                                                                           aCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGa                                                                                                                                                                                               >  6:310413/1‑148 (MQ=255)
                                                                                                                                                                                                           aCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGa                                                                                                                                                                                               <  5:35436/148‑1 (MQ=255)
                                                                                                                                                                                                           aCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGAc                                                                                                                                                                                              >  1:123056/1‑149 (MQ=255)
                                                                                                                                                                                                                ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATg                                                                                                                                                                                          >  5:893309/1‑148 (MQ=255)
                                                                                                                                                                                                                   cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg                                                                                                                                                                                       <  5:11907/148‑1 (MQ=255)
                                                                                                                                                                                                                    ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa                                                                                                                                                                                      <  6:688990/148‑1 (MQ=255)
                                                                                                                                                                                                                      cAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTa                                                                                                                                                                                   <  8:65536/149‑1 (MQ=255)
                                                                                                                                                                                                                         aGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGACGATTAAc                                                                                                                                                                                 <  6:6199/148‑1 (MQ=255)
                                                                                                                                                                                                                           gCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg                                                                                                                                                                               <  5:71955/148‑1 (MQ=255)
                                                                                                                                                                                                                              aGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGACAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATAAACTGAAGCGTGCCTTTNAGCATCTNACANTGTCGTNTa                                                                                                                                                                                   >  3:36164/1‑141 (MQ=255)
                                                                                                                                                                                                                                  tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTGCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg                                                                                                                                                                               <  1:36165/141‑1 (MQ=255)
                                                                                                                                                                                                                                   gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt                                                                                                                                                                              >  6:129083/1‑141 (MQ=255)
                                                                                                                                                                                                                                   gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCa                                                                                                                                                                       <  5:195855/148‑1 (MQ=255)
                                                                                                                                                                                                                                    agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCaa                                                                                                                                                                      >  5:283417/1‑148 (MQ=255)
                                                                                                                                                                                                                                     gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaaa                                                                                                                                                                                                                                                        >  8:579240/1‑65 (MQ=255)
                                                                                                                                                                                                                                     gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGGGaaa                                                                                                                                                                                                                                                        <  7:579240/65‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCaa                                                                                                                                                                      <  7:224180/144‑1 (MQ=255)
                                                                                                                                                                                                                                         ttCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt                                                                                                                                                                >  8:838686/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt                                                                                                                                                                <  6:530972/148‑1 (MQ=255)
                                                                                                                                                                                                                                            ctcgtcccccagcgTAAAGCCGGGCTTTCGGGGTTAAGCCCGAAAGAAAACCTTGAAATCCCGCCGCTTGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGGCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGGCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGGTGCCGATcc                               <  4:292380/263‑1 (MQ=255)
                                                                                                                                                                                                                                            atcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGaa                                                                                                                                                            <  2:176962/150‑1 (MQ=255)
                                                                                                                                                                                                                                             tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGa                                                                                                                                                             <  6:893309/148‑1 (MQ=255)
                                                                                                                                                                                                                                             tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGa                                                                                                                                                             <  5:487128/148‑1 (MQ=255)
                                                                                                                                                                                                                                                tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                                                                                                                                                                                       >  6:666053/1‑71 (MQ=255)
                                                                                                                                                                                                                                                tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                                                                                                                                                                                       <  5:666053/71‑1 (MQ=255)
                                                                                                                                                                                                                                                tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAagag                                                                                                                                                         <  7:838686/149‑1 (MQ=255)
                                                                                                                                                                                                                                                tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAagag                                                                                                                                                         >  8:184612/1‑149 (MQ=255)
                                                                                                                                                                                                                                                tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAagag                                                                                                                                                         <  7:184612/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAagag                                                                                                                                                         <  6:575678/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACACATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTAtt                                                                                    >  4:151926/1‑217 (MQ=255)
                                                                                                                                                                                                                                                   acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGcc                                                                                                                                                      >  8:613112/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTg                                                        >  4:384495/1‑243 (MQ=255)
                                                                                                                                                                                                                                                     gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACACATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTgg                                                                                <  2:151927/217‑1 (MQ=255)
                                                                                                                                                                                                                                                       acgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCaaa                                                                                                                                                   <  6:174263/148‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGca                                                                                                                                                >  8:501052/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGcac                                                                                                                                               >  7:643791/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGca                                                                                                                                                >  6:273478/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           aagaccgggttttCGGGGTTAAGCCCGAAAGAAAACCGGGAAATCCGGCCGATCGATTAACCGTTTGGGTATCAAATGAAGCGTGACTTTCCGCATCGGACAAGGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGTATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTgcg                <  4:9851/265‑1 (MQ=255)
                                                                                                                                                                                                                                                           aaaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGcaca                                                                                                                                              >  8:286785/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACAcc                                                                                                                                            >  5:316368/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCTTGAAAACCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCAt          <  2:552501/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCTGCCGCTCGATTATCCGCGTGTGTATCAACGGAAGCGTGACATTCCGCATCTGACAATGTCGCTGAACGGGGGTATAAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGCTGGTCGGGCGCGAGGCGTATACGAATCCGGGTATTCTGGCGGCGGTAGACGGGGAGAATCTTGTTTCCCCGGATACCGTTGCCGATCCGTTGGCGGTTGTgcgctgcctg         >  4:595338/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa                                                                                                                                       <  5:242666/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAAc                                                                                                                                      >  8:220933/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACat                                                                                                                                    >  6:720198/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      ttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtgg                                                                                                                                                                            >  6:745616/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                       tCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTgg                                                                                <  4:206026/199‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                                                                                                                 <  8:706808/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCTGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTTGGATTACGATGCCGATCCGGTGGTGGTCGTGCGCTCCATTTATCCGt  >  2:23479/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                            gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggatg                                                                                                                             <  7:88307/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                               
ATGATGAGATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGT  >  NC_000913/4261982‑4262521

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: