Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,937 T→G 5.6% intergenic (+22/+93) appA → / ← etk phosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,9370TG 5.6% 873.8 / 33.3 248intergenic (+22/+93)appA/etkphosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base T (114/120);  new base G (4/10);  total (118/130)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

Only 100 of 251 total aligned reads displayed.

GCCGCCAGGTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCC  >  NC_000913/1041663‑1042191
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                              
gccgccAGGTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGt                                                                                                                                                                                                                                                                <  2:239033/275‑1 (MQ=255)
                   gtgtTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAANCTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCNTGAATGAAGCCCGCATACCGGCGTGCAGTTTGTAATGCAAAAAAANANACATTACATTACCGCAATGCTCTTAGGCtgat                                                                                                                                                                                                                                         >  3:77260/1‑279 (MQ=255)
                                                        aGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACga                                                                                                                                                                                                                                >  4:327844/1‑251 (MQ=255)
                                                          tGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAACTGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                        >  1:180943/1‑225 (MQ=255)
                                                                           ggTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGAATTCAGTTACCTGAATGCTCTGGGGCTGATGACAAACGAAGAACTGTCAAATGGTTAGCCCGGAAAAGGCCTTCAcgccgaa                                                                                                                                                                                     >  4:170737/1‑273 (MQ=255)
                                                                                         tACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAACGAGCATTCAGTTACCTGAATGATCTGAGGCTGATGACAAACGAAGAACTGTCTACTGCGTAGACCGGAAAAGGCGTTACCGCCGCATCCGGCCACTTTc                                                                                                                                                                        >  4:522961/1‑274 (MQ=255)
                                                                                             gcagATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                                                                                                                                                   <  2:114707/275‑1 (MQ=255)
                                                                                                 aTGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAcc                                                                                                                                                                                                                                                            >  4:589691/1‑182 (MQ=255)
                                                                                                    cGTGATAAAACGCCGCTTTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAct                                                                                                                                                               <  2:202796/272‑1 (MQ=255)
                                                                                                        aTAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTACAGCCGCATCCGGCCACTTTCAggtctgtcccttat                                                                                                                                                         >  4:78562/1‑261 (MQ=255)
                                                                                                         tAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                                                                                                                                     <  1:85560/277‑1 (MQ=255)
                                                                                                                ccgcTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTATGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGt                                                                                                                                                   <  2:710516/272‑1 (MQ=255)
                                                                                                                   ctataaTTAAATACGNCGCCCGGGGAGGTGAAACTGACCCTGGCAGGNTGTGAAGANNGAAATGCGCAGGGCATGTGTTCGNTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                                                                                                                                      <  3:62218/181‑1 (MQ=255)
                                                                                                                   cTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAcc                                                                                                                                                                                                                                                            >  2:34836/1‑164 (MQ=255)
                                                                                                                   cTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                                                                                                                                      >  1:62218/1‑186 (MQ=255)
                                                                                                                      taaTTAAATACGCCGCCCGGGGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTa                                                                                                                                      <  1:77260/277‑1 (MQ=255)
                                                                                                                              tACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGt                                                                                                                                                                                                                                                                >  5:1688764/1‑149 (MQ=255)
                                                                                                                                 gccgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                                                                                                                             <  6:661444/149‑1 (MQ=255)
                                                                                                                                  ccgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTActct                                                                                                                                                             >  2:541873/1‑244 (MQ=255)
                                                                                                                                   cgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACNt                                                                                                                                                                                                                                                           <  7:102623/149‑1 (MQ=255)
                                                                                                                                       cGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCATGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAAt                                                                                                                                                                                                                                                       <  5:1116926/149‑1 (MQ=255)
                                                                                                                                       cGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAAGGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACTCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGGAATAGTTATAGCCg                                                                                                                         <  2:263342/275‑1 (MQ=255)
                                                                                                                                        ggAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATg                                                                                                                                                                                                                                                      >  5:967007/1‑149 (MQ=255)
                                                                                                                                           gagGTGAAACTGAACCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                                                                                                                                                   <  7:232145/149‑1 (MQ=255)
                                                                                                                                            agGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctct                                                                                                                                                                                                                                                  <  6:908920/149‑1 (MQ=255)
                                                                                                                                              gTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGa                                                                                                                                                                                                                                                >  6:1152255/1‑149 (MQ=255)
                                                                                                                                               tGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAg                                                                                                                                                                                                                                               >  6:1041907/1‑149 (MQ=255)
                                                                                                                                                gAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAgg                                                                                                                                                                                                                                              >  6:607646/1‑149 (MQ=255)
                                                                                                                                                gAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAgg                                                                                                                                                                                                                                              >  6:667096/1‑149 (MQ=255)
                                                                                                                                                    cTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGGGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtga                                                                                                                                                                                                                                          <  6:251480/149‑1 (MQ=255)
                                                                                                                                                     tGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgat                                                                                                                                                                                                                                         <  5:425591/149‑1 (MQ=255)
                                                                                                                                                       aCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgatga                                                                                                                                                                                                                                       <  7:631391/149‑1 (MQ=255)
                                                                                                                                                       aCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGCAAAGGCGTTCACGCCGCATCCGGCCCCTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTCACTGTAAGCGGTGc                                                                                                         >  2:257187/1‑275 (MQ=255)
                                                                                                                                                        cccTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                                                                                                                                      >  5:1383576/1‑149 (MQ=255)
                                                                                                                                                            ggCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGa                                                                                                                                                                                                                                                <  5:731898/135‑1 (MQ=255)
                                                                                                                                                                 gATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                                                                                                                             <  8:133069/149‑1 (MQ=255)
                                                                                                                                                                  aTGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaagaa                                                                                                                                                                                                                            <  5:768908/149‑1 (MQ=255)
                                                                                                                                                                        agtgCGAAAGGCGCATGGCATGTTTTCGTTGGAAGGTTTTAGGCAAAGCGTGAATGAAGCACGCATACCGGCTTGCAGTCTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTAACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATcaca                                                                                        <  4:257184/272‑1 (MQ=255)
                                                                                                                                                                        agagCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCt                                                                                                                                                                                                                      <  6:1688764/149‑1 (MQ=255)
                                                                                                                                                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                                                                                                                                                   <  8:712121/117‑1 (MQ=255)
                                                                                                                                                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                                                                                                                                      <  1:104652/130‑1 (MQ=255)
                                                                                                                                                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                                                                                                                                      >  3:104652/1‑130 (MQ=255)
                                                                                                                                                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGGGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAAt                                                                                                                                                                                                                   >  7:357568/1‑149 (MQ=255)
                                                                                                                                                                                    gcAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                                                                                                                                                   >  2:671057/1‑188 (MQ=255)
                                                                                                                                                                                      aGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCgg                                                                                                                                                                                                        >  6:293083/1‑149 (MQ=255)
                                                                                                                                                                                        ggCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAct                                                                                                                                                               <  4:671052/188‑1 (MQ=255)
                                                                                                                                                                                        ggCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATACCCGTAACTGTAAGCGGTTCTGGCCCGTTTAATACCACCATTTAGGATAgcg                                                                         >  4:529759/1‑274 (MQ=255)
                                                                                                                                                                                          cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGaaaa                                                                                                                                                                                                    <  5:311725/149‑1 (MQ=255)
                                                                                                                                                                                            tgtgTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAgg                                                                                                                                                                                                  <  7:950374/149‑1 (MQ=255)
                                                                                                                                                                                              tgtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCg                                                                                                                                                                                                <  7:1509796/149‑1 (MQ=255)
                                                                                                                                                                                               gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                                                                                                                               <  5:219816/149‑1 (MQ=255)
                                                                                                                                                                                               gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                                                                                                                               <  5:1560027/149‑1 (MQ=255)
                                                                                                                                                                                                ttCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                                                                                                                              >  6:575917/1‑149 (MQ=255)
                                                                                                                                                                                                ttCGTCGGCAGGTTTTACGCAAATCGCGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                                                                                                                              >  8:1440973/1‑149 (MQ=255)
                                                                                                                                                                                                  cGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACttt                                                                                                                                                                         <  1:285889/168‑1 (MQ=255)
                                                                                                                                                                                                    ttGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCTTAACTGTAAGCGGTTCTGGCGCGTTTAATCACCACATTGAGGATAGCCCCTTTAATATTg                                                             >  4:435570/1‑274 (MQ=255)
                                                                                                                                                                                                    ttGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATcaca                                                                                        <  2:520984/247‑1 (MQ=255)
                                                                                                                                                                                                           gTTTTACGCAAATCGGGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcctgc                                                     <  2:435574/275‑1 (MQ=255)
                                                                                                                                                                                                               tACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                               <  6:1254974/149‑1 (MQ=255)
                                                                                                                                                                                                                 cGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCa                                                                                                                                                                             >  7:803489/1‑149 (MQ=255)
                                                                                                                                                                                                                  gCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCAc                                                                                                                                                                            >  8:1256988/1‑149 (MQ=255)
                                                                                                                                                                                                                  gCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCAc                                                                                                                                                                            >  8:1396814/1‑149 (MQ=255)
                                                                                                                                                                                                                   cAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACt                                                                                                                                                                           >  8:1152701/1‑149 (MQ=255)
                                                                                                                                                                                                                      aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTc                                                                                                                                                                        >  7:808748/1‑149 (MQ=255)
                                                                                                                                                                                                                      aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                                                                                                                                                   >  2:509539/1‑154 (MQ=255)
                                                                                                                                                                                                                      aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTTGGAGTAACTATAACAGTAATAGTTATCGCCGTAACTGTAAAGGGTGCTGGCGCGTTTAATCACCCCATTGGGGATAGCGCATTAAATTCTGACGCCTACCTTTTCCAgcc                                          >  4:558876/1‑273 (MQ=255)
                                                                                                                                                                                                                       tCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCa                                                                                                                                                                       <  5:1316281/149‑1 (MQ=255)
                                                                                                                                                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                                                                                                                                    <  6:1535364/149‑1 (MQ=255)
                                                                                                                                                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                                                                                                                                    >  5:109929/1‑149 (MQ=255)
                                                                                                                                                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                                                                                                                                    <  8:117628/149‑1 (MQ=255)
                                                                                                                                                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAct                                                                                                                                                               <  4:509535/154‑1 (MQ=255)
                                                                                                                                                                                                                           gaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTg                                      <  2:558880/274‑1 (MQ=255)
                                                                                                                                                                                                                            aatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                                                                                                                               <  8:308137/120‑1 (MQ=255)
                                                                                                                                                                                                                             atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTa                                                                                                                                                                 >  6:1362234/1‑149 (MQ=255)
                                                                                                                                                                                                                              tgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAc                                                                                                                                                                <  6:348186/149‑1 (MQ=255)
                                                                                                                                                                                                                                  gCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                                                                                                                             <  7:1173607/84‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                                                                                                                                          >  6:766646/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aTACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCgg                                                                                                                                                      >  7:317868/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aTACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCg                                                                                                                                                       <  7:1245126/148‑1 (MQ=255)
                                                                                                                                                                                                                                         tACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                                                                                                                                     <  7:696539/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ggtatgCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                                                                                                                                     <  7:1402410/144‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTaa                                                                                                                                                 <  5:575917/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCg                                                                                                                                        >  6:182661/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        tGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGTTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtgt         >  2:317921/1‑274 (MQ=255)
                                                                                                                                                                                                                                                         gTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTaa                                                                                                                                     <  5:626825/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTaa                                                                                                                                     <  7:493958/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           aaTGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCGCGGAGTAACTATAACCGTAATa                                                                                                                                   <  5:1041907/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAg                                                                                                                                  >  8:960434/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTaa                                                                                                                                                                                                                    <  7:536268/61‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTaa                                                                                                                                                                                                                    >  8:536268/1‑61 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAg                                                                                                                            >  5:1433725/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGc                                                                                                                           >  7:1336080/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtct                                                                                                                                                         <  4:260303/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtct                                                                                                                                                         <  2:752432/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGcc                                                                                                                          <  7:899534/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                                                                                                                                          >  7:1569910/1‑115 (MQ=255)
                                                                                                                                                                                                                                                                      aaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                                                                                                                                          <  8:1569910/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                       agagCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                                                                                                                                     >  5:1265171/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                           cATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGttt                                                                                               >  4:717353/1‑169 (MQ=255)
                                                                                                                                                                                                                                                                           cATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                                                                                                                               >  6:1379687/1‑73 (MQ=39)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                              
GCCGCCAGGTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCC  >  NC_000913/1041663‑1042191

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