Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,210,134 A→G Y189C (TAT→TGT)  mcrA → McrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,210,1340AG83.5% 475.9 / 90.8 212Y189C (TAT→TGT) mcrAMcrA
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base G (86/91);  total (105/107)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 215 total aligned reads displayed.

ATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTTTTATTTAAAT  >  NC_000913/1210044‑1210225
                                                                                          |                                                                                           
atatGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTtgtg                                                                                           >  2:396212/1‑93 (MQ=255)
atatGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATg                                                                                           >  2:7076825/1‑93 (MQ=255)
 tatGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTgtgt                                                                                          <  1:4934520/93‑1 (MQ=255)
  atGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTc                                                                                         >  1:7168377/1‑93 (MQ=255)
  atGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTc                                                                                         >  1:3646966/1‑93 (MQ=255)
  atGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTc                                                                                         >  2:3726324/1‑93 (MQ=255)
   tGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTcc                                                                                        <  1:4609638/93‑1 (MQ=255)
   tGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTATcc                                                                                        >  1:5925176/1‑93 (MQ=255)
   tGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACACGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTcc                                                                                        >  1:7516936/1‑93 (MQ=255)
     aGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCga                                                                                      <  1:2954949/93‑1 (MQ=255)
     aGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCga                                                                                      <  2:3891142/93‑1 (MQ=255)
     aGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCga                                                                                      <  2:431289/93‑1 (MQ=255)
      gggTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCgag                                                                                     >  2:6151721/1‑93 (MQ=255)
       ggTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCgaga                                                                                    >  1:2691338/1‑93 (MQ=255)
       ggTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCgaga                                                                                    <  1:2292196/93‑1 (MQ=255)
       ggTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCgaga                                                                                    <  2:5095657/93‑1 (MQ=255)
         ttAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGAcc                                                                                  >  2:754474/1‑93 (MQ=255)
           aGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCg                                                                                >  1:1555872/1‑93 (MQ=255)
            gTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGa                                                                               >  1:5894324/1‑93 (MQ=255)
              aaaTTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATg                                                                             <  1:1449094/93‑1 (MQ=255)
              aaaTTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATg                                                                             <  2:1984788/93‑1 (MQ=255)
              aaaTTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATg                                                                             >  1:2748776/1‑93 (MQ=255)
               aattaattAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATgg                                                                            <  2:570401/93‑1 (MQ=255)
                attaattaAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGt                                                                           >  1:7318132/1‑93 (MQ=255)
                 ttaattaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTa                                                                          >  2:6330941/1‑93 (MQ=255)
                 ttaattaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTa                                                                          >  1:4002269/1‑93 (MQ=255)
                  taattaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTaa                                                                         <  2:2193549/93‑1 (MQ=255)
                  taattaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTaa                                                                         >  2:1413014/1‑93 (MQ=255)
                   aattaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTaaa                                                                        >  1:4157576/1‑93 (MQ=255)
                   aattaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGGAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTaaa                                                                        <  1:4420886/93‑1 (MQ=255)
                    attaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTaaaa                                                                       >  1:2304492/1‑93 (MQ=255)
                     ttaaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAg                                                                      >  2:5551988/1‑93 (MQ=255)
                       aaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCt                                                                    <  1:500322/93‑1 (MQ=255)
                       aaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCt                                                                    <  2:1943897/93‑1 (MQ=255)
                       aaGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCt                                                                    <  1:6081839/93‑1 (MQ=255)
                        aGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCtt                                                                   <  2:4615217/93‑1 (MQ=255)
                         gAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTg                                                                  <  2:6663947/93‑1 (MQ=255)
                          aaaaCGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTgg                                                                 <  2:2160251/93‑1 (MQ=255)
                            aaCGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGAt                                                               <  1:643592/93‑1 (MQ=255)
                             aCGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGAtt                                                              <  2:4149883/93‑1 (MQ=255)
                             aCGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGAtt                                                              >  2:2550538/1‑93 (MQ=255)
                              cGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGAttt                                                             >  1:2085821/1‑93 (MQ=255)
                                ttGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTa                                                           <  2:1725447/93‑1 (MQ=255)
                                 tGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTAc                                                          <  2:540682/93‑1 (MQ=255)
                                 tGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTAc                                                          <  2:719004/93‑1 (MQ=255)
                                  gAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACa                                                         <  2:4818415/93‑1 (MQ=255)
                                    gTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGc                                                       <  1:324221/93‑1 (MQ=255)
                                    gTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGc                                                       <  1:3093453/93‑1 (MQ=255)
                                       aaCCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCaaa                                                    <  2:6489310/93‑1 (MQ=255)
                                       aaCCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCaaa                                                    >  1:3550051/1‑93 (MQ=255)
                                        aCCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAg                                                   >  2:3017916/1‑93 (MQ=255)
                                           cGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTaa                                                >  2:6040621/1‑93 (MQ=255)
                                           cGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTaa                                                >  1:6665491/1‑93 (MQ=255)
                                            gAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTaaa                                               <  2:3516541/93‑1 (MQ=255)
                                            gAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTaaa                                               <  2:20075/93‑1 (MQ=255)
                                            gAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTaaa                                               >  2:7424842/1‑93 (MQ=255)
                                            gAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTaaa                                               <  1:1672739/93‑1 (MQ=255)
                                             aGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAg                                              <  1:1887407/93‑1 (MQ=255)
                                             aGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAg                                              >  2:1213263/1‑93 (MQ=255)
                                              gggCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAgg                                             >  2:1613578/1‑93 (MQ=255)
                                              gggCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAgg                                             >  2:2883585/1‑93 (MQ=255)
                                                gCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGta                                           <  1:3278212/93‑1 (MQ=255)
                                                gCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGta                                           <  2:3555686/93‑1 (MQ=255)
                                                  tCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGtata                                         <  1:4069155/93‑1 (MQ=255)
                                                    cAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATAtg                                       <  1:1035386/93‑1 (MQ=255)
                                                     aGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATAtgt                                      <  2:4167903/93‑1 (MQ=255)
                                                     aGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATAtgt                                      >  2:2259401/1‑93 (MQ=255)
                                                     aGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATAtgt                                      >  1:6241186/1‑93 (MQ=255)
                                                       gAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGa                                    <  1:3879147/93‑1 (MQ=255)
                                                        aaaCCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGaa                                   <  2:5143219/93‑1 (MQ=255)
                                                         aaCCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGaaa                                  >  1:5040197/1‑93 (MQ=255)
                                                          aCCAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAGGGTATATGTGaaaa                                 <  1:137516/93‑1 (MQ=255)
                                                            cAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACt                               <  2:4708337/93‑1 (MQ=255)
                                                            cAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACt                               <  1:1242384/93‑1 (MQ=255)
                                                            cAGTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACt                               >  1:7456133/1‑93 (MQ=255)
                                                              gTTGAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACtgt                             <  2:2434437/93‑1 (MQ=255)
                                                                 gAGGTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGt                          <  1:4282400/93‑1 (MQ=255)
                                                                  aGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTa                         <  1:5828446/93‑1 (MQ=255)
                                                                    gTAGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTaaa                       >  1:6941510/1‑93 (MQ=255)
                                                                     tAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTaaaa                      >  1:4238863/1‑93 (MQ=255)
                                                                      aGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTaaaaa                     >  2:3359667/1‑93 (MQ=255)
                                                                      aGAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTaaaaa                     <  2:2735456/93‑1 (MQ=255)
                                                                       gAAAGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAAt                    <  1:5158201/93‑1 (MQ=255)
                                                                       gAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAAt                    >  2:3131106/1‑93 (MQ=255)
                                                                        aaaGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAAta                   >  2:7374335/1‑92 (MQ=255)
                                                                         aaGACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGc                  >  1:450228/1‑93 (MQ=255)
                                                                           gACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTc                >  1:469162/1‑93 (MQ=255)
                                                                           gACTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTc                <  1:6806942/93‑1 (MQ=255)
                                                                            aCTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTcc               <  2:5256756/93‑1 (MQ=255)
                                                                             cTACAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCg              >  2:6735629/1‑93 (MQ=255)
                                                                               aCAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGtt            <  2:1519245/93‑1 (MQ=255)
                                                                                cAAAAAGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGttt           <  1:5696283/93‑1 (MQ=255)
                                                                                cAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGttt           >  2:7059419/1‑93 (MQ=255)
                                                                                  aaaaGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTttta         >  2:5475263/1‑93 (MQ=255)
                                                                                   aaaGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTtttat        >  2:459438/1‑93 (MQ=255)
                                                                                    aaGTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTtttatt       <  1:2638362/93‑1 (MQ=255)
                                                                                      gTTTGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTtttattta     >  2:4447617/1‑93 (MQ=255)
                                                                                       tttGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTTTTATTTaa    <  1:1396408/93‑1 (MQ=255)
                                                                                       tttGTGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTTTTATTTaa    >  2:1301973/1‑93 (MQ=255)
                                                                                         tATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTTTTATTTAAat  >  2:5899880/1‑93 (MQ=255)
                                                                                          |                                                                                           
ATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTATGTCCGAGACCCGATGGTAAAAGCTTGGATTTTACAGCAAAGTAAAGGTATATGTGAAAACTGTGGTAAAAATGCTCCGTTTTATTTAAAT  >  NC_000913/1210044‑1210225

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