Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 5.1% A26A (GCC→GCA yadK ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5220GT 5.1% 408.4 / 9.6 119A26A (GCC→GCAyadKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (56/57);  new base T (3/3);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 119 total aligned reads displayed.

AGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGC  >  NC_000913/151263‑151778
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                 
agCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATACCTTTCAGGAGTAACCTTACCTGAAAAATATCCCAACCCAACGACGCCAAGATAAATATTACCCTGGTTACTGATGCCGCTTTTGCAGGTGCTGTTTTTTACTTTGGCAACAAGTTCAATATCCTGTCCAGCAATGCCTGGACAGGCAAta                                                                                                                                                                                                                                                   >  1:101482/1‑274 (MQ=255)
    gtataagagacatccGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                                                                                                                                                                                                               <  1:18235/262‑1 (MQ=255)
     gcagcCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              >  1:144647/1‑275 (MQ=255)
           aTGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTTCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTat                                                                                                                                                                                                                                        >  1:44866/1‑274 (MQ=255)
                            tctAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGCTAACTTTCAGGAGTAACATTACCTGAAAAATATCACACCCCAACGACGCCAAGATCCATATTACCCTGGTTACTGCTTCCGCTTTTACAGGTGCTGTTTTTTACATTATCAACAAGGTCAATATCCTGTCCAGCAATGGCGCGGACGGGGATGCAAAACAGGAATTAAAGAAACAGAgtt                                                                                                                                                                                                                        >  4:17856/1‑274 (MQ=255)
                            tctAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTTATTCCGCTTTTGCCGGTGCTGTTTTTTACATTTGAAACAAGGTCAATTTCCTTTCCCGCAATGCCTGGACAGGTGGNGCCAANNAATAATAaccattttctttta                                                                                                                                                                                                                        >  3:18235/1‑262 (MQ=255)
                                     aaattttaacctaTTTAGCTAATCGCCATNNCCCTTTATTGCACAGTCGCATTCATTGATANTCACCTCTTTTCCACCTGGTTAATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGTCGCCANGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGGTGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                              <  3:68943/263‑1 (MQ=255)
                                                         caTCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          <  2:499927/275‑1 (MQ=255)
                                                            cGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGCATTATTCCCTTCAACAGCTTACTTTGAGGGGGATGc                                                                                                                                                                                       >  4:490575/1‑275 (MQ=255)
                                                                         tAATGCAAAGTCGGTTAACTTTATAGTGANNNCTTTTCCACCTGGTTCAAAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCCAAGACGCCAAGATAAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCAGTCCATCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACTTTGAGTGGGATGCATATATAGTTCCt                                                                                                                                                                          <  3:144647/275‑1 (MQ=255)
                                                                                 aGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAAAAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  2:136155/1‑247 (MQ=255)
                                                                                        tACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCACTGGCTGGACAGGCGATGCAAAACAGTACTGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGGGTGGGATGCCTATATAGTTCCAGAATCTTTAAGTCTg                                                                                                                                                           >  4:446902/1‑275 (MQ=255)
                                                                                        tACAGTGATAGTGAACTCGTATCCACCTGGTTGATAACGTTCAGGAGTAACATTACCTGAAAAATATCCCAGCCCAACGAAGCCAAGATCAATATTACCCGGGTTACTGATTCCTCTTTTACAGGTGCTGTTATTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAGTCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTg                                                                                                                                                           <  1:40814/275‑1 (MQ=255)
                                                                                         aCAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                          <  2:108294/275‑1 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGTCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                                                                                                                                   >  2:518258/1‑198 (MQ=255)
                                                                                             tGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  2:225491/1‑235 (MQ=255)
                                                                                               aTAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtctc                                                                                                                                                                                                                  >  4:153726/1‑213 (MQ=255)
                                                                                                tAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtat                                                                                                                                                   >  2:255631/1‑275 (MQ=255)
                                                                                                 aGTNAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTANCATTNCCTGAAAAATAGCCCACCCCAACGACGTCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                               <  4:518253/198‑1 (MQ=255)
                                                                                                 aGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagtta                                                                                                                                                    <  1:84330/273‑1 (MQ=255)
                                                                                                   tGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                                                                                                                                              <  2:153728/213‑1 (MQ=255)
                                                                                                       ctctTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTATGTTAGTTATAGGTTAc                                                                                                                                            >  2:68705/1‑275 (MQ=255)
                                                                                                                 ccTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                                                                                                                                                                                                                           >  4:81586/1‑186 (MQ=255)
                                                                                                                   tGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                                                                                                                                                                                              <  5:398645/149‑1 (MQ=255)
                                                                                                                    ggTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTTCTGTTTTTTACATCGGCAACACGGTCAATATCCTGTCCAGCAATGGCTTGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  2:403449/1‑212 (MQ=255)
                                                                                                                    ggTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGc                                                                                                                                <  2:490580/274‑1 (MQ=255)
                                                                                                                     gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                                                                                                                                                                                            <  6:68080/149‑1 (MQ=255)
                                                                                                                     gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                                                                                                                                       <  2:81587/186‑1 (MQ=255)
                                                                                                                        gcaagctgtCAGGCGTAACATTACCTGAACAAAACCCCACCCCAACGACGCCAAGATCACTATTACCCTGGTTACTGATTCCGTTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTACAGAAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTg                                                                                                                           <  4:85752/267‑1 (MQ=255)
                                                                                                                         aTAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGTACAgg                                                                                                                                                                                                                                                        >  6:455460/1‑149 (MQ=255)
                                                                                                                            aCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                                                                                                                                                                                                                     <  5:66667/149‑1 (MQ=255)
                                                                                                                              tttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                                                                                                                                                                                                                   >  8:454290/1‑149 (MQ=255)
                                                                                                                               ttCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa                                                                                                                     >  2:70206/1‑274 (MQ=255)
                                                                                                                                   ggAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              >  6:52076/1‑149 (MQ=255)
                                                                                                                                        aaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                                                                                                                                                                                                                            >  4:327820/1‑130 (MQ=255)
                                                                                                                                            ttACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAgg                                                                                                                                                                                                                                                        <  2:327824/130‑1 (MQ=255)
                                                                                                                                              aCCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGaca                                                                                                                                                                                                                                   <  5:28187/149‑1 (MQ=255)
                                                                                                                                                 tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaa                                                                                                                                                                                                                                >  6:103965/1‑149 (MQ=255)
                                                                                                                                                  gAAAAATATCCCACCCCAACGAAGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACAATGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                                                                                                                                                                                                               <  5:275425/149‑1 (MQ=255)
                                                                                                                                                      aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                                                                                                                                                                                                              >  8:213485/1‑130 (MQ=255)
                                                                                                                                                          tCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                                                                                                                                       <  6:266755/149‑1 (MQ=255)
                                                                                                                                                           cccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCACTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                                                                                                                                                                                                                      >  8:243942/1‑149 (MQ=255)
                                                                                                                                                            ccACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                                                                                                                                     <  5:52076/149‑1 (MQ=255)
                                                                                                                                                                 ccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCACAACAGTAATGACAGAAACAGAATTATTCTCTTCAACCGCTTACGTTGAGTGGGCTGCATATATAGTTCCTGAATCTGTAAGTCTTTTAGTTATAGGTTACTGTAAATGTTGCTGGTGCAGAAAACTTACCCGCAGTGAACTTGTTGACGGtttttttt                                                                                  >  4:85643/1‑275 (MQ=255)
                                                                                                                                                                           cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTa                                                                                                                                                                                                      >  5:330774/1‑149 (MQ=255)
                                                                                                                                                                               aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt                                                                                                                                                                                                  <  5:455460/149‑1 (MQ=255)
                                                                                                                                                                                tCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTg                                                                                                                                                                                                 >  8:148413/1‑149 (MQ=255)
                                                                                                                                                                                    tatTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa                                                                                                                                                                                                >  4:529703/1‑146 (MQ=255)
                                                                                                                                                                                     tttccccTGGTTACTTATTCCGCTTTTGCAGGGTCTTTTATTTACATGGGCAACAAGGTCAAAATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGGAATGACAGAACCAGAATTATTCTCTTCATCAGCTTACTGTGAGTGGGATGCATATATAGTTACTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTg                                                               <  4:78314/270‑1 (MQ=255)
                                                                                                                                                                                     atTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                                                                                                                            >  6:365767/1‑149 (MQ=255)
                                                                                                                                                                                      ttACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                                                                                                                                                                                           >  5:162021/1‑149 (MQ=255)
                                                                                                                                                                                      ttACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTggg                                                                                                                                                                                           >  7:340389/1‑149 (MQ=255)
                                                                                                                                                                                       tACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          <  8:308217/149‑1 (MQ=255)
                                                                                                                                                                                        aCCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                                                                                                                            <  2:529708/146‑1 (MQ=255)
                                                                                                                                                                                              gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                                                                                                                                      >  4:100541/1‑146 (MQ=255)
                                                                                                                                                                                                  cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGTGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatata                                                                                                                                                                                  <  2:100542/146‑1 (MQ=255)
                                                                                                                                                                                                  cTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAg                                                                                                                                                                               <  5:449173/149‑1 (MQ=255)
                                                                                                                                                                                                    gATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATAGATAGTTCCTGAATCTGTAAGTCAGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACCCATCGGataat                                                >  2:466700/1‑274 (MQ=255)
                                                                                                                                                                                                       tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCt                                                                                                                                                                          <  6:330774/149‑1 (MQ=255)
                                                                                                                                                                                                            ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATc                                                                                                                                                                     >  5:97472/1‑149 (MQ=255)
                                                                                                                                                                                                               tGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGt                                                                                                                                                                  <  7:454290/149‑1 (MQ=255)
                                                                                                                                                                                                                tcAGGGGGTGGTTTTTACATTGGAAAAAAGGGCAATATACTGTCAAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCAtt                                   <  4:70205/274‑1 (MQ=255)
                                                                                                                                                                                                                 cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                                                                                                                                                <  5:103965/149‑1 (MQ=255)
                                                                                                                                                                                                                 cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                                                                                                                                                <  7:148413/149‑1 (MQ=255)
                                                                                                                                                                                                                     tGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCt                                                                                                                                                            <  6:276058/149‑1 (MQ=255)
                                                                                                                                                                                                                       cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgt                                                                                                                                                          >  8:122584/1‑149 (MQ=255)
                                                                                                                                                                                                                         gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgtta                                                                                                                                                        <  7:300631/149‑1 (MQ=255)
                                                                                                                                                                                                                          ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTTAATCTGTAAGTCTgttag                                                                                                                                                       >  8:286528/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttgttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTgg                        <  2:446907/272‑1 (MQ=255)
                                                                                                                                                                                                                               tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                  <  7:56571/149‑1 (MQ=255)
                                                                                                                                                                                                                               tACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                  <  6:97472/149‑1 (MQ=255)
                                                                                                                                                                                                                                aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                                                                                                                                                 <  8:340389/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              >  7:234849/1‑103 (MQ=255)
                                                                                                                                                                                                                                 cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGt                                                                                                                                                                                              <  8:234849/103‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtata                                                                                                                                                  <  6:49453/143‑1 (MQ=255)
                                                                                                                                                                                                                                      gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACATAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACt                                                                                                                                           >  8:33343/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTg                                                                                                                                                           <  4:523848/131‑1 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCg                                                                                                            <  2:85779/178‑1 (MQ=255)
                                                                                                                                                                                                                                           aaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAATTATAGGTTACTGTaaa                                                                                                                                      >  5:52355/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAt                                                                                                                                     >  8:201257/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt                                                                                                                                   >  6:259473/1‑149 (MQ=255)
                                                                                                                                                                                                                                               taaaTATCCTGTCCAGCAATGGCTGGAAAGGCAATGAAAAAAAGTAATGCCAGAAACAGAATTATTCCCTTAATCAGCTTACTTGGAGTGGGAGGCATATATATTTCCTGAATTTGTAAGTCTGTAGTTATAGGTTACTGTAAATGTgg                                                                                                                                 <  7:286528/147‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGc  <  2:72182/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  atcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          <  7:149056/86‑1 (MQ=38)
                                                                                                                                                                                                                                                  atcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGa                                                                                                                                                                                          >  8:149056/5‑90 (MQ=38)
                                                                                                                                                                                                                                                      ccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGtt                                                                                                                                                                             >  7:82928/1‑99 (MQ=255)
                                                                                                                                                                                                                                                       cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                                                                                                                                        >  6:78178/1‑103 (MQ=255)
                                                                                                                                                                                                                                                       cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGa                                                                                                                                                                        <  5:78178/103‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGACTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                                                                                                                                 <  6:388674/108‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                                                                                                                                                 >  5:388674/1‑108 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTTCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGAc                                         >  4:326005/1‑228 (MQ=255)
                                                                                                                                                                                                                                                            cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                                                                                                                                            <  3:61639/78‑1 (MQ=255)
                                                                                                                                                                                                                                                            cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGga                       <  2:352657/243‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaaa                                                                                                                                      >  6:308389/1‑131 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaaa                                                                                                                                      <  5:308389/131‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGc                                                              <  1:40678/203‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCa                                     <  2:326009/228‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGGTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                                                                                   >  8:362692/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACtt                                                                                                                 >  7:86545/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa                                                                                                                <  5:336844/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                 
AGCCAGCAGCCATGACACCGCTAAGGGCTCTAAAATCAATATGTAACTGATTTAGCACATCGCCAGTTCCCTGTAATGCACAGTCGGATACAGTGATAGTGAACTCTTTTCCACCTGGTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCATCGGATAATAAACGACCGCATTGGCTGCAGTGGAGATTGTTACCTCTTGTGATGC  >  NC_000913/151263‑151778

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