Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,548,914 G→C 7.6% G36A (GGT→GCT)  fimH → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,9140GC 7.6% 646.4 / 35.0 184G36A (GGT→GCT) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base G (90/80);  new base C (6/8);  total (96/88)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-01

Only 100 of 185 total aligned reads displayed.

ATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCG  >  NC_000913/4548643‑4549182
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                            
atatcggCGCAACAAAAATATGTCGCGTGGATGATTCCCCACAATCGGTGAACTTCCCGTACAAGGTCGGGGATTTGACAGTACAGGGCGGAGCCACTCAGGGACCCATTCAGTAAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTgg                                                                                                                                                                                                                                                                           <  4:380567/270‑1 (MQ=255)
    tGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAgc                                                                                                                                                                                                                                                                       >  4:466453/1‑275 (MQ=255)
         cAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACTCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTTCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGGCTGTAAAAACGCCAATGGTACCGCTATCCCTATTTGCGGTgcgaccgcaca                                                                                                                                                                                                                                                                  >  2:380568/1‑265 (MQ=255)
         cAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGcc                                                                                                                                                                                                                                                                    >  2:284896/1‑273 (MQ=255)
                     ttCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAACCCGCCATTGGTACCGCTATCCCTCTTGGCGGTGGCCGCCCCAATGTTTATGTaaa                                                                                                                                                                                                                                                      >  4:421964/1‑275 (MQ=255)
                                                          gTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTc                                                                                                                                                                                                                 >  2:169606/1‑275 (MQ=255)
                                                                      agCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTc                                                                                                                                                                                                      >  2:79383/1‑274 (MQ=255)
                                                                                gTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCCATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCt                                                                                                                                                                                           >  1:620/1‑275 (MQ=255)
                                                                                           aGCCACTCAGGGAACCATTGCGGCAGTGAATAGAAGCACCTATCCCTACAGCTGAACCCGAAGCGAGGATGGGAATGAAACGAGTTATTACCCTGTTGGCTGTACTGCTGATGGGTTGGTCGGTAAATGCCTGGTCATTCGCCGGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATaa                                                                                                                                                                                 <  4:283417/274‑1 (MQ=255)
                                                                                             ccACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTACTGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCCATGGTACCGCTATCCCTATTGGCGGGGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTTGTCGTGGATCTTTCCACGCAACTCTTTTTCCATTAccg                                                                                                                                                                              >  4:5610/1‑273 (MQ=255)
                                                                                                         acATTCAGGCAGTGATTAGCATCACCTATAACTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGaaa                                                                                                                                                                   <  1:100746/273‑1 (MQ=255)
                                                                                                              cAGGAAGTGATTAGCATCACCGAACCCTAAAGCTGAACCCGAAGAGATGATTGGAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTa                                                                                                                                                             <  4:169605/275‑1 (MQ=255)
                                                                                                                 gCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCg                                                                                                                                                                       >  4:112122/1‑262 (MQ=255)
                                                                                                                     tGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTgg                                                                                                                                                                                                                                                                           <  3:134907/158‑1 (MQ=255)
                                                                                                                            aaTCAACTATACCTACAGCTAAAACCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCGGGAATGCTGAGGGGCTGGTCGGTAGATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                                                                                 <  2:25051/272‑1 (MQ=255)
                                                                                                                               cACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTTNTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCCCCCGTCGTGAATTGGGGGCAAAACCTGTTCGTGGATCTTTCGANGCAAATCTTTTGCCCTACCGAtttancg                                                                                                                                                                       >  3:96622/1‑243 (MQ=255)
                                                                                                                                 ccTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTg                                                                                                                                                                                                                                                 >  4:161865/1‑172 (MQ=255)
                                                                                                                                       ccTACAGCTGAACCCGAAGAGATGATTGTACTGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCa                                                                                                                                                                                                                                                                   <  5:733377/148‑1 (MQ=255)
                                                                                                                                         tACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAgcg                                                                                                                                                                                                                                                                      >  5:306467/1‑143 (MQ=255)
                                                                                                                                         tACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAAt                                                                                                                                                                                                                                                                 >  6:39871/1‑148 (MQ=255)
                                                                                                                                          aCAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTTCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCg                                                                                                                                                                       >  2:540534/1‑237 (MQ=255)
                                                                                                                                           cAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGTCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                                                                                  <  2:366090/273‑1 (MQ=255)
                                                                                                                                           cAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCg                                                                                                                                >  2:307797/1‑275 (MQ=255)
                                                                                                                                            aGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGtt                                                                                                                                                                                                                                                              <  5:39871/148‑1 (MQ=255)
                                                                                                                                                 aaCCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCtt                                                                                                                            >  2:100784/1‑273 (MQ=255)
                                                                                                                                                 aaCCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGGTTATGTAAACCTTGCGCCAGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTAt                                                                                                                          >  4:82632/1‑275 (MQ=255)
                                                                                                                                                                 tGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCcgtc                                                                                                                                                                                                                                         >  7:251399/1‑148 (MQ=255)
                                                                                                                                                                      tGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAg                                                                                                                                                                                                                                                                        >  6:288984/1‑112 (MQ=255)
                                                                                                                                                                      tGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGaa                                                                                                                                                                                                                                    >  7:170129/1‑148 (MQ=255)
                                                                                                                                                                      tGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCa                                                                                                                                                                                                                                                                   >  6:1031136/1‑117 (MQ=255)
                                                                                                                                                                         aaCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAt                                                                                                                                                                                                                                                                 >  8:629141/1‑116 (MQ=255)
                                                                                                                                                                         aaCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAt                                                                                                                                                                                                                                                                 <  7:629141/116‑1 (MQ=255)
                                                                                                                                                                         aaCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAt                                                                                                                                                                                                                                                                 <  5:356379/116‑1 (MQ=255)
                                                                                                                                                                         aaCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAt                                                                                                                                                                                                                                                                 >  6:356379/1‑116 (MQ=255)
                                                                                                                                                                          aCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTaa                                                                                                                                                                                                                                                       <  8:511931/125‑1 (MQ=255)
                                                                                                                                                                           cGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAAc                                                                                                                                                                                                                                                     >  7:1019923/1‑126 (MQ=255)
                                                                                                                                                                           cGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAAc                                                                                                                                                                                                                                                     <  8:1019923/126‑1 (MQ=255)
                                                                                                                                                                           cGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTgg                                                                                                                                                                                                                               <  8:251399/148‑1 (MQ=255)
                                                                                                                                                                            nggTTATTACCCTGTTGGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTNTTTCNCCTGTAAAACCGCCAATGGTACCNCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGaaa                                                                                                                                                                   <  4:540533/205‑1 (MQ=255)
                                                                                                                                                                             tgTTTTTTCCCTGGTTGCTTTACTGCTGATGGACTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATACGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                              <  4:307796/274‑1 (MQ=255)
                                                                                                                                                                              tttatatcccTGTTGGCTGTCTTGCTGATGGGCTGGTCGGTAAATGCCTGGTCAGTCGCCTGTAAAACCGCCAATGGTACCGTGATCCCTATTGGCGGTGGCAGCGACAATGTTTATGGAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGTAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                              <  4:383612/267‑1 (MQ=255)
                                                                                                                                                                                 attaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCg                                                                                           <  2:317873/274‑1 (MQ=255)
                                                                                                                                                                                  ttaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTa                                                                                                                                                                                                                                                            <  5:243792/112‑1 (MQ=255)
                                                                                                                                                                                  ttaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTa                                                                                                                                                                                                                                                            >  6:243792/1‑112 (MQ=255)
                                                                                                                                                                                  ttaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAAtgtg                                                                                                                                                                                                                                <  8:598703/140‑1 (MQ=255)
                                                                                                                                                                                  ttaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAAtgtg                                                                                                                                                                                                                                >  7:598703/1‑140 (MQ=255)
                                                                                                                                                                                    aCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAt                                                                                                                                                                                                                                                                 >  5:686030/1‑105 (MQ=255)
                                                                                                                                                                                     cccTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTCTATGTAAACCCTGCGCCCGTCGTGAATGTGGGGCAAAACCt                                                                                                                                                                                                                     >  6:61474/1‑148 (MQ=255)
                                                                                                                                                                                      ccTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTg                                                                                                                                                                                                                    <  8:192900/148‑1 (MQ=255)
                                                                                                                                                                                         gTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCcgtcg                                                                                                                                                                                                                                        >  5:214642/1‑125 (MQ=255)
                                                                                                                                                                                         gTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCcgtcg                                                                                                                                                                                                                                        <  6:214642/125‑1 (MQ=255)
                                                                                                                                                                                          tttGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                  >  1:122731/1‑274 (MQ=255)
                                                                                                                                                                                           ttGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGa                                                                                                                                                                                                                                     >  5:31467/1‑126 (MQ=255)
                                                                                                                                                                                           ttGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCg                                                                                                                                                                                                                <  7:222939/147‑1 (MQ=255)
                                                                                                                                                                                            tGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCATCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCg                                                                                                                                                                       >  2:444338/1‑187 (MQ=255)
                                                                                                                                                                                              cTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCCTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGa                                                                                                                                                                                                            <  5:61474/148‑1 (MQ=255)
                                                                                                                                                                                               tGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTCTGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGCCCGTCAAATATAGTGGc                                                                             >  2:148881/1‑274 (MQ=255)
                                                                                                                                                                                                gTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATc                                                                                                                                                                                                          >  7:300246/1‑148 (MQ=255)
                                                                                                                                                                                                gTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAta                                                                                                                                                                                                                                                                <  7:658826/94‑2 (MQ=255)
                                                                                                                                                                                                gTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGCTGCCACCGCCAAta                                                                                                                                                                                                                                                                >  8:658826/1‑93 (MQ=255)
                                                                                                                                                                                                  aCTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAAtgtg                                                                                                                                                                                                                                <  5:826473/124‑1 (MQ=255)
                                                                                                                                                                                                  aCTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAAtgtg                                                                                                                                                                                                                                >  6:826473/1‑124 (MQ=255)
                                                                                                                                                                                                      ctgATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTATCTAATTTTTCCGGGACCGTaa                                                                                        >  2:195116/1‑255 (MQ=255)
                                                                                                                                                                                                        gATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCAt                                                                                                                                                                                   <  4:247283/163‑1 (MQ=255)
                                                                                                                                                                                                         aTGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAACTATAGTGGCAGTAGCTATcc                                                                  >  2:493183/1‑275 (MQ=255)
                                                                                                                                                                                                           gggCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGTGACCGTAAAATATAGTGGCAGTTGCTATCCAt                                                                >  2:132863/1‑275 (MQ=255)
                                                                                                                                                                                                              cTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCt                                                                                                                                                                                                                     <  7:572369/123‑1 (MQ=255)
                                                                                                                                                                                                                ggTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCtt                                                                                                                                                                                          <  8:300246/148‑1 (MQ=255)
                                                                                                                                                                                                                     gTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTaccac                                                      <  2:421965/275‑1 (MQ=255)
                                                                                                                                                                                                                          tGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGtt                                                                                                                                                                                                                                                              >  6:803969/1‑70 (MQ=255)
                                                                                                                                                                                                                          tGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGtt                                                                                                                                                                                                                                                              <  5:803969/70‑1 (MQ=255)
                                                                                                                                                                                                                           gCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACg                                                                                                                                                                               <  6:306872/148‑1 (MQ=255)
                                                                                                                                                                                                                            ccTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCtt                                                                                                                                                                                                        <  6:717180/122‑1 (MQ=255)
                                                                                                                                                                                                                              ttgTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCTCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGa                                                                                                                                                                              <  7:1112305/144‑1 (MQ=255)
                                                                                                                                                                                                                              tGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCaa                                                                                                                                                                                                                                                                  <  5:1065718/62‑1 (MQ=255)
                                                                                                                                                                                                                              tGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGa                                                                                                                                                                              >  8:1112305/1‑146 (MQ=255)
                                                                                                                                                                                                                                gTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCGATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTgg                                                                              >  2:341252/1‑240 (MQ=255)
                                                                                                                                                                                                                                    ttCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCGATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGt                                                                          <  4:341251/240‑1 (MQ=255)
                                                                                                                                                                                                                                     tCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTTATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTAACCATTTCCTACCACAAACAAAACGCcgcgc                                      >  4:260408/1‑275 (MQ=255)
                                                                                                                                                                                                                                      ggCCGGTAAAACCGCCAATGGTACCGCTATCCCAATTGGCGGTGGCAGCGCAAATGTTTACGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                     <  2:476292/274‑1 (MQ=255)
                                                                                                                                                                                                                                      cGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                     <  2:5610/275‑1 (MQ=255)
                                                                                                                                                                                                                                        ccTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGa                                                                                                                                                                     >  7:715796/1‑145 (MQ=255)
                                                                                                                                                                                                                                              aaaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTAc                                                                                                                                                            <  8:726802/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTAGTGAATGTGGGGCAAAACCTGGTCGTGGATCATTCGACGCAAATCTTTTTCCATAAAGATTATCCGGAAACCATTACAGa                                                                                                                                                         >  7:914293/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        tGGTACCGCTATCCCTATTGGCGGTGGAAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                   <  2:25570/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCa                                                                                                                                                                >  6:657536/1‑129 (MQ=255)
                                                                                                                                                                                                                                                              cGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCa                                                                                                                                                                >  7:371713/1‑128 (MQ=255)
                                                                                                                                                                                                                                                              cGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCa                                                                                                                                                                <  8:371713/128‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGc                                                                                                                                           >  8:708490/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     ccTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTAc                                                                                                                                                            <  5:716034/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATg                                                                                                                                                    <  5:108736/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATg                                                                                                                                                    >  6:108736/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                         ttGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCAtt                                                                                                                                                              <  6:714654/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACg                                                                                                                                                                                                  >  5:1067006/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                          tGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACg                                                                                                                                                                                                  <  6:1067006/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATa                                                                                                                                                                                  >  6:404188/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                          tGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATa                                                                                                                                                                                  <  5:404188/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCg                                                                                                                                                                                                                >  7:44204/1‑66 (MQ=255)
                                                                                                                                                                                                                                                                            gCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCg                                                                                                                                                                                                                <  8:44204/66‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCt                                                                                                                             >  7:781384/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                            
ATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCG  >  NC_000913/4548643‑4549182

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