Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,246,188 A→G E213E (GAA→GAG lamB → phage lambda receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,246,1880AG99.5% 680.1 / NA 215E213E (GAA→GAGlamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (114/100);  total (114/101)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 223 total aligned reads displayed.

CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246098‑4246277
                                                                                          |                                                                                         
 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                        >  2:3325811/1‑93 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                       >  1:2111594/1‑93 (MQ=255)
   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      >  1:368787/1‑93 (MQ=255)
     tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                    >  2:3181200/1‑93 (MQ=255)
     tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                    >  2:2221052/1‑93 (MQ=255)
     tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                    >  1:2233819/1‑93 (MQ=255)
      ttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtgt                                                                                   >  1:1775388/1‑93 (MQ=255)
      ttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtgt                                                                                   <  2:87679/93‑1 (MQ=255)
         tGACTATGCAAACGAAACCGCGAACGCCGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCGCATTAGAGCTGGGTGTCGa                                                                                <  2:3576492/93‑1 (MQ=255)
         tGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGa                                                                                <  2:1391495/93‑1 (MQ=255)
         tGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGa                                                                                <  1:2432311/93‑1 (MQ=255)
         tGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGa                                                                                >  1:1424731/1‑93 (MQ=255)
          gACTATACCAACGAAACCGCGAACGACGTTTTCGATTTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGAc                                                                               >  2:2584820/1‑93 (MQ=255)
            ccataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTa                                                                             <  1:3737104/91‑1 (MQ=255)
            cTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACtc                                                                             >  1:488527/1‑92 (MQ=255)
            cTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTa                                                                             >  1:1851265/1‑93 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            <  2:930854/93‑1 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            <  1:1191698/93‑1 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            >  2:1095793/1‑93 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            >  1:4409204/1‑93 (MQ=255)
               taCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACgg                                                                          >  1:1633917/1‑93 (MQ=255)
                aCCAACGAAACCTCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                         <  1:1089872/93‑1 (MQ=255)
                 ccAACGAAGCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        <  2:3445407/93‑1 (MQ=255)
                  cAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCg                                                                       >  1:4270532/1‑93 (MQ=255)
                   aaCGAAACCGCGAACGACGTTTTCGTTGTGCGGTTAGCGCAGATAGAAATCAACCCGGGCGGCAAATTAGAGCTGGGTGTCGACTACgggcag                                                                      >  2:4310559/1‑89 (MQ=255)
                   aaCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                      >  1:999874/1‑93 (MQ=255)
                    aCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTg                                                                     <  2:972840/93‑1 (MQ=255)
                     cGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGc                                                                    >  2:1882498/1‑93 (MQ=255)
                        aaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCaa                                                                 >  1:1983056/1‑93 (MQ=255)
                         aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAAc                                                                <  2:2625273/93‑1 (MQ=255)
                         aCCGCGAACAAAGTTTTCAATGTGCGTTTAGCGCAGATGGCAATCAACCCGGACGGCACATTAGAGCTGGGTGGCGACTACGGGCGTGCCaat                                                                >  2:1519846/1‑92 (MQ=255)
                          ccGGGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               <  1:592955/93‑1 (MQ=255)
                           cgcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACtt                                                              <  1:3430849/93‑1 (MQ=255)
                           cgcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACtt                                                              <  1:4601063/93‑1 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             >  2:1146626/1‑93 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             >  2:2492264/1‑93 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             <  2:2741525/93‑1 (MQ=255)
                             cgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGc                                                            <  2:2047280/93‑1 (MQ=255)
                              gAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCg                                                           <  2:1704843/93‑1 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  2:3143015/1‑93 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          <  2:625668/93‑1 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  2:552980/1‑93 (MQ=255)
                                acgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTg                                                         >  2:902461/1‑93 (MQ=255)
                                   acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                      >  1:4179502/1‑93 (MQ=255)
                                   acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                      >  1:3549112/1‑93 (MQ=255)
                                   acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGGGTCGACTACGGTCGTGCCACCTTGCGTGCTa                                                      >  1:1621711/1‑93 (MQ=255)
                                     gTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAAc                                                    <  2:273918/93‑1 (MQ=255)
                                     gTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAAc                                                    >  1:953808/1‑93 (MQ=255)
                                     gTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAAc                                                    >  1:1399047/1‑93 (MQ=255)
                                       tttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTa                                                  <  1:3645968/93‑1 (MQ=255)
                                        ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTAt                                                 <  2:3614475/93‑1 (MQ=255)
                                           gATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGt                                              <  1:2827837/93‑1 (MQ=255)
                                            aTGTGCGTTTAGCGCAGATGGAAGTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  1:1907331/93‑1 (MQ=255)
                                            aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  1:2885745/93‑1 (MQ=255)
                                            aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             >  1:4146509/1‑93 (MQ=255)
                                            aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  1:1492303/93‑1 (MQ=255)
                                            aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             >  2:4269508/1‑93 (MQ=255)
                                              gtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTg                                           <  2:355889/93‑1 (MQ=255)
                                               tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                          <  2:500039/93‑1 (MQ=255)
                                               tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                          >  1:451430/1‑93 (MQ=255)
                                                gCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGt                                         <  2:3634569/93‑1 (MQ=255)
                                                 cGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGtt                                        >  2:2182511/1‑93 (MQ=255)
                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       >  2:4535992/1‑93 (MQ=255)
                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGGCTGGTTg                                       >  2:4062748/1‑93 (MQ=255)
                                                   tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTgt                                      >  1:3945421/1‑92 (MQ=255)
                                                   tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      >  1:3756754/1‑93 (MQ=255)
                                                   tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      <  2:1287664/93‑1 (MQ=255)
                                                    ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  2:2786451/93‑1 (MQ=255)
                                                    ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  1:3777311/93‑1 (MQ=255)
                                                     tAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATg                                    >  2:167395/1‑93 (MQ=255)
                                                     tAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATg                                    <  2:3898784/93‑1 (MQ=255)
                                                      aGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATgg                                   >  1:3188844/1‑93 (MQ=255)
                                                       gcgcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgc                                  >  1:3535992/1‑93 (MQ=255)
                                                         gcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgcgc                                >  1:949649/1‑93 (MQ=255)
                                                          cAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCa                               >  2:2416326/1‑93 (MQ=255)
                                                          cAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCa                               >  2:77681/1‑93 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             <  2:2748083/93‑1 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             >  1:3037813/1‑93 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             >  2:4341861/1‑93 (MQ=255)
                                                              tGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGa                           >  1:2313737/1‑93 (MQ=255)
                                                              tGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGa                           <  2:4174776/93‑1 (MQ=255)
                                                                  aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                       <  1:1036045/93‑1 (MQ=255)
                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                      <  1:542625/93‑1 (MQ=255)
                                                                    tCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGatg                     <  2:2959140/93‑3 (MQ=255)
                                                                      aaCCCGTGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   <  2:2108856/93‑1 (MQ=255)
                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   <  2:1816464/93‑1 (MQ=255)
                                                                         ccGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTgg                >  1:4474955/1‑93 (MQ=255)
                                                                         ccGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTgg                >  1:2750359/1‑93 (MQ=255)
                                                                          cGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGt               <  2:1741389/93‑1 (MQ=255)
                                                                            ggcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTa             <  2:3795106/93‑1 (MQ=255)
                                                                             gcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTAt            <  1:3370989/93‑1 (MQ=255)
                                                                              cggcACATTAGAGCTGGGTGTGGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTAtt           <  1:4573687/93‑1 (MQ=255)
                                                                                 cacaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAc        <  1:154050/93‑1 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       <  2:1018982/93‑1 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       <  2:4452326/93‑1 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       <  1:1788039/93‑1 (MQ=255)
                                                                                   caTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActg      <  1:619267/93‑1 (MQ=255)
                                                                                    aTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgc     >  1:3826478/1‑93 (MQ=255)
                                                                                     ttAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgct    >  1:281634/1‑93 (MQ=255)
                                                                                       agagCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGACGCATCGAAAGACGGCTGGTTATTCACTGCTGa  <  1:3154028/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246098‑4246277

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