Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 5.1% I27M (ATT→ATG)  yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5190AC 5.1% 406.1 / 11.5 118I27M (ATT→ATG) yadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (56/56);  new base C (3/3);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 120 total aligned reads displayed.

AGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGC  >  NC_000913/151263‑151778
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                    
agCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATACCTTTCAGGAGTAACCTTACCTGAAAAATATCCCAACCCAACGACGCCAAGATAAATATTACCCTGGTTACTGATGCCGCTTTTGCAGGTGCTGTTTTTTACTTTGGCAACAAGTTCAATATCCTGTCCAGCAATGCC‑TGGACAGGCAAta                                                                                                                                                                                                                                                   >  1:101482/1‑274 (MQ=255)
    gtataagagacatccGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCaaa                                                                                                                                                                                                                                               <  1:18235/262‑1 (MQ=255)
     gcagcCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              >  1:144647/1‑275 (MQ=255)
                            tctAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGCTAACTTTCAGGAGTAACATTACCTGAAAAATATCACACCCCAACGACGCCAAGATCCATATTACCCTGGTTACTGCTTCCGCTTTTACAGGTGCTGTTTTTTACATTATCAACAAGGTCAATATCCTGTCCAGCAATGGCGCGGAC‑GGGGATGCAAAACAGGAATTAAAGAAACAGAgtt                                                                                                                                                                                                                        >  4:17856/1‑274 (MQ=255)
                            tctAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCCGGTGCTGTTTTTTACATTTGAAACAAGGTCAATTTCCTTTCCCGCAATGCC‑TGGACAGGTGGNGCCAANNAATAATAaccattttctttta                                                                                                                                                                                                                        >  3:18235/1‑262 (MQ=255)
                                     aaattttaacctaTTTAGCTAATCGCCATNNCCCTTTATTGCACAGTCGCATTCATTGATANTCACCTCTTTTCCACCTGGTTAATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGTCGCCANGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGGTGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                              <  3:68943/263‑1 (MQ=255)
                                                            cGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGCATTATTCCCTTCAACAGCTTACTTTGAGGGGGATGc                                                                                                                                                                                       >  4:490575/1‑275 (MQ=255)
                                                                         tAATGCAAAGTCGGTTAACTTTATAGTGANNNCTTTTCCACCTGGTTCAAAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCCAAGACGCCAAGATAAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCAGTCCATCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACTTTGAGTGGGATGCATATATAGTTCCt                                                                                                                                                                          <  3:144647/275‑1 (MQ=255)
                                                                                 aGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCACTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAGGGTCAATATCCTGTCCAGCAATGGC‑TGTACAGGCGTTGCAAAACAGTAATGACACACAAAGAATTATCCTCTTCATAAGCTTACGTTGTGTGGGCTACATATCTTGTTCTTGActcttt                                                                                                                                                                  >  4:72181/1‑273 (MQ=255)
                                                                                 aGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAAAAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  2:136155/1‑247 (MQ=255)
                                                                                        tACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCACTGGC‑TGGACAGGCGATGCAAAACAGTACTGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGGGTGGGATGCCTATATAGTTCCAGAATCTTTAAGTCTg                                                                                                                                                           >  4:446902/1‑275 (MQ=255)
                                                                                        tACAGTGATAGTGAACTCGTATCCACCTGGTTGATAACGTTCAGGAGTAACATTACCTGAAAAATATCCCAGCCCAACGAAGCCAAGATCAATATTACCCGGGTTACTGATTCCTCTTTTACAGGTGCTGTTATTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGAGTCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTg                                                                                                                                                           <  1:40814/275‑1 (MQ=255)
                                                                                         aCAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                          <  2:108294/275‑1 (MQ=255)
                                                                                          cAGGGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtt                                                                                                                                                         <  2:17856/275‑1 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGTCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                                                                                                                                   >  2:518258/1‑198 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  2:225491/1‑235 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTCTTGCAGGTGCTGTTCTTTACATTGGNAACAAGGCCAATATCCTGTTCCGCAATGGC‑TGGACAGGCGATGCAAAACCGTAATGAAAGAAACAGAAAAATACGCTCAANCAGACTACCATGAGTGGGAGGCAAANAAAGTTCTCGAATCTGGAAGtaggtgtgt                                                                                                                                                      >  3:84330/1‑267 (MQ=255)
                                                                                               aTAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtctc                                                                                                                                                                                                                  >  4:153726/1‑213 (MQ=255)
                                                                                                tAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtat                                                                                                                                                   >  2:255631/1‑275 (MQ=255)
                                                                                                 aGTNAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTANCATTNCCTGAAAAATAGCCCACCCCAACGACGTCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                               <  4:518253/198‑1 (MQ=255)
                                                                                                 aGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagtta                                                                                                                                                    <  1:84330/273‑1 (MQ=255)
                                                                                                       ctctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTATGTTAGTTATAGGTTAc                                                                                                                                            >  2:68705/1‑275 (MQ=255)
                                                                                                             tCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCaa                                                                                                                                                                                                                                                                     <  7:4659/149‑1 (MQ=255)
                                                                                                              ccACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAAt                                                                                                                                                                                                                                                                    >  7:308217/1‑149 (MQ=255)
                                                                                                                 ccTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                                                                                                                                                                                           >  4:81586/1‑186 (MQ=255)
                                                                                                                   tGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑Tg                                                                                                                                                                                                                                                              <  5:333711/149‑1 (MQ=255)
                                                                                                                   tGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑Tg                                                                                                                                                                                                                                                              <  5:398645/149‑1 (MQ=255)
                                                                                                                    ggTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTTCTGTTTTTTACATCGGCAACACGGTCAATATCCTGTCCAGCAATGGC‑TTGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  2:403449/1‑212 (MQ=255)
                                                                                                                    ggTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                                                                                                                                <  2:490580/274‑1 (MQ=255)
                                                                                                                     gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGa                                                                                                                                                                                                                                                            <  6:68080/149‑1 (MQ=255)
                                                                                                                     gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                                                                                                                                       <  2:81587/186‑1 (MQ=255)
                                                                                                                        gcaagctgtCAGGCGTAACATTACCTGAACAAAACCCCACCCCAACGACGCCAAGATCACTATTACCCTGGTTACTGATTCCGTTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGAAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg                                                                                                                           <  4:85752/267‑1 (MQ=255)
                                                                                                                        gATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAg                                                                                                                                                                                                                                                         >  6:66667/1‑149 (MQ=255)
                                                                                                                            aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGa                                                                                                                                                                                                                                                     <  5:66667/149‑1 (MQ=255)
                                                                                                                              tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATg                                                                                                                                                                                                                                                   >  8:454290/1‑149 (MQ=255)
                                                                                                                               ttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa                                                                                                                     >  2:70206/1‑274 (MQ=255)
                                                                                                                                   ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              >  6:52076/1‑149 (MQ=255)
                                                                                                                                        aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGa                                                                                                                                                                                                                                                            >  4:327820/1‑130 (MQ=255)
                                                                                                                                        aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt                                                                                                                               >  2:85753/1‑255 (MQ=255)
                                                                                                                                            ttACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAgg                                                                                                                                                                                                                                                        <  2:327824/130‑1 (MQ=255)
                                                                                                                                                 tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGacagaa                                                                                                                                                                                                                                >  6:103965/1‑149 (MQ=255)
                                                                                                                                                  gAAAAATATCCCACCCCAACGAAGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACAATGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                               <  5:275425/149‑1 (MQ=255)
                                                                                                                                                      aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              >  8:213485/1‑130 (MQ=255)
                                                                                                                                                      aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              <  7:213485/130‑1 (MQ=255)
                                                                                                                                                       atatCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                                                                                                                                          >  6:336844/1‑149 (MQ=255)
                                                                                                                                                          tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                                                                                                                                       <  6:266755/149‑1 (MQ=255)
                                                                                                                                                           cccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCACTGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                                                                                                                                                                                                      >  8:243942/1‑149 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                                                                                                                                     <  5:52076/149‑1 (MQ=255)
                                                                                                                                                                           cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTa                                                                                                                                                                                                      >  5:330774/1‑149 (MQ=255)
                                                                                                                                                                               aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                                                                                                                                                                                  <  5:455460/149‑1 (MQ=255)
                                                                                                                                                                                tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTg                                                                                                                                                                                                 >  8:148413/1‑149 (MQ=255)
                                                                                                                                                                                    tatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                                                                                                                                                                                >  4:529703/1‑146 (MQ=255)
                                                                                                                                                                                     tttccccTGGTTACTTATTCCGCTTTTGCAGGGTCTTTTATTTACATGGGCAACAAGGTCAAAATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGGAATGACAGAACCAGAATTATTCTCTTCATCAGCTTACTGTGAGTGGGATGCATATATAGTTACTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTg                                                               <  4:78314/270‑1 (MQ=255)
                                                                                                                                                                                     atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                                                                                                                            >  6:365767/1‑149 (MQ=255)
                                                                                                                                                                                      ttACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                                                                                                                                                                                           >  7:340389/1‑149 (MQ=255)
                                                                                                                                                                                      ttACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                                                                                                                                                                                           >  5:162021/1‑149 (MQ=255)
                                                                                                                                                                                       tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          <  8:308217/149‑1 (MQ=255)
                                                                                                                                                                                        aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                                                                                                                            <  2:529708/146‑1 (MQ=255)
                                                                                                                                                                                              gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                                                                                                                                      >  4:100541/1‑146 (MQ=255)
                                                                                                                                                                                                  cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGTGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata                                                                                                                                                                                  <  2:100542/146‑1 (MQ=255)
                                                                                                                                                                                                    gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATAGATAGTTCCTGAATCTGTAAGTCAGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACCCATCGGataat                                                >  2:466700/1‑274 (MQ=255)
                                                                                                                                                                                                       tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt                                                                                                                                                                          <  6:330774/149‑1 (MQ=255)
                                                                                                                                                                                                            ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATc                                                                                                                                                                     >  5:97472/1‑149 (MQ=255)
                                                                                                                                                                                                               tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGt                                                                                                                                                                  <  7:454290/149‑1 (MQ=255)
                                                                                                                                                                                                                tcAGGGGGTGGTTTTTACATTGGAAAAAAGGGCAATATACTGTCAAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCAtt                                   <  4:70205/274‑1 (MQ=255)
                                                                                                                                                                                                                 cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                                                                                                                                                <  7:148413/149‑1 (MQ=255)
                                                                                                                                                                                                                 cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                                                                                                                                                <  5:103965/149‑1 (MQ=255)
                                                                                                                                                                                                                       cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                          >  8:122584/1‑149 (MQ=255)
                                                                                                                                                                                                                        tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCt                                                                                                                                                            <  5:109371/146‑1 (MQ=255)
                                                                                                                                                                                                                        tGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCt                                                                                                                                                            >  6:109371/1‑146 (MQ=255)
                                                                                                                                                                                                                          ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTTAATCTGTAAGTCTgttag                                                                                                                                                       >  8:286528/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttgttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTgg                        <  2:446907/272‑1 (MQ=255)
                                                                                                                                                                                                                               tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                  <  7:56571/149‑1 (MQ=255)
                                                                                                                                                                                                                                aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                                                                                                                                 <  8:340389/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              <  8:234849/103‑1 (MQ=255)
                                                                                                                                                                                                                                 cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  7:234849/1‑103 (MQ=255)
                                                                                                                                                                                                                                    tGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAg                                                                                                                                                               >  2:523853/1‑131 (MQ=255)
                                                                                                                                                                                                                                     ggCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                  <  6:49453/143‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                  >  5:49453/1‑143 (MQ=255)
                                                                                                                                                                                                                                      gCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACATAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACt                                                                                                                                           >  8:33343/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCg                                                                                                            <  2:85779/178‑1 (MQ=255)
                                                                                                                                                                                                                                           aaGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAATTATAGGTTACTGTaaa                                                                                                                                      >  5:52355/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAt                                                                                                                                     >  8:201257/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt                                                                                                                                   >  6:259473/1‑149 (MQ=255)
                                                                                                                                                                                                                                               taaaTATCCTGTCCAGCAATGGC‑TGGAAAGGCAATGAAAAAAAGTAATGCCAGAAACAGAATTATTCCCTTAATCAGCTTACTTGGAGTGGGAGGCATATATATTTCCTGAATTTGTAAGTCTG‑TAGTTATAGGTTACTGTAAATGTgg                                                                                                                                 <  7:286528/147‑1 (MQ=255)
                                                                                                                                                                                                                                                  atcgccTGTCCAGCCATTGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          <  7:149056/86‑1 (MQ=38)
                                                                                                                                                                                                                                                  atcgccTGTCCAGCCATTGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          >  8:149056/5‑90 (MQ=38)
                                                                                                                                                                                                                                                      ccTGTCCAGCCATTGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                                                                                                                                                                             >  7:82928/1‑99 (MQ=255)
                                                                                                                                                                                                                                                       cTGTCCAGCCATTGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                                                                                                                                        >  6:78178/1‑103 (MQ=255)
                                                                                                                                                                                                                                                        tGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                                                                                                                                                                                >  1:61639/1‑78 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGACTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                                                                                                                                 <  6:388674/108‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                                                                                                                                            >  2:369301/1‑97 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                                                                                                                                 >  5:388674/1‑108 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                                                                                                                                        <  4:369296/97‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaaa                                                                                                                                      >  6:308389/1‑131 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaaa                                                                                                                                      <  5:308389/131‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGc                                                              <  1:40678/203‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCa                                     <  2:326009/228‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGGTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                                                                                   >  8:362692/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aaTGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACtt                                                                                                                 >  7:86545/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                    
AGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGC‑TGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGC  >  NC_000913/151263‑151778

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