Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 998,433 A→G intergenic (+19/‑6) ycbR → / → ycbS putative chaperone/putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913998,4330AG99.3% 472.2 / NA 149intergenic (+19/‑6)ycbR/ycbSputative chaperone/putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (95/53);  total (96/53)
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 158 total aligned reads displayed.

AATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTTGTGTC  >  NC_000913/998343‑998521
                                                                                          |                                                                                        
aaTAGCCGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGc                                                                                        >  2:882323/1‑93 (MQ=255)
 aTAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCg                                                                                       <  2:6135063/93‑1 (MQ=255)
  tAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCga                                                                                      >  1:4170481/1‑93 (MQ=255)
  tAGCAGCGTAACGTATCGAACCATCAATCACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCga                                                                                      <  1:1163771/93‑1 (MQ=255)
    gcagcGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTAAAGAGGGCgaga                                                                                    <  2:1757938/93‑1 (MQ=255)
     cagcGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGAt                                                                                   <  1:4051695/93‑1 (MQ=255)
       gcGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGt                                                                                 <  2:2525902/93‑1 (MQ=255)
        cGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGta                                                                                >  2:5581803/1‑93 (MQ=255)
        cGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGta                                                                                >  1:2585231/1‑93 (MQ=255)
          tAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGtata                                                                              >  2:3443938/1‑93 (MQ=255)
          tAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGtata                                                                              <  2:3667037/93‑1 (MQ=255)
           aaCGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAACTCGTATGCCGTCTTCTGCTTGGAAAGCGAGATGTATAg                                                                             >  1:4023506/1‑93 (MQ=255)
           aaCGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAg                                                                             >  2:3225248/1‑93 (MQ=255)
            aCGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGa                                                                            >  1:5210158/1‑93 (MQ=255)
                aTCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTc                                                                        >  2:4921472/1‑93 (MQ=255)
                 tCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCa                                                                       >  1:3082469/1‑93 (MQ=255)
                 tCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCa                                                                       >  1:5981069/1‑93 (MQ=255)
                    aaCCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCg                                                                    >  2:3970121/1‑93 (MQ=255)
                    aaCCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCg                                                                    >  1:5790904/1‑93 (MQ=255)
                     aCCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGa                                                                   >  2:2662174/1‑93 (MQ=255)
                     aCCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGa                                                                   >  2:1214933/1‑93 (MQ=255)
                     aCCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGa                                                                   >  1:3697237/1‑93 (MQ=255)
                       cATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGaca                                                                 >  2:3800718/1‑93 (MQ=255)
                           aaTGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCTGCTAGGAGGGAGAGATGTATAGAACTCATCGACAacac                                                             >  1:1683378/1‑93 (MQ=255)
                            aTGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAAcaca                                                            >  1:4853130/1‑93 (MQ=255)
                               aCTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCACAAACACAGcc                                                         >  2:3946894/1‑93 (MQ=255)
                                cTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCt                                                        >  1:3834945/1‑93 (MQ=255)
                                cTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCt                                                        <  2:447204/93‑1 (MQ=255)
                                 aaCGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTg                                                       <  1:971694/92‑1 (MQ=255)
                                    ggggCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTa                                                    >  2:926633/1‑93 (MQ=255)
                                    ggggCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTa                                                    >  2:3837440/1‑93 (MQ=255)
                                     gggCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTaa                                                   <  1:2968985/93‑1 (MQ=255)
                                       gCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGc                                                 >  2:3668347/1‑93 (MQ=255)
                                       gCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGc                                                 <  2:6254697/93‑1 (MQ=255)
                                         ggAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCTGCTAGGAGGGCGAGATGTATAGAACTCATAGACAACACAGCCTGTTAAGctt                                               >  1:4238686/1‑92 (MQ=255)
                                         ggAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctc                                               >  1:1573725/1‑93 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCTGCTTGGAGGGCGAAATGCATAGAACTCATCGACAACACAGCCTGTTAAGctcg                                              >  1:5260086/1‑92 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctgt                                              >  1:2480421/1‑91 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctct                                              >  2:3507035/1‑93 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctct                                              >  1:2109884/1‑93 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctct                                              >  2:2831555/1‑93 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctct                                              <  1:4533429/93‑1 (MQ=255)
                                          gAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGctct                                              <  2:2452044/93‑1 (MQ=255)
                                           aaaCGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtg                                             <  2:414351/93‑1 (MQ=255)
                                            aaCGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCTGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtgg                                            >  1:1101525/1‑93 (MQ=255)
                                            aaCGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtgg                                            >  1:3503934/1‑93 (MQ=255)
                                             aCGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtggt                                           <  1:5268734/93‑1 (MQ=255)
                                              cgccgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtggtg                                          <  2:5452161/93‑1 (MQ=255)
                                               gccgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtggtgg                                         <  2:2216630/93‑1 (MQ=255)
                                               cccgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCtggtgg                                         <  1:2300920/92‑1 (MQ=255)
                                                ccgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGa                                        >  2:4286734/1‑93 (MQ=255)
                                                ccgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGa                                        >  2:5775261/1‑93 (MQ=255)
                                                ccgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGa                                        <  2:4726087/93‑1 (MQ=255)
                                                 cgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAg                                       >  1:5014100/1‑93 (MQ=255)
                                                 cgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAg                                       >  2:3386455/1‑93 (MQ=255)
                                                 cgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAg                                       >  2:1512028/1‑93 (MQ=255)
                                                 cgcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAg                                       >  1:1528629/1‑93 (MQ=255)
                                                  gcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGt                                      <  1:3565579/93‑1 (MQ=255)
                                                  gcAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGt                                      >  2:510260/1‑93 (MQ=255)
                                                   caactcaactGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTg                                     <  1:4987852/93‑1 (MQ=255)
                                                   caactcaactGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTg                                     >  2:3790465/1‑93 (MQ=255)
                                                    aactcaactGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGc                                    >  1:1282781/1‑93 (MQ=255)
                                                      ctcaactGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCa                                  >  1:788805/1‑93 (MQ=255)
                                                        taactGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATc                                <  1:3182432/92‑1 (MQ=255)
                                                         aactGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCg                               >  1:4972446/1‑93 (MQ=255)
                                                          actGTAAATCGTAAGCTGTCTGCAGTTAATAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGt                              <  2:1074962/93‑1 (MQ=255)
                                                          actGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGt                              <  2:806673/93‑1 (MQ=255)
                                                            tGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGttt                            >  1:270148/1‑93 (MQ=255)
                                                            tGTAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGttt                            >  2:2335918/1‑93 (MQ=255)
                                                              tAAATCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTAt                          >  1:1363935/1‑93 (MQ=255)
                                                                aaTCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATtg                        >  1:5029933/1‑93 (MQ=255)
                                                                 aTCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATtgg                       >  2:728959/1‑93 (MQ=255)
                                                                 aTCGTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGAGCCATCGTTTATtgg                       >  1:2355548/1‑93 (MQ=255)
                                                                    gTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATtggtgg                    <  1:4680575/93‑1 (MQ=255)
                                                                    gTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATtggtgg                    <  2:6198088/93‑1 (MQ=255)
                                                                    gTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATtggtgg                    <  1:5309941/93‑1 (MQ=255)
                                                                    gTAAGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATtggtgg                    >  1:866381/1‑93 (MQ=255)
                                                                     tAAGCCGTCTTCAGCTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTggg                   >  1:841103/1‑93 (MQ=255)
                                                                       aGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCt                 >  1:6185835/1‑93 (MQ=255)
                                                                       aGCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCt                 >  2:5329028/1‑93 (MQ=255)
                                                                        gCCGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtg                >  2:698884/1‑93 (MQ=255)
                                                                         ccGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtgg               >  1:4778766/1‑93 (MQ=255)
                                                                         ccGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtgg               <  1:5685440/93‑1 (MQ=255)
                                                                         ccGTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtgg               >  2:2932398/1‑93 (MQ=255)
                                                                           gTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtggtg             >  2:508073/1‑93 (MQ=255)
                                                                           gTCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtggtg             <  2:2972555/93‑1 (MQ=255)
                                                                            tatacAGTTAAGAGGGCGAGATTTATACAACTCTTCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCttgtgg            <  2:5720777/89‑1 (MQ=255)
                                                                            tCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtggtgg            >  1:4415465/1‑93 (MQ=255)
                                                                            tCTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCtggtgg            <  1:588933/93‑1 (MQ=255)
                                                                             cTTCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTggtggt           >  1:2691812/1‑93 (MQ=255)
                                                                              ttCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGgtggtg          >  1:4215338/1‑93 (MQ=255)
                                                                              ttCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGgtggtg          >  2:3023864/1‑93 (MQ=255)
                                                                               tCAGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGgtgt         <  1:5151655/93‑1 (MQ=255)
                                                                                 aGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGttt       >  1:6267662/1‑93 (MQ=255)
                                                                                 aGTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGttt       >  1:5648793/1‑93 (MQ=255)
                                                                                  gTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTtg      <  1:2128708/93‑1 (MQ=255)
                                                                                  gTTAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTtg      <  2:2630064/93‑1 (MQ=255)
                                                                                   ttAAGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTtgt     >  1:2784164/1‑93 (MQ=255)
                                                                                     aaGAGGGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTTgtgt   >  2:4491836/1‑93 (MQ=255)
                                                                                      agagggCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTTGTGTc  >  2:4644643/1‑93 (MQ=255)
                                                                                          |                                                                                        
AATAGCAGCGTAACGTATCGAACCATCAATGACTACGGGGCGGAAACGCCGCAACTCAACTGTAAATCGTAAGCCGTCTTCAGTTAAGAGAGCGAGATGTATAGAACTCATCGACAACACAGCCTGTTAAGCTCTGGTGGAGTGCCATCGTTTATTGGTGGGCTGGTGGTGTTTGTGTC  >  NC_000913/998343‑998521

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