Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,210,786 A→G intergenic (+384/+117) mcrA → / ← ycgW McrA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,210,7860AG99.6% 806.7 / NA 236intergenic (+384/+117)mcrA/ycgWMcrA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (108/126);  total (109/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGC  >  NC_000913/1210696‑1210875
                                                                                          |                                                                                         
gcAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGatcgtc                                                                                         >  1:3821622/1‑93 (MQ=255)
 caatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGAtcgtcg                                                                                        >  2:8424247/1‑93 (MQ=255)
  aatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGa                                                                                       >  1:7736337/1‑93 (MQ=255)
  aatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGa                                                                                       >  1:2192136/1‑93 (MQ=255)
    tcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGaaa                                                                                     <  1:1373626/93‑1 (MQ=255)
     caatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGaaaa                                                                                    <  2:6400584/93‑1 (MQ=255)
     caatGCCAAGACCATAACTTATGACTTCGACCGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGaaaa                                                                                    <  2:6645526/93‑1 (MQ=255)
       atGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACa                                                                                  <  1:2699101/93‑1 (MQ=255)
       atGCCAAGACCGTAACTTATGACTTCGAACGTCGGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACa                                                                                  <  1:8486275/93‑1 (MQ=255)
        tGCCAAGACCGTAACTTATGCCTTCTAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACAt                                                                                 <  2:3647742/93‑1 (MQ=255)
          ccAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATga                                                                               >  1:1156461/1‑92 (MQ=255)
          ccAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGt                                                                               <  2:2226696/93‑1 (MQ=255)
           cAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTa                                                                              <  1:8441298/93‑1 (MQ=255)
            aaGACCGTAACTTATTACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTaa                                                                             <  1:7348412/93‑1 (MQ=255)
              gACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTTAtc                                                                           <  2:8626175/93‑1 (MQ=255)
              gACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAAtc                                                                           >  1:1024194/1‑93 (MQ=255)
              aaCCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAAtc                                                                           <  2:6441164/92‑1 (MQ=255)
                acGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAAtctc                                                                         <  1:5714099/92‑1 (MQ=255)
                   tAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTcc                                                                      >  1:1815314/1‑93 (MQ=255)
                    aaCTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCa                                                                     >  2:2045108/1‑93 (MQ=255)
                     aCTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCAt                                                                    >  1:8730385/1‑93 (MQ=255)
                     aCTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCAt                                                                    >  1:4383270/1‑93 (MQ=255)
                     aCTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCAt                                                                    >  2:3556710/1‑93 (MQ=255)
                      cTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCAtg                                                                   <  2:2934909/93‑1 (MQ=255)
                       ttATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCAtgt                                                                  >  2:7734145/1‑93 (MQ=255)
                        tATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATCTAATCTCTCCAtgtg                                                                 <  2:8709559/93‑1 (MQ=255)
                          tGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGtt                                                               >  2:8392784/1‑93 (MQ=255)
                          tGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGtt                                                               <  1:4026885/93‑1 (MQ=255)
                             cTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTaaa                                                            >  2:1363558/1‑93 (MQ=255)
                              ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAat                                                           >  2:4071243/1‑93 (MQ=255)
                              ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAat                                                           <  1:5777766/93‑1 (MQ=255)
                                 gAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATAtt                                                        >  2:1888227/1‑93 (MQ=255)
                                 gAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATAtt                                                        >  1:1197402/1‑93 (MQ=255)
                                 gAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATAtt                                                        >  1:744941/1‑93 (MQ=255)
                                 gAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATAtt                                                        >  1:7295118/1‑93 (MQ=255)
                                  aaCGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTg                                                       <  1:3988943/93‑1 (MQ=255)
                                  aaCGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTg                                                       <  2:6042128/93‑1 (MQ=255)
                                   aCGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGa                                                      >  2:7866628/1‑93 (MQ=255)
                                   aCGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGa                                                      >  1:93618/1‑93 (MQ=255)
                                    cGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGaa                                                     >  1:6905702/1‑93 (MQ=255)
                                    cGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGaa                                                     <  1:811524/93‑1 (MQ=255)
                                    cGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGaa                                                     >  1:3655718/1‑93 (MQ=255)
                                     gTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTTaaa                                                    >  2:309814/1‑93 (MQ=255)
                                     gTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGaaa                                                    <  2:7566213/93‑1 (MQ=255)
                                     gTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGaaa                                                    <  2:2105482/93‑1 (MQ=255)
                                      tCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAAc                                                   >  2:8600400/1‑93 (MQ=255)
                                      tCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAAc                                                   <  2:1542911/93‑1 (MQ=255)
                                         gATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACggg                                                <  2:6248213/93‑1 (MQ=255)
                                          aTGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGc                                               <  2:116766/93‑1 (MQ=255)
                                            ggAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGt                                             >  2:3714743/1‑93 (MQ=255)
                                              aaGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtat                                           >  2:4944790/1‑93 (MQ=255)
                                              aaGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtat                                           >  1:3012882/1‑93 (MQ=255)
                                               aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                          <  2:5242641/93‑1 (MQ=255)
                                               aGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGtata                                          >  1:2039953/1‑93 (MQ=255)
                                                 gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAac                                        <  1:3505827/93‑1 (MQ=255)
                                                 gcgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAac                                        <  2:3865451/93‑1 (MQ=255)
                                                   gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                      >  1:6152289/1‑93 (MQ=255)
                                                   gcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAacac                                      <  1:2361411/93‑1 (MQ=255)
                                                    cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACg                                     <  2:5081188/93‑1 (MQ=255)
                                                    cTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACg                                     >  1:780132/1‑93 (MQ=255)
                                                     tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGc                                    <  1:6059035/93‑1 (MQ=255)
                                                     tAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGc                                    >  2:4232342/1‑93 (MQ=255)
                                                       aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                  <  2:5434016/93‑1 (MQ=255)
                                                       aGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCGTCCAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                  <  2:6315772/93‑1 (MQ=255)
                                                         ctgctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgt                                <  1:967889/93‑1 (MQ=255)
                                                          tgctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgtt                               <  2:3513858/93‑1 (MQ=255)
                                                          tgctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgtt                               >  1:2445379/1‑93 (MQ=255)
                                                          tgctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgtt                               >  2:8174607/1‑93 (MQ=255)
                                                           gctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                              >  1:6016314/1‑93 (MQ=255)
                                                           gctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                              <  1:6412293/93‑1 (MQ=255)
                                                            ctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                             <  1:6664016/93‑1 (MQ=255)
                                                            ctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                             <  2:6590084/93‑1 (MQ=255)
                                                            ctgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                             <  2:1701611/93‑1 (MQ=255)
                                                             tgAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgtt                            >  2:7105716/1‑93 (MQ=255)
                                                              gAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                           <  1:463606/93‑1 (MQ=255)
                                                              gAAATGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                           >  2:4727091/1‑93 (MQ=255)
                                                                aaTGTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTttta                         <  1:3750964/93‑1 (MQ=255)
                                                                   gTTCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                      <  2:8798543/93‑1 (MQ=255)
                                                                     tCAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                    <  1:4265890/93‑1 (MQ=255)
                                                                      cAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtg                   <  2:334851/93‑1 (MQ=255)
                                                                      cAGAGTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtg                   <  2:2642347/93‑1 (MQ=255)
                                                                       agagTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgt                  <  2:5917089/93‑1 (MQ=255)
                                                                       agagTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgt                  <  1:5156536/93‑1 (MQ=255)
                                                                        gagTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                 >  1:1274692/1‑93 (MQ=255)
                                                                         agTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTgg                <  2:2388128/93‑1 (MQ=255)
                                                                          gTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGa               <  1:3775200/93‑1 (MQ=255)
                                                                          gTTTGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGa               >  1:7587956/1‑93 (MQ=255)
                                                                           tttGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGat              >  2:3880528/1‑93 (MQ=255)
                                                                            ttGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGata             <  2:6465839/93‑1 (MQ=255)
                                                                            ttGGTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGata             <  1:5299872/93‑1 (MQ=255)
                                                                              tgTGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatt           <  1:4500910/92‑1 (MQ=255)
                                                                                tGAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtattat         <  1:5390579/93‑1 (MQ=255)
                                                                                 gAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt        >  2:7260348/1‑93 (MQ=255)
                                                                                 gAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt        <  1:8254403/93‑1 (MQ=255)
                                                                                 gAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt        >  1:8723256/1‑93 (MQ=255)
                                                                                 gAAGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt        >  1:3186496/1‑93 (MQ=255)
                                                                                   aGCGATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa      >  1:7184929/1‑93 (MQ=255)
                                                                                     cGATCGTCGANAACATGTAATNTNTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa    <  1:8563478/93‑1 (MQ=255)
                                                                                      gATCGTCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg   <  2:4702416/93‑1 (MQ=255)
                                                                                       atcgtcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTGTTTATGTGGATGTTATTAATAga  >  1:4851748/1‑92 (MQ=255)
                                                                                          |                                                                                         
GCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGC  >  NC_000913/1210696‑1210875

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