Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,209,664 (G)5→6 coding (96/834 nt) mcrA → McrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,209,6591.G99.3% 475.4 / ‑3.8 147coding (91/834 nt)mcrAMcrA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (80/66);  total (81/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 189 total aligned reads displayed.

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCG  >  NC_000913/1209567‑1209746
                                                                                             |                                                                                       
ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                          >  1:2310732/1‑93 (MQ=255)
ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                          >  1:598288/1‑93 (MQ=255)
 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTg                                                                                        <  2:2392386/93‑1 (MQ=255)
 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTg                                                                                        >  2:4072234/1‑93 (MQ=255)
 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTg                                                                                        <  2:1059549/93‑1 (MQ=255)
  atGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTgg                                                                                       >  2:4695788/1‑93 (MQ=255)
  atGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTgg                                                                                       <  1:1992822/93‑1 (MQ=255)
   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggg                                                                                      <  1:4476936/93‑1 (MQ=255)
   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggg                                                                                      <  1:1529154/93‑1 (MQ=255)
   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGAGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggg                                                                                      >  1:1897647/1‑93 (MQ=255)
    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTgggg                                                                                     <  2:6019765/93‑1 (MQ=255)
    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTgggg                                                                                     <  2:1598570/93‑1 (MQ=255)
    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTgggg                                                                                     <  1:1544905/93‑1 (MQ=255)
     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGGTAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggg                                                                                    >  2:2080842/1‑93 (MQ=255)
     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggg                                                                                    <  2:3323732/93‑1 (MQ=255)
     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggg                                                                                    <  2:186431/93‑1 (MQ=255)
     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggg                                                                                    <  2:819080/93‑1 (MQ=255)
     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggg                                                                                    <  2:943955/93‑1 (MQ=255)
      aTGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTgggggg                                                                                   >  2:5270327/1‑92 (MQ=255)
       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggggc                                                                                  >  1:3051652/1‑91 (MQ=255)
       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggggc                                                                                  >  1:2683152/1‑91 (MQ=255)
        gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggggcc                                                                                 >  1:3770223/1‑90 (MQ=255)
         tttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTggggggccg                                                                                <  1:108583/93‑5 (MQ=255)
            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCgggg                                                                              >  2:5055000/1‑92 (MQ=255)
            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGt                                                                              <  2:3738632/93‑1 (MQ=255)
             tGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTg                                                                             >  1:5430469/1‑93 (MQ=255)
              gATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGa                                                                            <  1:1704303/93‑1 (MQ=255)
              gATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGa                                                                            >  2:5634715/1‑93 (MQ=255)
                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACa                                                                          <  2:4469654/93‑1 (MQ=255)
                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACa                                                                          >  1:407044/1‑93 (MQ=255)
                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACa                                                                          >  2:2485405/1‑93 (MQ=255)
                 aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAg                                                                         >  2:4564565/1‑93 (MQ=255)
                 aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAg                                                                         <  1:5323898/93‑1 (MQ=255)
                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGa                                                                        >  2:1423658/1‑93 (MQ=255)
                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGa                                                                        >  1:2917890/1‑93 (MQ=255)
                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaa                                                                       >  1:2580274/1‑93 (MQ=255)
                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaa                                                                       >  2:3539865/1‑93 (MQ=255)
                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaa                                                                       <  1:1645956/93‑1 (MQ=255)
                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaaa                                                                      >  1:1176608/1‑93 (MQ=255)
                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaaa                                                                      >  1:5622118/1‑93 (MQ=255)
                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaaa                                                                      <  1:1338469/93‑1 (MQ=255)
                     atGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAAc                                                                     >  2:228242/1‑93 (MQ=255)
                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACa                                                                    >  1:2458646/1‑93 (MQ=255)
                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACa                                                                    >  1:3462107/1‑93 (MQ=255)
                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaa                                                                   <  2:4965415/93‑1 (MQ=255)
                         aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAg                                                                 >  1:3719738/1‑93 (MQ=255)
                         aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAg                                                                 <  1:6099628/93‑1 (MQ=255)
                           ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGAt                                                               >  1:2497237/1‑93 (MQ=255)
                           ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGAt                                                               >  2:782107/1‑93 (MQ=255)
                            tGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGAtt                                                              >  1:3249414/1‑93 (MQ=255)
                             taaCTGAAATCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTa                                                             <  1:1616215/92‑1 (MQ=255)
                              aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTAc                                                            >  2:6191516/1‑93 (MQ=255)
                              aaCTGAAAGCTGAGTGTGCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTAc                                                            <  2:5987115/93‑1 (MQ=255)
                               aCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACa                                                           <  1:1668845/93‑1 (MQ=255)
                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACaa                                                          >  2:5395156/1‑93 (MQ=255)
                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAat                                                         <  2:3179857/93‑1 (MQ=255)
                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAat                                                         >  2:5560463/1‑93 (MQ=255)
                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAat                                                         <  2:2926725/93‑1 (MQ=255)
                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAata                                                        <  1:1690289/93‑1 (MQ=255)
                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAata                                                        <  1:3911022/93‑1 (MQ=255)
                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAata                                                        <  1:2236513/93‑1 (MQ=255)
                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAatatt                                                      >  2:4405580/1‑92 (MQ=255)
                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATc                                                      <  2:5669910/93‑1 (MQ=255)
                                     aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCg                                                     <  1:5477765/93‑1 (MQ=255)
                                     aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCg                                                     >  1:4666079/1‑93 (MQ=255)
                                      gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGc                                                    <  2:3242076/93‑1 (MQ=255)
                                      gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGc                                                    >  2:3820716/1‑93 (MQ=255)
                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  >  1:3311679/1‑93 (MQ=255)
                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  <  2:4182569/93‑1 (MQ=255)
                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  >  2:6258088/1‑93 (MQ=255)
                                         gAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctct                                                 >  1:4139075/1‑93 (MQ=255)
                                         gAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctct                                                 >  2:3873424/1‑93 (MQ=255)
                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                              >  1:5198671/1‑93 (MQ=255)
                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                              >  2:504787/1‑93 (MQ=255)
                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                              <  1:4220160/93‑1 (MQ=255)
                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                              >  1:6283722/1‑93 (MQ=255)
                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                            <  2:1224256/93‑1 (MQ=255)
                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGGCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                            >  1:4479199/1‑93 (MQ=255)
                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                           <  1:5951331/93‑1 (MQ=255)
                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                           <  1:769926/93‑1 (MQ=255)
                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtata                                         <  2:1125003/93‑1 (MQ=255)
                                                  aTAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTatat                                        <  1:5880738/93‑1 (MQ=255)
                                                   tAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                       >  2:5247490/1‑93 (MQ=255)
                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATCTCATTGaa                                 >  2:1536821/1‑93 (MQ=255)
                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                                 <  1:5568124/93‑1 (MQ=255)
                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             >  2:514237/1‑93 (MQ=255)
                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            <  2:3989387/93‑1 (MQ=255)
                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            <  1:1454799/93‑1 (MQ=255)
                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            <  1:4787763/93‑1 (MQ=255)
                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        <  2:6307215/93‑1 (MQ=255)
                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        <  1:2047076/93‑1 (MQ=255)
                                                                  gtggTTATTGGCTAATCCTTGAGGCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        <  1:3154246/93‑1 (MQ=255)
                                                                  gggtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        >  2:2593902/3‑93 (MQ=255)
                                                                       tATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTc                   >  1:1844084/1‑93 (MQ=255)
                                                                         tGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTg                 >  1:4624299/1‑93 (MQ=255)
                                                                           gTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgt               <  2:5343431/93‑1 (MQ=255)
                                                                           gTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgt               >  2:508266/1‑93 (MQ=255)
                                                                               aaTCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTAt           >  1:3032138/1‑93 (MQ=255)
                                                                                    ttGAGTCGTGGGGGGCCGGGTGACAGAGACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa      >  1:4014357/1‑93 (MQ=255)
                                                                                    ttGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa      >  1:3005703/1‑93 (MQ=255)
                                                                                             |                                                                                       
TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCG  >  NC_000913/1209567‑1209746

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