Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,693 G→T 5.2% T395K (ACG→AAG)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,318,6930GT 5.2% 584.1 / 14.1 166T395K (ACG→AAG) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (76/81);  new base T (5/4);  total (81/85)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 170 total aligned reads displayed.

ATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACC  >  NC_000913/2318419‑2318956
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                         
 tACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCg                                                                                                                                                                                                                                                                           <  2:251930/274‑1 (MQ=255)
    cgccgcgTGCGCAAGTCATGTATGACGGTGGCAAGGACCCTCGATTGTGACTGCCTCGCCTGTTCCGCTGCGTGTGCCCTCCCCTGCAACGACTCATCCATCTTCACGCGCGAAGAAACCACCACCAGCACACAATTGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAg                                                                                                                                                                                                                                                                       <  4:123191/271‑1 (MQ=255)
     gcggcgTTAGCAGCTCATTCCTGACGGTGGCAAGCAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCATCACGCGCGAAGCAACATCCACCCGTACACAAATTGCCACGTGTTAATTACGATAGCGCGAATAGACTAAGCTCATTCGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGGCCACAGATAATTTGCGTCGGTCGTTGGCGGTCCTCATGCGTAAGc                                                                                                                                                                                                                                                                      <  4:384750/275‑1 (MQ=255)
      cggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGTTTGTGCCATCTCCTGCAACGACTCTTCCTTCTTCACGCACCAAGAAACCTCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGTCGATGCTGATTTGCAGATCGCTATTGTTGCTGGTCAGGACATTAACAAAATTGAACTGCTGCCCACAGACAATTGGCGGAAGTAGTTGGAGGTCCCCATGCTTAAGCa                                                                                                                                                                                                                                                                     >  4:191361/1‑275 (MQ=255)
       ggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCat                                                                                                                                                                                                                                                                    <  2:263769/275‑1 (MQ=255)
        gcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCCCACGCGTAAGCat                                                                                                                                                                                                                                                                    >  1:126538/1‑274 (MQ=255)
         cgTTGGCAGCTCATGCCTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCatat                                                                                                                                                                                                                                                                  <  4:346249/275‑1 (MQ=255)
         ccttCGCAGCTCATGACTGACGGTGTCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCatat                                                                                                                                                                                                                                                                  <  2:126398/273‑1 (MQ=255)
             cGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                                                                                                                                              >  2:524121/1‑275 (MQ=255)
               cagcaaagaactgacGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑t                                                                                                                                                                                                                                                           <  2:125376/266‑1 (MQ=255)
                            acGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGACTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATTCGTCAGCCTATTGAGA‑TCGGTATgcgcacg                                                                                                                                                                                                                                              >  1:5765/1‑272 (MQ=255)
                                tGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGATGCTTGTGCCATCGCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCg                                                                                                                                                                                                                                          <  2:278962/275‑1 (MQ=255)
                                          atgcggttttTACCGGCTCGCCTGTGCCGCTGCTGTGGCCCTCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACCAATGGCCACGTTTTCATTACGATAGAGCGAATGGACGAAGCTGATTTGAAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                                                                                                                                                 <  4:332456/268‑1 (MQ=255)
                                                                                        tGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCTACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACCCAAAATGCAAATACCCACTgc                                                                                                                                                                                     >  4:15350/1‑271 (MQ=255)
                                                                                        tGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgaga                                                                                                                                                                                                                                                             <  2:109498/201‑1 (MQ=255)
                                                                                         gCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGGTTTGCAAATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGcgcg                                                                                                                                                                                 <  2:372430/275‑1 (MQ=255)
                                                                                               aCTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTTACGCCATCAGCGGAACGCAAAATGCAAATACCCAGTGGCCCGGAggac                                                                                                                                                                           >  4:529995/1‑273 (MQ=255)
                                                                                                 tctTCCATCTTCACGCCCGAAGAAACCTCCACCCGCCCCCAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTTATTTGCAGATTGGTATTGTGGCTGGTAAGGACATCAACAAAATTGACCTGCTGCCCACAAATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGAc                                                                                                                                                                         <  4:139025/275‑1 (MQ=255)
                                                                                                        tCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCtt                                                                                                                                                                   >  4:475924/1‑274 (MQ=255)
                                                                                                                           aTCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAATCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGtt                                                                                                                                               <  2:33309/275‑1 (MQ=255)
                                                                                                                           aTCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATTATTTGCATCATACGTTGGCGGTCCTCATGCGTAAGCATATTGATA‑TAGGTATTCCCCAGTTCGTTACTTAAACTATTCCCGCCATCAGCGGTACGCACATTGCCGAAACCCACCGGCCAGGGAGCGGCAATTTTGCGATTAACCTGCCCAttct                                                                                                                                               >  4:428202/1‑272 (MQ=255)
                                                                                                                             ccaccaGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCAAATTGAGA‑TAGGTATGCCCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCAACTGGCCGGGACGCAACAATCTGGCGATTGAAATGCCCATTttc                                                                                                                                              >  4:529627/1‑274 (MQ=255)
                                                                                                                               accaGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCg                                                                                                                                                                                                                                                                           >  6:554112/1‑148 (MQ=255)
                                                                                                                                  aGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGt                                                                                                                                                                                                                                                                          >  7:6698/1‑146 (MQ=255)
                                                                                                                                    cacacaAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGTCGCGGATGCGACAATCTTGCGATTGAACTGCTTATGTTCTTCCAGTc                                                                                                                                      >  2:514866/1‑275 (MQ=255)
                                                                                                                                               ccACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑Ta                                                                                                                                                                                                                                                          >  6:674946/1‑148 (MQ=255)
                                                                                                                                                 aCGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGc                                                                                                                                                                                                                                                                      >  7:7613/1‑135 (MQ=255)
                                                                                                                                                 aCGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                                                                                                                                              >  5:95145/1‑143 (MQ=255)
                                                                                                                                                    ttttCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAg                                                                                                                                                                                                                                                                       >  8:47403/1‑131 (MQ=255)
                                                                                                                                                          ttACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCat                                                                                                                                                                                                                                                                    <  6:517479/128‑1 (MQ=255)
                                                                                                                                                           tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTa                                                                                                                                                                                                                                                                         <  5:388291/122‑1 (MQ=255)
                                                                                                                                                           tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTa                                                                                                                                                                                                                                                                         >  6:388291/1‑122 (MQ=255)
                                                                                                                                                           tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑Ta                                                                                                                                                                                                                                                          <  8:708485/136‑1 (MQ=255)
                                                                                                                                                           tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAg                                                                                                                                                                                                                                              >  8:495448/1‑148 (MQ=255)
                                                                                                                                                                     cgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTa                                                                                                                                                                                                                                                      >  5:265001/1‑130 (MQ=255)
                                                                                                                                                                     cgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTa                                                                                                                                                                                                                                                      <  6:265001/130‑1 (MQ=255)
                                                                                                                                                                      gAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACg                                                                                                    <  2:475925/275‑1 (MQ=255)
                                                                                                                                                                       aaTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑t                                                                                                                                                                                                                                                           >  5:546466/1‑123 (MQ=255)
                                                                                                                                                                        aTGGACGAAGCTGATGTTCGGATTGGTATGGTTGCTGGTCAGTACATCAACAAAATGGCCCTGCTGCCCCCAGATAATTTGCGGCAGTCGTTTGCGGGCATCATGCGTAAGCATATTGAGA‑TAGTGATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTc                                                                                                  <  4:207857/275‑1 (MQ=255)
                                                                                                                                                                             cGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatg                                                                                                                                                                                                                                                                   >  7:1036694/1‑109 (MQ=255)
                                                                                                                                                                               aaGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGa                                                                                                                                                                                                                          <  7:495448/148‑1 (MQ=255)
                                                                                                                                                                                    gATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCa                                                                                                                                                                                                                                                                     <  8:937962/101‑1 (MQ=255)
                                                                                                                                                                                       ttGCAGATTGGGATTGTTGGTGGTCACGCCCTCAACAAAGTTTACCTGCTGCCCAATAATCACTTTCGTCAGGAGTTGGCAGTCCTCATGCGTAAGCATATT‑AGAGTAGGGATGCTCCAGTTCGTTACTTAAGAGATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCACGGAGGCGACAAGCTTGCGATTGAACTGATCATGTTCTTCCAGTCGCAGGGCGTCGCATTCCGCAGGAATGAAAATACGTCGCGCGT‑ACATCCGTg                                                                                  <  4:125252/275‑1 (MQ=255)
                                                                                                                                                                                          cAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCTCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAAATGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCga                                                                               >  4:219214/1‑275 (MQ=255)
                                                                                                                                                                                            gATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                                                                                                                                                 <  8:535211/128‑1 (MQ=255)
                                                                                                                                                                                            gATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                                                                                                                                                 >  7:535211/1‑128 (MQ=255)
                                                                                                                                                                                              ttGGGATTGTTGCTGGTCAGGACATCAACAAAAGTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGt                                                                           <  2:529996/275‑1 (MQ=255)
                                                                                                                                                                                               tGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGT‑TAGGTATGCGCCAGTTCGtt                                                                                                                                                                                                                                        >  6:83740/1‑118 (MQ=255)
                                                                                                                                                                                                  tATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                                                                                                                                                 <  7:1065556/122‑1 (MQ=255)
                                                                                                                                                                                                  tATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                                                                                                                                                 >  8:1065556/1‑122 (MQ=255)
                                                                                                                                                                                                    ttgttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTa                                                                                                                                                                                                                                   >  8:826463/1‑118 (MQ=255)
                                                                                                                                                                                                    ttgttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCaaaa                                                                                                                                                                                                     <  5:674946/148‑1 (MQ=255)
                                                                                                                                                                                                       ttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatga                                                                                                                                                                                                                                                                  <  5:945490/85‑3 (MQ=255)
                                                                                                                                                                                                       ttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatga                                                                                                                                                                                                                                                                  >  6:945490/1‑83 (MQ=255)
                                                                                                                                                                                                        tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatga                                                                                                                                                                                                                                                                  <  7:427202/84‑3 (MQ=255)
                                                                                                                                                                                                        tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGcc                                                                                                                                                                                                 >  6:986806/1‑147 (MQ=255)
                                                                                                                                                                                                        tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCa                                                                                                                                                                                                 <  5:986806/148‑1 (MQ=255)
                                                                                                                                                                                                         gCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTTCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTAATGTTCTTACAGTCGCAGGGCGTCGCTTTCCCCCGGAATGAAAATCCGTCGCCCGT‑ACATCCGTTCGGGAGAAAACAATgaaa                                                                >  4:161740/1‑272 (MQ=255)
                                                                                                                                                                                                         gCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAg                                                                >  2:201677/1‑275 (MQ=255)
                                                                                                                                                                                                          cTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                                                                                                                                                            <  7:241396/135‑1 (MQ=255)
                                                                                                                                                                                                          cTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                                                                                                                                                            >  8:241396/1‑135 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATgcgca                                                                                                                                                                                                                                                <  8:1023351/97‑2 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                                                                                                                                                  <  5:227397/143‑1 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                                                                                                                                                  >  6:227397/1‑143 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                                                                                                                                             >  6:131473/1‑148 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                                                                                                                                             >  8:774222/1‑148 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATAGAACTGCTCATGTTCTTCCAGTCGCAGGGCGGCGCTTTCCGCCGGAATGAAAATACGTCCCACGG‑ACATCCGAGCGGGAGTAAACCATGAAGccc                                                             >  4:116492/1‑273 (MQ=255)
                                                                                                                                                                                                            ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGAAGCCATCAGCGGTACGCAAAATGCAGATCCCCACTGGCGCGGAGGCGAACATCTTTCGTTTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCCTTCCGCCGGAATGGAAATAAGGCCCGCGTAACACCGGCGCGAGAG‑AAAAAATGCAGccc                                                             >  4:458823/1‑273 (MQ=255)
                                                                                                                                                                                                              tCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                                                                                                                                                  >  7:537493/1‑141 (MQ=255)
                                                                                                                                                                                                               cAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACt                                                                                                                                                                                                                                     <  5:799324/105‑1 (MQ=255)
                                                                                                                                                                                                               cAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                                                                                                                                                  <  5:204291/140‑1 (MQ=255)
                                                                                                                                                                                                               cAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                                                                                                                                                  >  6:204291/1‑140 (MQ=255)
                                                                                                                                                                                                                 ggACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑Ta                                                                                                                                                                                                                                                          >  8:544460/1‑82 (MQ=255)
                                                                                                                                                                                                                     aTCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTAATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGtt                                                                                                                                                                                                                                        >  5:234444/1‑96 (MQ=255)
                                                                                                                                                                                                                       caacaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                                                                                                                                             <  5:488982/137‑1 (MQ=255)
                                                                                                                                                                                                                        tacaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGcaca                                                  <  2:15350/273‑1 (MQ=255)
                                                                                                                                                                                                                         acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTg                                                                                                                                                                                                                                                                <  8:792802/69‑1 (MQ=255)
                                                                                                                                                                                                                         acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTg                                                                                                                                                                                                                                                                >  7:792802/1‑69 (MQ=255)
                                                                                                                                                                                                                         acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCg                                                                                                                                                                                                                                          >  6:714198/1‑90 (MQ=255)
                                                                                                                                                                                                                          caaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATg                                                                                                                                                                                                                                                    <  6:972613/79‑1 (MQ=255)
                                                                                                                                                                                                                           aaaaTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAg                                                                                                                                                                              <  7:1023107/148‑1 (MQ=255)
                                                                                                                                                                                                                             aaTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGc                                            >  2:397597/1‑275 (MQ=255)
                                                                                                                                                                                                                              aTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGcc                                                                                                                                                                                                                                                >  7:819796/1‑79 (MQ=255)
                                                                                                                                                                                                                              aTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                                                                                                                                                  >  7:9747/1‑125 (MQ=255)
                                                                                                                                                                                                                                  aCCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAAt                                       <  2:219215/275‑1 (MQ=255)
                                                                                                                                                                                                                                  aCCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAAt                                       >  2:422412/1‑275 (MQ=255)
                                                                                                                                                                                                                                            cacaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGt                                                                                                                                                                                                                                         <  7:257071/72‑1 (MQ=255)
                                                                                                                                                                                                                                            cacaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGt                                                                                                                                                                                                                                         >  8:257071/1‑72 (MQ=255)
                                                                                                                                                                                                                                             acaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTg                                                                                                                                                                  >  6:886470/1‑142 (MQ=255)
                                                                                                                                                                                                                                             acaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTg                                                                                                                                                                  <  5:886470/142‑1 (MQ=255)
                                                                                                                                                                                                                                              caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTaa                                                                                                                                                                                                                                                                        >  8:1003310/1‑40 (MQ=255)
                                                                                                                                                                                                                                              caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCg                                                                                                                                                                           >  8:350196/1‑132 (MQ=255)
                                                                                                                                                                                                                                              caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCg                                                                                                                                                                           <  7:350196/132‑1 (MQ=255)
                                                                                                                                                                                                                                              caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCCGCCCCGGCAAGCCCATCCAGCAAAATTGCGTTAAGGAt                           >  4:284213/1‑275 (MQ=255)
                                                                                                                                                                                                                                               aGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                                                                                                                                                            <  7:93392/98‑1 (MQ=255)
                                                                                                                                                                                                                                                        gCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGc                                                                                                                                                                 <  6:262035/132‑1 (MQ=255)
                                                                                                                                                                                                                                                            cAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCa                                                                                                                                                                                         >  5:708484/1‑104 (MQ=255)
                                                                                                                                                                                                                                                              gTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAg                                                                                                                                                                              <  5:1026817/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tgggtcccctggagTAGGCATATTGGGA‑TAGGTATGCGGCAGTTCGTTACTTAAAATATTGAGGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCa    <  4:343934/262‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGc                                                                                                                               >  8:602454/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                         
ATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGA‑TAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGT‑ACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACC  >  NC_000913/2318419‑2318956

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