Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,971 C→G 6.0% intergenic (+20/+11) ybaY → / ← ybaZ outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,9710CG 6.0% 762.3 / 30.6 217intergenic (+20/+11)ybaY/ybaZouter membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand):  ref base C (95/109);  new base G (5/8);  total (100/117)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-03

Only 100 of 222 total aligned reads displayed.

AGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGA  >  NC_000913/475702‑476240
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
agtgttttaGCCGAACTCGCGTTATCTTTTGAGGTCGGCTATTACCATGAATGACAAACTGGTCTTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                                                                                                                                                                                                                                                          <  1:185494/267‑1 (MQ=255)
 gATGTTCAGCCGAACGCGCGTATTCTGGTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGgag                                                                                                                                                                                                                                                                         <  2:133438/275‑1 (MQ=255)
    gTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACACACTGGTATTTACCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGAACTGACATTGGTGCAGGTACAGCAAACCGCCGTGACGGTTCAGGCCAGCGGGTGCGCAACGACTACCGTACCTTCGACGACACCAACTCAGGTGAAGCGGACGTCGTCAGTTCCAGCCCCTACACCATATTAACGCGACTACACCCTCCCCGCCCGGAgggg                                                                                                                                                                                                                                                                      >  4:183381/1‑275 (MQ=255)
          cctgaCGCGCGTGTTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTAGTTATCACCGATACCGTTCAGCCGGTGCTCAACCAGGGGGGAACTAAAGCCGACCTCACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                                                                                                                                                                                                                                                                <  2:619115/271‑1 (MQ=255)
              aCGCGCGTATTCTGTGGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGGACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGtt                                                                                                                                                                                                                                                            <  1:210717/275‑1 (MQ=255)
                                           agcGTGAATCAAAAAATGGTATTGAGCACCGGATCCGTTTAGCGTGTGCTCAACCAGGGCGAAACTAAAGCCGCCCTGCCATTGGTGCCGCTACAGCAAACCGCCGGGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCGACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGcc                                                                                                                                                                                                                               <  4:309909/273‑1 (MQ=255)
                                              gTGAATGACAAACTGGTATGTATTACCTATACCGTGAATCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCAGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                                                                                                                                                                                             <  2:263233/274‑1 (MQ=255)
                                              gTGAATGACAAACTGGGATTTATAACCGATAACGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                                                                                                                                                                                             <  4:484772/274‑1 (MQ=255)
                                                      cAAACTGGTATGTATCACCCATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTGTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGa                                                                                                                                                                                                                     <  2:288742/274‑1 (MQ=255)
                                                              tATTTATCACCGATACCGGTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTAAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCGTCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATca                                                                                                                                                                                                             <  2:141777/274‑1 (MQ=255)
                                                                   aTCACCGATCCGGTCCAGGCGGGGCTCAACCATGGAGCAAGTAAAGCGGACCTGATATTGGGGCCGGTACAGAAAACCCCCGTCCCGGTTCAGGCCAGCGGGGGCGCAAGGACTACCGTACCTTCGACTTCCCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCGCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCt                                                                                                                                                                                                        <  4:355861/274‑1 (MQ=255)
                                                                   aTCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTAACCAACACAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGCGAGTGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTTCACGCTTCCCGATCCCATCACACCt                                                                                                                                                                                                        >  2:453017/1‑274 (MQ=255)
                                                                     cACCGATACCGTTCAGTCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACACTGCAAGTCGATTTGCCCGCTTCACGATACCAGCACACCTTCt                                                                                                                                                                                                     >  1:88452/1‑275 (MQ=255)
                                                                                ttCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTaa                                                                                                                                                                                             >  2:186092/1‑272 (MQ=255)
                                                                                ttCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                                                                                                                                                          >  4:485996/1‑275 (MQ=255)
                                                                                  cAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGccc                                                                                                                                                                                                                                                                             >  2:222193/1‑190 (MQ=255)
                                                                                     ccGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGAAAGGCGATTTg                                                                                                                                                                                                                                 >  2:190842/1‑231 (MQ=255)
                                                                                      cGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGgag                                                                                                                                                                                                                                                                         <  4:222193/190‑1 (MQ=255)
                                                                                         tGAGCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGGGCCGTGACAGCAAACCGCCGTGCCGGTTCAGGCCCGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                                                                                                                                                                                             <  4:190842/231‑1 (MQ=255)
                                                                                          caacccccgggggggCACTAAAGCCCACTNNCGTTTGGGGCAGGTCCAGCACCGCGCCGGCNCGGTTCAGGCCAGCGGTGGCGGCACGCCTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCANTTCCAGCTCCTACGCAGTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                                                                                                                                                <  3:95095/266‑1 (MQ=255)
                                                                                                         aaCTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTg                                                                                                                                                                                                                                             >  4:650910/1‑199 (MQ=255)
                                                                                                             aaaGCCGCCCTGACATTGGTTCCGGTNCAGCAANCCGNCGTGCCGNTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCgg                                                                                                                                                                                                                                                                           <  3:29159/165‑1 (MQ=255)
                                                                                                             aaaGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAg                                                                                                                                                                                                                                         <  2:650910/199‑1 (MQ=255)
                                                                                                                                 gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGagag                                                                                                                                                                                                                                                                       <  5:850707/149‑1 (MQ=255)
                                                                                                                                 gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGagag                                                                                                                                                                                                                                                                       >  8:975295/1‑149 (MQ=255)
                                                                                                                                  ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                                                                                                                                                      <  5:871306/149‑1 (MQ=255)
                                                                                                                                   cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                                                                                                                                                     <  6:1093990/149‑1 (MQ=255)
                                                                                                                                   cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                                                                                                                                                     <  5:742113/149‑1 (MQ=255)
                                                                                                                                   cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                                                                                                                                                     >  7:744384/1‑149 (MQ=255)
                                                                                                                                         aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                                                                                                                                                               <  6:748619/149‑1 (MQ=255)
                                                                                                                                           cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTagtagt                                                                                                                                                                                                                                                             <  5:654464/149‑1 (MQ=255)
                                                                                                                                              aCCGGCGGGCCGGTTCCGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCAACAACTCAGGAGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc                                                                                                                             <  2:268069/274‑1 (MQ=255)
                                                                                                                                                     tGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCa                                                                                                                     <  2:485996/275‑1 (MQ=255)
                                                                                                                                                     ggCCGGTTCAGGCCAGCGGTGGCCCAACGACTACCTTACCTTCGACTTCACCAACTCAGGTGAATCCGACTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGcc                                                                                                                      <  4:438541/273‑1 (MQ=255)
                                                                                                                                                       ccGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAAc                                                                                                                                                                                                                                                 <  7:629638/149‑1 (MQ=255)
                                                                                                                                                            tCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                                                                                                                                                            <  6:875150/149‑1 (MQ=255)
                                                                                                                                                                gCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGCCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCCCACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTACGGTTAGCGACATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGtg                                                                                                           >  2:66516/1‑274 (MQ=255)
                                                                                                                                                                   aGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGa                                                                                                                                                                                                                                     >  5:380109/1‑149 (MQ=255)
                                                                                                                                                                      ggTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                                                                                                                                                  <  7:816761/149‑1 (MQ=255)
                                                                                                                                                                        tGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCTCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGc                                                                                                                                                                                                                                <  8:488295/149‑1 (MQ=255)
                                                                                                                                                                            gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                                                                                                                                                            >  7:908996/1‑149 (MQ=255)
                                                                                                                                                                               aCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGgaga                                                                                            <  1:117159/274‑1 (MQ=255)
                                                                                                                                                                                   cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGa                                                                                                                                                                                                                     <  7:550255/149‑1 (MQ=255)
                                                                                                                                                                                    tACCGTACCTTCGACTTCACAAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAt                                                                                                                                                                                                                    <  5:302404/149‑1 (MQ=255)
                                                                                                                                                                                           ccTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATca                                                                                                                                                                                                             <  7:410177/149‑1 (MQ=255)
                                                                                                                                                                                              tCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGc                                                                                                                                                                                                                                >  1:37758/1‑127 (MQ=255)
                                                                                                                                                                                                   ttCACCAACTCCGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                                                                                                                                                     >  7:613904/1‑149 (MQ=255)
                                                                                                                                                                                                      aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGcc                                                                                                                                                                                                  >  5:787825/1‑149 (MQ=255)
                                                                                                                                                                                                           cTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTaa                                                                                                                                                                                             >  6:474986/1‑149 (MQ=255)
                                                                                                                                                                                                              aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGgaga                                                                                                                                                                                                                                                                        <  8:128638/71‑1 (MQ=39)
                                                                                                                                                                                                              aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGgaga                                                                                                                                                                                                                                                                        >  7:128638/1‑71 (MQ=39)
                                                                                                                                                                                                               ggTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGcc                                                                                                                                                                                         <  8:613904/149‑1 (MQ=255)
                                                                                                                                                                                                                gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                                                                                                                                                     >  6:722879/1‑72 (MQ=255)
                                                                                                                                                                                                                 tGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTg                                                                                                                                                                                       >  6:308283/1‑149 (MQ=255)
                                                                                                                                                                                                                    aTCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCg                                                                                                                                                                                    >  8:253472/1‑149 (MQ=255)
                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCTGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGgcgcg                                                                <  2:577262/264‑1 (MQ=255)
                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc                                                                                                                                                                                   <  7:728923/149‑1 (MQ=255)
                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc                                                                                                                                                                                   >  5:1080373/1‑149 (MQ=255)
                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                                                                                                                                                >  3:169819/1‑184 (MQ=255)
                                                                                                                                                                                                                      ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCt                                                                                                                                                                                  <  6:593313/149‑1 (MQ=255)
                                                                                                                                                                                                                      acGTCTTCGGCAGTTCCAACTCCTAGACAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCt                                                                                                                                                                                  >  6:392549/2‑149 (MQ=255)
                                                                                                                                                                                                                       cGTCTTCGGCAATTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGc                                                                                                                                                                                                                                               >  7:613151/1‑87 (MQ=255)
                                                                                                                                                                                                                         tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAc                                                                                                                                                                               >  8:756316/1‑149 (MQ=255)
                                                                                                                                                                                                                          cTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGcc                                                                                                                                                                                                                               <  3:40638/100‑1 (MQ=255)
                                                                                                                                                                                                                               gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                                                                                                                                                  <  6:862481/92‑1 (MQ=255)
                                                                                                                                                                                                                               gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                                                                                                                                                         >  5:101280/1‑149 (MQ=255)
                                                                                                                                                                                                                               gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCAGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                                                                                                                                                         >  6:803065/1‑149 (MQ=255)
                                                                                                                                                                                                                                 aGTTCCAGCTCCTACGCAATATTTAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATATCGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTaaa                                                                                                                                                                       <  5:721246/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ttCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                                                                                                                                                               <  5:386432/107‑1 (MQ=255)
                                                                                                                                                                                                                                     ccAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccg                                                                                                                                                                   <  5:876679/149‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCATTAATGCCTGTCGCTGACGCTGTAAAtccgg                                                                                                                                                                  >  8:642271/1‑149 (MQ=255)
                                                                                                                                                                                                                                        gCTCCTCCGCAATATTAAGCGCCCGTAACCCTCTCCGCCCGGGGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCca                                  <  4:28699/275‑1 (MQ=255)
                                                                                                                                                                                                                                          tCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                                                                                                                                                 >  5:642629/1‑98 (MQ=255)
                                                                                                                                                                                                                                                gCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGGTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAATTCCGGTCCGGTTAGc                                                                                                                                                        >  5:945881/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 cAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                                                                                                                                                        <  5:1061651/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                                                                                                                                                        >  6:1061651/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    tatTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaaa                                                                                                                                                    <  8:554836/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    tatTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAgg                                                                                                >  4:703333/1‑201 (MQ=255)
                                                                                                                                                                                                                                                      ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                                                                                                                                                  <  6:992933/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCACGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                                                                                                                                                  >  8:159345/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aaGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTGATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                                                                                                                                                <  7:469312/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCg                                                                                                                                                       <  2:242566/142‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                                                                                                                                                >  7:953125/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aaGCGACGTGAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAgtg                    <  2:312654/273‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgtg                                                                                                                                               >  5:1028581/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGc                                                                                                                                              <  7:303644/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGt                                                                                                                                           <  8:541418/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTg                                                                                                                                                                                                                                 >  3:100692/1‑62 (MQ=255)
                                                                                                                                                                                                                                                               cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCt                                                                                                                                                                                                        >  4:34233/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                gTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTgg                                                                                                                                         >  6:366484/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc                                                                                                                                        <  6:18176/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCg                                                                                                                                       >  7:5677/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTACCGAAATTGTGCCGTGTCGa                                                                                                                                      >  8:78890/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                                                                                                                                                                                             <  1:100692/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTg                                                                                                                                                                                                    <  2:34233/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACATTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGTAATTGTGCCGTGGCGAt                                                                                                                                     >  8:305027/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAtt                                                                                                                                    >  5:436053/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAtt                                                                                                                                    >  8:78133/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCCCGGGGCGATCCCGCCCGTTTCGCCACATCACCTGAAGTGGGGACATAGCCCTAgggg  >  4:583957/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                             gCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGgaga                                                                                            <  4:250895/181‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
AGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGA  >  NC_000913/475702‑476240

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