Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,474,878 T→C 17.1% intergenic (+22/‑106) 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,8780TC17.1% 263.3 / 46.9 122intergenic (+22/‑106)yfdQ/yfdRCPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (54/47);  new base C (10/11);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 123 total aligned reads displayed.

AACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGA  >  NC_000913/2474779‑2474976
                                                                                                   |                                                                                                  
aaCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCt                                                                                                    >  1:837807/1‑100 (MQ=255)
 aCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTg                                                                                                   <  2:679744/100‑1 (MQ=255)
   gAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCg                                                                                                 >  2:1001316/1‑100 (MQ=255)
     aTTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggg                                                                                               >  2:1832629/1‑100 (MQ=255)
      tttCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgggg                                                                                              <  2:1699177/100‑1 (MQ=255)
      tttCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgggg                                                                                              <  2:424633/100‑1 (MQ=255)
       ttCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                             <  2:1460208/100‑1 (MQ=255)
       tccgatctCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                             <  2:2562728/93‑1 (MQ=255)
        tCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGc                                                                                            >  2:2789753/1‑100 (MQ=255)
        tCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGc                                                                                            >  2:2733950/1‑100 (MQ=255)
           tGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGACACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCAtt                                                                                         <  1:1672434/100‑1 (MQ=255)
           tGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtt                                                                                         <  1:544218/100‑1 (MQ=255)
             aTCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttg                                                                                       <  1:1256495/100‑2 (MQ=255)
              tCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTAt                                                                                      >  2:1778954/1‑100 (MQ=255)
              tCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgc                                                                                      <  2:2462720/100‑3 (MQ=255)
               cagCTTGTTGTGAAATTTATAGCCAGCAAAGTAGAAACCTTTATTTGTACTTTCACCGCCTCATTTCATTGCTGGTAATGCCCCTGCAGGGGCAtttgc                                                                                      <  1:1418365/97‑3 (MQ=255)
               cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGACACCTTTATTGGTACTTTCACCGCCTGATTTCATCACTGCAAATGCCCCTGCGCGGACATTTCTg                                                                                     >  1:2004947/1‑100 (MQ=255)
               cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgggggagatcgga                                                                                     >  1:2562713/1‑92 (MQ=255)
               cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATg                                                                                     <  2:938080/100‑1 (MQ=255)
               cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATg                                                                                     <  2:837810/100‑1 (MQ=255)
               cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATg                                                                                     <  1:228166/100‑1 (MQ=255)
                 gCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAAGCTTTATTGGTACGGTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgcag                                                                                    <  1:2859964/99‑5 (MQ=38)
                 gCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGa                                                                                   <  2:2378587/100‑1 (MQ=255)
                 gCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgcagt                                                                                   <  2:178148/100‑6 (MQ=39)
                 gCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgcagt                                                                                   <  2:1825939/100‑6 (MQ=39)
                   ttgttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGaaa                                                                                 >  2:1222885/1‑100 (MQ=255)
                   ttgttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtttgcagtaa                                                                                 <  1:510071/100‑8 (MQ=38)
                     gttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACg                                                                               <  1:1952385/100‑1 (MQ=255)
                      ttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGt                                                                              <  2:1102933/100‑1 (MQ=255)
                       tgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTa                                                                             >  2:2590527/1‑100 (MQ=255)
                        gagaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTaa                                                                            >  2:485780/1‑100 (MQ=255)
                        gagaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTaa                                                                            >  1:1262380/1‑100 (MQ=255)
                         agaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                           <  1:129117/100‑1 (MQ=255)
                          gaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAtt                                                                          >  2:189310/1‑100 (MQ=255)
                           aaaTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTa                                                                         <  1:2864704/100‑1 (MQ=255)
                           aaaTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTa                                                                         <  2:556387/100‑1 (MQ=255)
                            aaTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTaa                                                                        >  2:2663559/1‑100 (MQ=255)
                                cAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCa                                                                    >  2:2335018/1‑100 (MQ=255)
                                cAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCa                                                                    >  2:961733/1‑100 (MQ=255)
                                 aaaGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaa                                                                   <  2:626346/100‑1 (MQ=255)
                                    gACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaataa                                                                >  1:1483507/1‑100 (MQ=255)
                                    gACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaataa                                                                <  2:1861310/100‑1 (MQ=255)
                                      cAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                              <  2:366847/100‑1 (MQ=255)
                                        gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                            >  2:1118835/1‑100 (MQ=255)
                                        gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                            >  1:2073308/1‑100 (MQ=255)
                                        gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                            >  2:1477803/1‑100 (MQ=255)
                                        gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                            >  1:2790466/1‑100 (MQ=255)
                                         cAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                            >  1:680858/1‑99 (MQ=255)
                                                   aCCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAg                                                 <  2:1061546/100‑1 (MQ=255)
                                                    ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAgg                                                <  1:1284884/100‑1 (MQ=255)
                                                    ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAgg                                                >  1:595286/1‑100 (MQ=255)
                                                     cTTTATTGGTACTTTCACCGACTGATTTCATTACTGCAAATGCTCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAggg                                               >  1:1128116/1‑100 (MQ=255)
                                                       ttATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCATGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCt                                             <  1:665117/100‑1 (MQ=255)
                                                         aTTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTc                                           <  1:1487326/100‑1 (MQ=255)
                                                          ttGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTcc                                          >  1:891286/1‑100 (MQ=255)
                                                           tGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                         >  2:1257776/1‑100 (MQ=255)
                                                            ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCtt                                        <  1:218642/100‑1 (MQ=255)
                                                            ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCtt                                        >  2:492701/1‑100 (MQ=255)
                                                            ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCtt                                        >  2:2269873/1‑100 (MQ=255)
                                                             gTACCTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATCGTGAGGGCTTCCttt                                       <  1:192328/100‑1 (MQ=255)
                                                                cTTTCACCGCCTGATTTCATTACTGCAAAGGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAc                                    <  2:1177851/100‑1 (MQ=255)
                                                                 tttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                   <  2:2100594/100‑1 (MQ=255)
                                                                 tttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                   <  2:1882826/100‑1 (MQ=255)
                                                                 tttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                   <  1:783647/100‑1 (MQ=255)
                                                                 tttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACCCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                   >  1:636645/1‑100 (MQ=255)
                                                                   tCACCTGCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGTTGAGGGCTTCCTTTTAccc                                  <  1:2495135/99‑1 (MQ=255)
                                                                     aCCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCaaa                               >  2:2018662/1‑100 (MQ=255)
                                                                       cGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACt                             >  2:208351/1‑100 (MQ=255)
                                                                        gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTc                            <  2:1794274/100‑1 (MQ=255)
                                                                        gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTc                            >  2:717666/1‑100 (MQ=255)
                                                                              tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGcgcg                      >  1:2033910/1‑100 (MQ=255)
                                                                              tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGcgcg                      >  2:290238/1‑100 (MQ=255)
                                                                               ttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCgg                     <  1:1222880/100‑1 (MQ=255)
                                                                                 cATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTg                   <  2:2150650/100‑1 (MQ=255)
                                                                                  aTTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTACTTAACTCAATAATCACCGGATGGTGAGAGCTTCCTTTTACCCAAACTCAGCGCGGt                    >  1:600850/1‑98 (MQ=255)
                                                                                   tttccataaaTGCAACCGCAGGGACATTTGTGGAAACGTAATTAACTCAATAATCCACGGATGGTCAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCa                 >  1:2744867/8‑100 (MQ=255)
                                                                                    taatGCAAATGCCCCTGCGGGGGCATTGATGAATACGTAATTAACTCAATTATCACCAGATGCAAAGGGCAGCCCTAGACGCAAACTCAGAGAGGTGCAg                >  1:2557027/4‑100 (MQ=255)
                                                                                     tccataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgc               >  2:736652/6‑100 (MQ=39)
                                                                                     aCTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgc               >  2:2860637/1‑100 (MQ=255)
                                                                                      cTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATCATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCgg                     >  1:2726953/1‑93 (MQ=255)
                                                                                      cTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg              >  2:381101/1‑100 (MQ=255)
                                                                                      cTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg              >  2:162653/1‑100 (MQ=255)
                                                                                      cTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcg              >  2:1549845/1‑100 (MQ=255)
                                                                                       cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTCCCCAAACTCAGCGCGGTGCAgcgc             >  1:848455/4‑100 (MQ=255)
                                                                                       cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCACACTCAGCGCGGTGCAgcgc             >  1:2563967/4‑100 (MQ=255)
                                                                                       cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc             >  2:951890/4‑100 (MQ=255)
                                                                                       cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc             >  1:1151723/4‑100 (MQ=255)
                                                                                       cataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAgcgc             >  1:104552/4‑100 (MQ=255)
                                                                                        gCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa            <  2:592815/100‑1 (MQ=255)
                                                                                        ataaaTGCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCa            >  1:10590/3‑100 (MQ=255)
                                                                                         cAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCat           <  1:2273969/100‑1 (MQ=255)
                                                                                         cAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCat           <  2:1698787/100‑1 (MQ=255)
                                                                                           aaTGCCCCTGCGGGGGCCTTTATGGACCCCTAATTAACTCACAAATCACCAGATGGTAAAGGCTAACTTGTACCGACACGCAGcgcgccgacttcc             >  1:2616987/1‑87 (MQ=255)
                                                                                             tGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATAc       >  2:1234708/1‑100 (MQ=255)
                                                                                             tGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATAc       >  2:1226701/1‑100 (MQ=255)
                                                                                              gCCCCCGCAGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACg      >  2:100066/1‑100 (MQ=255)
                                                                                               ccccTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGt     <  2:1942574/100‑1 (MQ=255)
                                                                                               ccccTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGt     >  1:259743/1‑100 (MQ=255)
                                                                                                 ccTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTgg   <  1:2208472/100‑1 (MQ=255)
                                                                                                 ccTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTgg   <  2:1262385/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGA  >  NC_000913/2474779‑2474976

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