Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 963,286 A→G H79R (CAC→CGC)  himD → integration host factor beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913963,2860AG98.3% 542.4 / ‑1.5 172H79R (CAC→CGC) himDintegration host factor beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (92/77);  total (93/79)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 174 total aligned reads displayed.

GCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCG  >  NC_000913/963196‑963375
                                                                                          |                                                                                         
 cAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCtt                                                                                        >  1:629913/1‑93 (MQ=255)
     ttCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTaaa                                                                                    <  1:2288607/93‑1 (MQ=255)
     ttCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTaaa                                                                                    >  1:1731778/1‑93 (MQ=255)
      tctcTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAAc                                                                                   <  1:4039561/93‑1 (MQ=255)
       ctctTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAAcc                                                                                  >  1:2618025/1‑93 (MQ=255)
        tctTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCt                                                                                 >  1:4204029/1‑93 (MQ=255)
        tctTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCt                                                                                 >  1:4582627/1‑93 (MQ=255)
              cACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTaaa                                                                           >  2:1945155/1‑93 (MQ=255)
               aCTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAg                                                                          >  2:2315112/1‑93 (MQ=255)
               aCTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAg                                                                          >  1:844614/1‑93 (MQ=255)
                 tACCGCGCACCACGTACCGGACGTAATCCGATGACTGGCGATAAAGTAGAACTGGAAGTAAAATACGTTCCTCTCTTTATACTTGGTaatgta                                                                        >  2:4354957/1‑89 (MQ=255)
                 tACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGaa                                                                        <  1:2782568/93‑1 (MQ=255)
                 tACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGaa                                                                        >  1:2650889/1‑93 (MQ=255)
                   ccGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACt                                                                      >  1:2030348/1‑93 (MQ=255)
                    cgcgCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTg                                                                     <  2:858476/93‑1 (MQ=255)
                        caccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGcgcg                                                                 >  1:1513493/1‑93 (MQ=255)
                         cccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                <  2:4499070/92‑1 (MQ=255)
                         accacGTACGGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                >  2:3145426/1‑93 (MQ=255)
                         accacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                <  1:3954280/93‑1 (MQ=255)
                         accacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                >  2:1457354/1‑93 (MQ=255)
                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               <  2:1294033/93‑1 (MQ=255)
                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               <  2:2752575/93‑1 (MQ=255)
                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               <  1:755919/93‑1 (MQ=255)
                            acGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATcg                                                             >  2:2486966/1‑93 (MQ=255)
                            acGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATcg                                                             >  2:4292609/1‑93 (MQ=255)
                              gTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATcgcg                                                           >  1:2401705/1‑93 (MQ=255)
                               tACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCgcgc                                                          >  2:1438494/1‑93 (MQ=255)
                                    gACGTAATGCGAAGACAGACGATAAAGGAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAata                                                     <  1:3173983/93‑1 (MQ=255)
                                    gACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAata                                                     >  2:3084551/1‑93 (MQ=255)
                                      cGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATAtt                                                   <  2:723750/93‑1 (MQ=255)
                                       gTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATAttt                                                  >  1:2588151/1‑93 (MQ=255)
                                       gTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATAttt                                                  >  1:4723391/1‑93 (MQ=255)
                                        tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTa                                                 <  2:3074730/93‑1 (MQ=255)
                                        tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTa                                                 >  2:550320/1‑93 (MQ=255)
                                        tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGAAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTa                                                 <  1:197174/93‑1 (MQ=255)
                                         aaTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTAc                                                >  1:3432130/1‑93 (MQ=255)
                                          aTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACg                                               <  2:786189/93‑1 (MQ=255)
                                           tCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACgg                                              <  1:1132788/93‑1 (MQ=255)
                                             cGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGtt                                            <  1:3574178/93‑1 (MQ=255)
                                             cGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGtt                                            >  1:330961/1‑93 (MQ=255)
                                             cGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGtt                                            >  2:235479/1‑93 (MQ=255)
                                              gAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTa                                           >  1:2854517/1‑93 (MQ=255)
                                              gAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTa                                           <  1:2473702/93‑1 (MQ=255)
                                               aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTaa                                          >  2:1669338/1‑93 (MQ=255)
                                               aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTaa                                          <  2:888161/93‑1 (MQ=255)
                                               aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTaa                                          <  2:3560514/93‑1 (MQ=255)
                                               aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTaa                                          >  1:1775364/1‑93 (MQ=255)
                                                aGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAg                                         <  1:29992/93‑1 (MQ=255)
                                                aGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAg                                         >  1:1266664/1‑93 (MQ=255)
                                                aGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAg                                         <  1:29991/93‑1 (MQ=255)
                                                aGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAg                                         >  1:2741586/1‑93 (MQ=255)
                                                 gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGt                                        <  1:2823005/93‑1 (MQ=255)
                                                 gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGt                                        >  2:3984944/1‑93 (MQ=255)
                                                    tGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtttt                                     >  1:103747/1‑93 (MQ=255)
                                                     ggCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttttt                                    >  1:3940645/1‑93 (MQ=255)
                                                        gATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTAc                                 <  1:2319612/93‑1 (MQ=255)
                                                        gATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTAc                                 >  2:101374/1‑93 (MQ=255)
                                                           aaaGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCa                              >  2:3363442/1‑93 (MQ=255)
                                                            aaGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaa                             <  2:796732/93‑1 (MQ=255)
                                                             aGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaaa                            >  1:148683/1‑93 (MQ=255)
                                                             aGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaaa                            >  2:208251/1‑93 (MQ=255)
                                                              gTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAAc                           >  2:2574335/1‑93 (MQ=255)
                                                              gTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAAc                           >  1:4404244/1‑93 (MQ=255)
                                                               tAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACt                          >  2:3482500/1‑93 (MQ=255)
                                                                aGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACtt                         <  1:4191811/93‑1 (MQ=255)
                                                                aGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACtt                         <  2:2841817/93‑1 (MQ=255)
                                                                aGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACtt                         >  2:4640748/1‑93 (MQ=255)
                                                                 gAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTg                        >  2:1068027/1‑93 (MQ=255)
                                                                 gAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTg                        >  2:3360619/1‑93 (MQ=255)
                                                                   aCTGGAAGGAAACTACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGaa                      <  2:572172/93‑1 (MQ=255)
                                                                   aCTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGaa                      >  2:185124/1‑93 (MQ=255)
                                                                    cTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAAc                     <  1:3143725/93‑1 (MQ=255)
                                                                    cTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAAc                     <  1:2374750/93‑1 (MQ=255)
                                                                     tGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACg                    <  1:4242536/93‑1 (MQ=255)
                                                                     tGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACg                    >  1:4173662/1‑93 (MQ=255)
                                                                      ggaaggaaTATTCGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGGCAATATTTACGGTGAAGTTATAGACTCAAACTATAACga                   >  2:4225313/1‑93 (MQ=255)
                                                                       gaaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgag                  >  2:3222406/1‑93 (MQ=255)
                                                                        aaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTATTCAAACTTGAACgaga                 <  1:2257626/93‑1 (MQ=255)
                                                                        aaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgaga                 <  2:619510/93‑1 (MQ=255)
                                                                          ggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagaga               >  1:3734209/1‑93 (MQ=255)
                                                                          ggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagaga               >  2:4063334/1‑93 (MQ=255)
                                                                           gaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGaa              >  2:3435703/1‑93 (MQ=255)
                                                                                tACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCa         <  2:2519540/93‑1 (MQ=255)
                                                                                tACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCa         <  2:3677084/93‑1 (MQ=255)
                                                                                  cGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAcc       <  2:3821774/93‑1 (MQ=255)
                                                                                  cGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAcc       <  2:1753421/93‑1 (MQ=255)
                                                                                  cGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAcc       <  2:4710130/93‑1 (MQ=255)
                                                                                   gTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt      >  1:2560556/1‑93 (MQ=255)
                                                                                   gTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt      >  2:4701346/1‑93 (MQ=255)
                                                                                   gTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt      >  2:2163365/1‑93 (MQ=255)
                                                                                   gTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt      >  2:1213738/1‑93 (MQ=255)
                                                                                    ttCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGTGAAAAGCACCTg     >  1:652871/1‑93 (MQ=255)
                                                                                    ttCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTg     >  2:4335786/1‑93 (MQ=255)
                                                                                    ttCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTg     <  2:2624590/93‑1 (MQ=255)
                                                                                    ttCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTg     <  1:1745298/93‑1 (MQ=255)
                                                                                     tCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGt    >  1:1437943/1‑93 (MQ=255)
                                                                                     tCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGt    <  2:535698/93‑1 (MQ=255)
                                                                                     tCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGt    >  2:517205/1‑93 (MQ=255)
                                                                                      ccTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTc   >  1:291723/1‑93 (MQ=255)
                                                                                       cTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCg  >  1:3385573/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCG  >  NC_000913/963196‑963375

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