Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,474,881 G→A 16.1% intergenic (+25/‑103) yfdQ → / → yfdR CPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,474,8810GA16.1% 463.1 / 68.2 192intergenic (+25/‑103)yfdQ/yfdRCPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (94/67);  new base A (27/4);  total (121/71)
Fisher's exact test for biased strand distribution p-value = 2.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 202 total aligned reads displayed.

AATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAAC  >  NC_000913/2474783‑2474980
                                                                                                  |                                                                                                   
  tttCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgggg                                                                                                  <  2:138262/100‑1 (MQ=255)
   ttCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                                 <  1:453833/100‑1 (MQ=255)
   tgcGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAATCCTTTATTGGTACTTTCACCGCCTGATTTCATTACTTCAAATGCCCCTGCggggg                                                                                                 <  1:1840404/98‑1 (MQ=255)
      gTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCAt                                                                                              >  2:4009111/1‑100 (MQ=255)
      gTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCAt                                                                                              >  1:760002/1‑100 (MQ=255)
       tGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTGACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCAtt                                                                                             >  1:2995012/1‑100 (MQ=255)
        gATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCAttt                                                                                            >  1:1555997/1‑100 (MQ=255)
           cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATg                                                                                         <  1:2126567/100‑1 (MQ=255)
             gCTTGTTGAGAAATTCAAAGCCAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgcagt                                                                                       <  2:3699741/100‑6 (MQ=38)
             gCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGa                                                                                       >  2:3120500/1‑100 (MQ=255)
              cTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgcagta                                                                                      <  1:2189533/100‑7 (MQ=38)
                 gttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACg                                                                                   >  1:2046297/1‑100 (MQ=255)
                   tgcgaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTACGGGGGCATTTATGGAAACGTa                                                                                 <  1:1523417/97‑1 (MQ=255)
                    gagaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTaa                                                                                >  1:4156210/1‑100 (MQ=255)
                    gagaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTaa                                                                                >  2:548516/1‑100 (MQ=255)
                    gagaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTaa                                                                                >  2:3544136/1‑100 (MQ=255)
                    gagaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTaa                                                                                >  2:2083060/1‑100 (MQ=255)
                     agaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                               >  2:445246/1‑100 (MQ=255)
                     agaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                               >  2:1613424/1‑100 (MQ=255)
                     agaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                               <  2:1719494/100‑1 (MQ=255)
                      gaAATTCGAAGCCAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCACATGCCCCTGCGGGGGCATTTATGGAAACGTAAtt                                                                              >  1:2238231/1‑100 (MQ=255)
                      gaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTTATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACTTAAt                                                                               >  1:1849432/1‑99 (MQ=255)
                      gaAATTCAAAGACAGCAAAGTAGAAACCTTTAATGGTACTTTCACAGCCTGATTTCATTACTGCAATTGCCCCTGCGGGGGCATTTATGGACACGTAAtt                                                                              <  1:1421284/100‑1 (MQ=255)
                        aaTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTaa                                                                            >  1:699455/1‑100 (MQ=255)
                         aTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAAc                                                                           >  1:724947/1‑100 (MQ=255)
                         aTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAAc                                                                           <  1:3229830/100‑1 (MQ=255)
                          ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACt                                                                          <  2:256307/100‑1 (MQ=255)
                          ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACt                                                                          <  1:171585/100‑1 (MQ=255)
                            cAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCa                                                                        >  1:2425581/1‑100 (MQ=255)
                              aaGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaat                                                                      <  1:3787480/100‑1 (MQ=255)
                                  cAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                                  >  1:2744417/1‑100 (MQ=255)
                                  cAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                                  <  2:4259187/100‑1 (MQ=255)
                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                >  1:2389502/1‑100 (MQ=255)
                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                >  1:4076102/1‑100 (MQ=255)
                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                >  1:3794637/1‑100 (MQ=255)
                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                >  2:3249194/1‑100 (MQ=255)
                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                >  2:1452219/1‑100 (MQ=255)
                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAACCAc                                                                >  1:585732/1‑100 (MQ=255)
                                     cAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAcc                                                               <  1:3267549/100‑1 (MQ=255)
                                           aGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTCATTAACTCAATAATCACCGGATgg                                                         >  1:1990858/1‑100 (MQ=255)
                                           aGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATgg                                                         >  2:2407745/1‑100 (MQ=255)
                                                 cTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAggg                                                   <  1:1166934/100‑1 (MQ=255)
                                                 cTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAggg                                                   <  1:2542775/100‑1 (MQ=255)
                                                 cTTTATTGGTACGTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAggg                                                   >  2:2816440/1‑100 (MQ=255)
                                                    tATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCtt                                                <  2:4244243/100‑1 (MQ=255)
                                                    tATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCtt                                                >  1:2357353/1‑100 (MQ=255)
                                                       tGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                             <  2:3612074/100‑1 (MQ=255)
                                                        ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCtt                                            >  2:2138415/1‑100 (MQ=255)
                                                        ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCtt                                            >  1:673441/1‑100 (MQ=255)
                                                         gTACTTTCCGTGCCTGATATCATTACTGGAAATGCCCCTGCGGGGGCATTTATGAAAACGTAATTAACGCAATAATCACCCGATGGTGAGAGCTTCCgtt                                           <  1:1222043/100‑1 (MQ=255)
                                                            cTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAc                                        <  1:3839548/100‑1 (MQ=255)
                                                            cTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAc                                        >  1:3331563/1‑100 (MQ=255)
                                                             tttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                       <  2:1422772/100‑1 (MQ=255)
                                                             tttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                       <  2:4115474/100‑1 (MQ=255)
                                                                  ccGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAAc                                  <  2:634455/100‑1 (MQ=255)
                                                                  ccGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAAc                                  >  2:1933668/1‑100 (MQ=255)
                                                                   cGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAAc                                  <  1:2002847/99‑1 (MQ=255)
                                                                    gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTc                                >  2:1515639/1‑100 (MQ=255)
                                                                    gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTc                                >  2:2867774/1‑100 (MQ=255)
                                                                      cTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAg                              <  1:2838764/100‑1 (MQ=255)
                                                                      cTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAg                              <  1:2183864/100‑1 (MQ=255)
                                                                       tGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgc                             >  2:844397/1‑100 (MQ=255)
                                                                        gATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgcg                            <  1:3209574/100‑1 (MQ=255)
                                                                         aTTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgcgc                           <  1:258718/100‑1 (MQ=255)
                                                                          tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgcgc                           >  2:4169783/1‑99 (MQ=255)
                                                                           ttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCgg                         <  1:2861632/100‑1 (MQ=255)
                                                                                tACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                    <  2:916266/100‑1 (MQ=255)
                                                                                tACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGATTCCTTTTACCCAAACTCAGCGCGGTGCAg                    >  1:2695083/1‑100 (MQ=255)
                                                                                  ccataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg                  >  2:2301935/5‑100 (MQ=255)
                                                                                  ccataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg                  >  2:2261486/5‑100 (MQ=39)
                                                                                  ccataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg                  >  2:2956158/5‑100 (MQ=39)
                                                                                  cTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg                  <  2:2589879/100‑1 (MQ=255)
                                                                                  cTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                    <  1:3953240/98‑1 (MQ=255)
                                                                                   tGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:1037232/1‑100 (MQ=255)
                                                                                   cttaaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:2892986/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  1:1015527/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:510684/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  1:3535067/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:1256449/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:1787413/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:2725873/4‑100 (MQ=255)
                                                                                   cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc                 >  2:1710594/4‑100 (MQ=255)
                                                                                    gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                <  1:2611435/100‑1 (MQ=255)
                                                                                    gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                >  1:489955/1‑100 (MQ=255)
                                                                                    gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                >  1:3927663/1‑100 (MQ=255)
                                                                                    gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                >  2:173814/1‑100 (MQ=255)
                                                                                    gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                >  1:3079244/1‑100 (MQ=255)
                                                                                    gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTCTTACCCAAACTCAGCGCGGTGCAGCGCa                <  1:686439/100‑1 (MQ=255)
                                                                                    ataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                >  2:764217/3‑100 (MQ=255)
                                                                                    ataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa                >  2:3910385/3‑100 (MQ=255)
                                                                                       aaTGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCatat             >  2:1151218/1‑100 (MQ=255)
                                                                                       aaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCatat             >  2:1814101/1‑100 (MQ=255)
                                                                                          gCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACg          >  2:3085981/1‑100 (MQ=255)
                                                                                           ccccTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATTGTGAGGGCTTCCTTTTACCCAACCTCAGCGCGGTGCAGCGCATATACGt         >  1:2908901/1‑100 (MQ=255)
                                                                                             ccTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTgg       <  2:473057/100‑1 (MQ=255)
                                                                                              cTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTggn      <  2:3774837/100‑2 (MQ=255)
                                                                                               tGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGgag     >  2:75244/1‑100 (MQ=255)
                                                                                               tGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGgag     >  2:300578/1‑100 (MQ=255)
                                                                                               tGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGgag     <  2:1652971/100‑1 (MQ=255)
                                                                                                  gggggCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAAc  >  1:4029659/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAAC  >  NC_000913/2474783‑2474980

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