Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,931,803 T→C V250A (GTA→GCA)  ygdG → flap endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,931,8030TC99.7% 1189.0 / NA 282V250A (GTA→GCA) ygdGflap endonuclease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (137/144);  total (137/145)
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 282 total aligned reads displayed.

ATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTC  >  NC_000913/2931705‑2931900
                                                                                                  |                                                                                                 
aTTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCAc                                                                                                 >  1:4649391/1‑101 (MQ=255)
  tAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACgg                                                                                               >  2:4649394/1‑101 (MQ=255)
   aGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTCGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGt                                                                                              >  1:4649396/1‑101 (MQ=255)
     aaaCCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTaa                                                                                            <  1:4649145/101‑1 (MQ=255)
     aaaCCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTaa                                                                                            >  1:4649400/1‑101 (MQ=255)
        cccATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACgg                                                                                         >  2:4649405/1‑101 (MQ=255)
        cccATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACgg                                                                                         >  1:4649404/1‑101 (MQ=255)
         ccATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGc                                                                                        >  2:4649408/1‑101 (MQ=255)
          cATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCg                                                                                       >  1:4649411/1‑101 (MQ=255)
           aTAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      <  1:4649024/101‑1 (MQ=255)
           aTAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649415/1‑101 (MQ=255)
            tAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649417/1‑100 (MQ=255)
            tAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAg                                                                                     <  2:4649207/101‑1 (MQ=255)
            tAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAg                                                                                     <  1:4649080/101‑1 (MQ=255)
               agagATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649422/1‑97 (MQ=255)
               agagATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649421/1‑97 (MQ=255)
                 agaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649427/1‑95 (MQ=255)
                  gaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAAc                                                                                           >  1:4649432/1‑89 (MQ=255)
                  gaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649430/1‑94 (MQ=255)
                    tGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATAc                                                                             <  2:4648912/101‑1 (MQ=255)
                      gCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATAcgc                                                                           <  1:4648643/101‑1 (MQ=255)
                         tttCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649437/1‑87 (MQ=255)
                          ttCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649438/1‑86 (MQ=255)
                             tgtgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649442/1‑83 (MQ=255)
                             tgtgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAAc                                                                    <  2:4649044/101‑1 (MQ=255)
                             tgtgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAAc                                                                    <  1:4648874/101‑1 (MQ=255)
                             tgtgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAAc                                                                    <  1:4648938/101‑1 (MQ=255)
                              gtgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAAcg                                                                   <  2:4649164/101‑1 (MQ=255)
                               tgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCg                                                                                       >  2:4649445/1‑80 (MQ=255)
                               tgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAAcgt                                                                  <  2:4648898/101‑1 (MQ=255)
                                gCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649448/1‑80 (MQ=255)
                                 ccGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649450/1‑79 (MQ=255)
                                  cgcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649454/1‑78 (MQ=255)
                                   gcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649457/1‑77 (MQ=255)
                                   gcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTa                                                              <  1:4648923/101‑1 (MQ=255)
                                     gATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAg                                                                                     >  1:4649460/1‑76 (MQ=255)
                                     gATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATc                                                            <  2:4648997/101‑1 (MQ=255)
                                       tatTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649463/1‑73 (MQ=255)
                                       tatTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCg                                                          <  1:4648890/101‑1 (MQ=255)
                                        atTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649465/1‑72 (MQ=255)
                                          tGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCg                                                                                       >  1:4649468/1‑69 (MQ=255)
                                           gCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCAt                                                      <  1:4649033/101‑1 (MQ=255)
                                            cccGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649472/1‑68 (MQ=255)
                                            cccGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649473/1‑68 (MQ=255)
                                              cGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649477/1‑66 (MQ=255)
                                              cGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTAt                                                   <  1:4649133/101‑1 (MQ=255)
                                              cGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTAt                                                   <  1:4649140/101‑1 (MQ=255)
                                               gCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649478/1‑65 (MQ=255)
                                               gCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649483/1‑65 (MQ=255)
                                                cTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATca                                                 <  2:4649003/101‑1 (MQ=255)
                                                 ttACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649486/1‑63 (MQ=255)
                                                 ttACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATcac                                                <  2:4649180/101‑1 (MQ=255)
                                                  tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACgg                                                                                         >  2:4649492/1‑59 (MQ=255)
                                                  tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCg                                                                                       >  1:4649491/1‑61 (MQ=255)
                                                  tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649488/1‑62 (MQ=255)
                                                  tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATcaca                                               <  2:4648975/101‑1 (MQ=255)
                                                     aaaCCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCa                                            <  2:4649248/101‑1 (MQ=255)
                                                      aaCCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649496/1‑58 (MQ=255)
                                                       aCCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649501/1‑57 (MQ=255)
                                                       aCCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649500/1‑57 (MQ=255)
                                                        ccGATTTTCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTAt                                                   <  1:4649077/91‑1 (MQ=255)
                                                        ccGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAgcg                                         <  2:4648803/101‑1 (MQ=255)
                                                          gATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649506/1‑54 (MQ=255)
                                                          gATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCa                                       <  2:4649015/101‑1 (MQ=255)
                                                           aTTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649507/1‑53 (MQ=255)
                                                           aTTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCAt                                      <  2:4648892/101‑1 (MQ=255)
                                                             ttGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCAttt                                    <  1:4649111/101‑1 (MQ=255)
                                                              tGCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTa                                   <  1:4649334/101‑1 (MQ=255)
                                                               gCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCa                                                                       <  1:4648958/64‑1 (MQ=255)
                                                               gCATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTAc                                  <  1:4648969/101‑1 (MQ=255)
                                                                cATATCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGc                                                                                        >  2:4649511/1‑46 (MQ=255)
                                                                 atatCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649514/1‑47 (MQ=255)
                                                                 atatCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCa                                <  2:4649056/101‑1 (MQ=255)
                                                                  tatCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649515/1‑46 (MQ=255)
                                                                    tCGACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGc                             <  2:4649064/101‑1 (MQ=255)
                                                                      gACGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649520/1‑42 (MQ=255)
                                                                       aCGGCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGt                          <  1:4649095/101‑1 (MQ=255)
                                                                         ggCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649523/1‑39 (MQ=39)
                                                                         ggCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTAt                        <  2:4648866/101‑1 (MQ=255)
                                                                          gCAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATg                       <  2:4649029/101‑1 (MQ=255)
                                                                           cAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649525/1‑37 (MQ=38)
                                                                           cAATTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATgg                      <  1:4648941/101‑1 (MQ=255)
                                                                            aaTTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  2:4649529/1‑36 (MQ=38)
                                                                            aaTTTACAGCAATTGCGGTTGGCACGGTAACGGCGa                                                                                      >  1:4649527/1‑36 (MQ=38)
                                                                            aaTTTACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGt                     <  1:4649156/101‑1 (MQ=255)
                                                                              tttACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTaa                   <  1:4648814/101‑1 (MQ=255)
                                                                               ttACAGCAATTGCGGTTGGCGCGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTaaa                  <  1:4648988/101‑1 (MQ=255)
                                                                               ttACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTa                                   <  1:4649012/84‑1 (MQ=255)
                                                                               ttACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTaaa                  <  2:4649104/101‑1 (MQ=255)
                                                                                tACAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAg                 <  1:4649057/101‑1 (MQ=255)
                                                                                  cAGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGct               <  2:4648836/101‑1 (MQ=255)
                                                                                   aGCAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGctc              <  1:4649120/101‑1 (MQ=255)
                                                                                     cAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTa            <  1:4649185/101‑1 (MQ=255)
                                                                                     cAATTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTa            <  2:4649189/101‑1 (MQ=255)
                                                                                       aTTGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACa          <  2:4649086/101‑1 (MQ=255)
                                                                                         tGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAat        <  2:4648925/101‑1 (MQ=255)
                                                                                         tGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAat        <  2:4648952/101‑1 (MQ=255)
                                                                                         tGCGGTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAat        <  1:4649004/101‑1 (MQ=255)
                                                                                             gTTGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATca                                                 <  1:4649046/56‑1 (MQ=255)
                                                                                               tGGCACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTc  <  1:4649269/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
ATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTC  >  NC_000913/2931705‑2931900

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