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.7% 977.4 / NA 298E213E (GAA→GAGlamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (153/143);  total (153/145)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246098‑4246277
                                                                                          |                                                                                         
caaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCAGGGCGGCACATTAGAGCt                                                                                         >  2:759818/1‑93 (MQ=255)
 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                        >  2:5339695/1‑93 (MQ=255)
 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                        >  1:1836410/1‑93 (MQ=255)
 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGCGCGGCACATTAGAGCTg                                                                                        >  1:2408278/1‑93 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                       >  1:5081347/1‑93 (MQ=255)
   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      >  1:3317252/1‑93 (MQ=255)
    atTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgt                                                                                     <  2:1661505/93‑1 (MQ=255)
      ttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtgt                                                                                   >  1:4556704/1‑93 (MQ=255)
      ttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtgt                                                                                   <  2:2503330/93‑1 (MQ=255)
       tATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTc                                                                                  <  1:2374897/93‑1 (MQ=255)
        aTGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCg                                                                                 <  2:5265633/93‑1 (MQ=255)
          gACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGAc                                                                               >  1:3910023/1‑93 (MQ=255)
           aCTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACt                                                                              >  1:4839259/1‑93 (MQ=255)
           aCTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACt                                                                              >  1:1554101/1‑93 (MQ=255)
            cTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTa                                                                             <  1:2121202/93‑1 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            <  2:6598271/93‑1 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            <  2:4780087/93‑1 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            >  1:309892/1‑93 (MQ=255)
              ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                           >  1:2623333/1‑93 (MQ=255)
              ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                           >  2:4073698/1‑93 (MQ=255)
                aCCAACGAAACCGCGAACGACGTTTTCGCTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTCCGGt                                                                         >  1:4211878/1‑93 (MQ=255)
                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACggc                                                                         >  1:221076/1‑92 (MQ=255)
                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                         <  1:1746067/93‑1 (MQ=255)
                 ccAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        >  1:2933598/1‑93 (MQ=255)
                  cAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCg                                                                       <  2:5586880/93‑1 (MQ=255)
                  cAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCg                                                                       <  1:96410/93‑1 (MQ=255)
                   caCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                      <  2:6437940/92‑1 (MQ=255)
                   aaCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                      >  2:3061102/1‑93 (MQ=255)
                    aCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTg                                                                     >  2:1242875/1‑93 (MQ=255)
                    aCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTg                                                                     >  1:1305961/1‑93 (MQ=255)
                    aCGAAACCGCGAACGACGTTTTCGATGCGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTg                                                                     >  1:4328165/1‑93 (MQ=255)
                     cGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGc                                                                    >  2:6507191/1‑93 (MQ=255)
                     cGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGc                                                                    <  2:4156766/93‑1 (MQ=255)
                      gAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGcc                                                                   <  2:5661381/93‑1 (MQ=255)
                        aaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGAAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCaa                                                                 >  2:2382050/1‑93 (MQ=255)
                         aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAAc                                                                >  1:5747736/1‑93 (MQ=255)
                          ccGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               <  2:4851977/93‑1 (MQ=255)
                          ccGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               <  2:5185207/93‑1 (MQ=255)
                          ccGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               >  2:1831162/1‑93 (MQ=255)
                           cgcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACtt                                                              <  2:5035431/93‑1 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             >  1:4127065/1‑93 (MQ=255)
                             cgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGc                                                            <  1:6013902/93‑1 (MQ=255)
                             cgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACAGTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGc                                                            >  2:3359486/1‑93 (MQ=255)
                              gAACGACGTTTTCGATGTGCGTTTAGCGCCGATGGAAATCAACCCGGGCGGCACAGTAGAGCTGGGTGTCGACTACGGCCGCGCCAACTTGCg                                                           >  2:1334004/1‑93 (MQ=255)
                              gAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCg                                                           >  1:1065816/1‑93 (MQ=255)
                              gAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCg                                                           >  1:3521926/1‑93 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          <  2:3844369/93‑1 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  1:2806487/1‑93 (MQ=255)
                                acgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTg                                                         >  1:4766388/1‑93 (MQ=255)
                                 cgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGa                                                        <  2:4543494/93‑1 (MQ=255)
                                     gTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAAc                                                    <  1:2016706/93‑1 (MQ=255)
                                      ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                   >  1:3814864/1‑93 (MQ=255)
                                       tttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTa                                                  >  2:3505606/1‑93 (MQ=255)
                                        ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTAt                                                 <  1:1913501/93‑1 (MQ=255)
                                         tCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATc                                                <  2:4201151/93‑1 (MQ=255)
                                           gATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGt                                              >  2:6580651/1‑93 (MQ=255)
                                           gATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGt                                              >  1:6417454/1‑93 (MQ=255)
                                             tgtgCGTTTAGCGCAGATGGAAATCAGCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                            <  2:540734/93‑1 (MQ=255)
                                              gtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTg                                           <  2:5773392/93‑1 (MQ=255)
                                                gCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGt                                         <  2:594063/93‑1 (MQ=255)
                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       <  1:2180801/93‑1 (MQ=255)
                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       >  2:2775004/1‑93 (MQ=255)
                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       >  1:1517237/1‑93 (MQ=255)
                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       >  1:6486289/1‑93 (MQ=255)
                                                   tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      >  1:6205253/1‑93 (MQ=255)
                                                    ttAGCGCAGATGGAAATCAACCCGNGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  1:688913/93‑1 (MQ=255)
                                                    ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  1:4488603/93‑1 (MQ=255)
                                                       gcgcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgc                                  <  1:141292/93‑1 (MQ=255)
                                                         gcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgcgc                                >  1:1239066/1‑93 (MQ=255)
                                                          gaGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCa                               <  1:1993248/92‑1 (MQ=255)
                                                          cAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCa                               <  1:5325071/93‑1 (MQ=255)
                                                           aGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCAt                              >  1:2738500/1‑93 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             <  2:684455/93‑1 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             <  1:5688160/93‑1 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             >  1:869558/1‑93 (MQ=255)
                                                                gAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaaa                         <  2:3723406/93‑1 (MQ=255)
                                                                gAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaaa                         <  2:4941426/93‑1 (MQ=255)
                                                                 aaaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAg                        >  1:1354733/1‑93 (MQ=255)
                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                      <  2:28887/93‑1 (MQ=255)
                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                      <  1:376152/93‑1 (MQ=255)
                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   <  2:3634702/93‑1 (MQ=255)
                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   >  2:1149225/1‑93 (MQ=255)
                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   >  2:2592212/1‑93 (MQ=255)
                                                                        cccGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTg                 <  1:5355931/93‑1 (MQ=255)
                                                                          cGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGt               >  1:516378/1‑93 (MQ=255)
                                                                          cGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGt               <  2:2636650/93‑1 (MQ=255)
                                                                          cGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGt               <  1:5494737/93‑1 (MQ=255)
                                                                            ggcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTa             >  2:5283237/1‑93 (MQ=255)
                                                                             gcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTAt            <  2:3007279/93‑1 (MQ=255)
                                                                              cggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTAtt           <  1:908178/93‑1 (MQ=255)
                                                                                gcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCa         <  2:2302171/93‑1 (MQ=255)
                                                                                gcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCa         <  2:6195613/93‑1 (MQ=255)
                                                                                 cacaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAc        <  1:5803862/93‑1 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       <  1:4966991/93‑1 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       >  2:6259823/1‑93 (MQ=255)
                                                                                   caTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActg      >  1:4918236/1‑93 (MQ=255)
                                                                                     ttAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgct    >  2:2093692/1‑93 (MQ=255)
                                                                                      tAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgctg   <  2:5705894/93‑1 (MQ=255)
                                                                                      tAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgctg   <  2:631366/93‑1 (MQ=255)
                                                                                       agagCTGGGTGTCGACTACGGTCGTGCCAACTTGGGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGa  >  2:5480251/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246098‑4246277

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