Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,013,772 +C intergenic (‑53/‑246) ysgA ← / → udp hypothetical protein/uridine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,013,7701.C99.2% 781.8 / ‑2.0 239intergenic (‑51/‑248)ysgA/udphypothetical protein/uridine phosphorylase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (123/114);  total (124/115)
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 267 total aligned reads displayed.

GCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTG  >  NC_000913/4013679‑4013857
                                                                                              |                                                                                       
gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                                                                          <  2:496419/93‑1 (MQ=255)
 cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAc                                                                                        <  2:1686892/93‑1 (MQ=255)
  gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc                                                                                       <  2:5914123/93‑1 (MQ=255)
  gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc                                                                                       >  1:2535241/1‑93 (MQ=255)
  gtcgtcCGGGGTCTGAACGATGGGCGAAGCGAGAGGCGATGCAGCAGGGGCAAATCCAGATTGTTGTGTTGTTGCCATGGGATTTTCCGGAcc                                                                                       >  1:101316/1‑93 (MQ=255)
    cgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccct                                                                                     <  1:1655825/93‑3 (MQ=255)
    cgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccct                                                                                     <  2:4716581/93‑3 (MQ=255)
     gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt                                                                                    <  1:644660/93‑4 (MQ=255)
      tcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctta                                                                                   <  1:1894735/93‑5 (MQ=255)
      tcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctta                                                                                   <  1:348145/93‑5 (MQ=255)
        cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTGCGTACCCTtata                                                                                  >  2:4326285/1‑93 (MQ=255)
        cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata                                                                                  >  1:543048/1‑93 (MQ=255)
          gggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCCAATCCAGATTGTTGCGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa                                                                                >  2:3379471/1‑93 (MQ=255)
           ggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa                                                                               >  2:4241568/1‑93 (MQ=255)
            gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                                              <  1:2926599/93‑1 (MQ=255)
            gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                                              >  1:3869899/1‑93 (MQ=255)
            gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                                              <  2:5866842/93‑1 (MQ=255)
            gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                                              <  2:5402621/93‑1 (MQ=255)
             tCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTTCCATGGTATTCTCCGTACCCTTATAAAAAt                                                                             <  2:5620132/93‑1 (MQ=255)
              cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg                                                                            <  2:2626422/93‑1 (MQ=255)
               tGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGt                                                                           >  2:2996493/1‑93 (MQ=255)
                gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGtt                                                                          >  2:2350480/1‑93 (MQ=255)
                gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTATCCGTACCCTTATAAAAATGtt                                                                          >  2:3135134/1‑93 (MQ=255)
                 aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCTTACCCTTATACAAATGTTg                                                                         <  1:6521471/93‑1 (MQ=255)
                  aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgc                                                                        <  1:1078415/93‑1 (MQ=255)
                  aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgc                                                                        <  2:3912716/93‑1 (MQ=255)
                  aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgc                                                                        >  1:2771640/1‑93 (MQ=255)
                    gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgcgc                                                                      >  1:1958790/1‑93 (MQ=255)
                      tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCaa                                                                    <  1:844795/93‑1 (MQ=255)
                       ggTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAAt                                                                   <  1:6347980/93‑1 (MQ=255)
                          cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt                                                                <  2:103879/93‑1 (MQ=255)
                           gAAGCGAGAGGCGATGCAGCAGGTGCACATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTCCCCTTATAAAAATGTTGCGCAATGTTa                                                               <  1:3752282/93‑1 (MQ=255)
                           gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa                                                               >  2:168113/1‑93 (MQ=255)
                           gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa                                                               <  2:979924/93‑1 (MQ=255)
                            aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa                                                              <  1:2729266/93‑1 (MQ=255)
                            aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa                                                              >  1:6063114/1‑93 (MQ=255)
                             aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAAc                                                             <  2:5157815/93‑1 (MQ=255)
                              gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt                                                            >  2:1324402/1‑93 (MQ=255)
                              gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt                                                            <  2:6652698/93‑1 (MQ=255)
                              gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt                                                            >  1:4669633/1‑93 (MQ=255)
                               cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta                                                           <  1:5228500/93‑1 (MQ=255)
                               cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta                                                           <  2:6010213/93‑1 (MQ=255)
                                 agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtata                                                         <  2:3329705/93‑1 (MQ=255)
                                    ggcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTc                                                      >  1:4240871/5‑93 (MQ=255)
                                    ggcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTc                                                      <  2:224120/89‑1 (MQ=255)
                                     gcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCa                                                     <  2:707998/90‑1 (MQ=255)
                                     gcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTAAAGTCa                                                     >  1:5322435/4‑93 (MQ=255)
                                      cgctGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAg                                                    <  1:5540291/90‑1 (MQ=255)
                                       gaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGc                                                   >  2:6154035/2‑93 (MQ=255)
                                        aTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCa                                                  <  1:4572031/93‑1 (MQ=255)
                                         tGCAGCAGGTGCAAATCCGGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                                                 <  2:4807616/93‑1 (MQ=255)
                                         tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                                                 <  1:709419/93‑1 (MQ=255)
                                         tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                                                 <  1:5695989/93‑1 (MQ=255)
                                         tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                                                 <  2:136053/93‑1 (MQ=255)
                                          gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                                                <  2:418377/93‑1 (MQ=255)
                                          gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                                                <  1:2412701/93‑1 (MQ=255)
                                            agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGca                                              <  2:5427997/93‑1 (MQ=255)
                                            agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGca                                              >  1:6273813/1‑93 (MQ=255)
                                            agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGca                                              >  1:1461814/1‑93 (MQ=255)
                                             gcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaa                                             >  1:4077316/1‑93 (MQ=255)
                                               agGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaaca                                           >  2:4029324/1‑93 (MQ=255)
                                                    cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATca                                      <  1:4138868/93‑1 (MQ=255)
                                                    cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATca                                      >  1:31356/1‑93 (MQ=255)
                                                    cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATca                                      <  2:5089027/93‑1 (MQ=255)
                                                     aaaTCCAGCTTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATcac                                     >  1:6572512/1‑93 (MQ=255)
                                                     aaaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATcac                                     >  1:4266254/1‑93 (MQ=255)
                                                       aTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAt                                   >  2:578272/1‑93 (MQ=255)
                                                         ccAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTg                                 <  1:3071762/93‑1 (MQ=255)
                                                            gATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCt                              >  2:5762022/1‑93 (MQ=255)
                                                            gATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCt                              <  2:5555913/93‑1 (MQ=255)
                                                                gttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAAt                          <  2:3624241/93‑1 (MQ=255)
                                                                gttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAAt                          <  1:3843270/93‑1 (MQ=255)
                                                                 ttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc                         >  1:5084556/1‑93 (MQ=255)
                                                                 ttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc                         <  2:2565019/93‑1 (MQ=255)
                                                                 ttgggttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGCTGCGCAATGCTAACTATAGTCAGCATGCAACAAATCACATTGTCTGAATc                         >  2:2192897/1‑93 (MQ=255)
                                                                  tgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCg                        <  2:5939586/93‑1 (MQ=255)
                                                                   gtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCgg                       >  1:5863333/1‑93 (MQ=255)
                                                                      ttgttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTc                    >  1:228409/1‑93 (MQ=255)
                                                                        gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCAt                  >  2:2511026/1‑93 (MQ=255)
                                                                        gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCAt                  >  2:4550580/1‑93 (MQ=255)
                                                                        gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCAt                  <  1:212947/93‑1 (MQ=255)
                                                                          tgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCt                >  2:3733434/1‑93 (MQ=255)
                                                                           gCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCtt               >  2:348442/1‑93 (MQ=255)
                                                                            ccATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCttt              >  2:1161225/1‑93 (MQ=255)
                                                                             cATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCtttt             >  2:4351087/1‑93 (MQ=255)
                                                                              ctGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTa            <  1:3611467/92‑1 (MQ=255)
                                                                               tGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTAt           >  1:4905978/1‑93 (MQ=255)
                                                                               tGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTAt           >  1:5782643/1‑93 (MQ=255)
                                                                                 gTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGc         >  1:1321873/1‑93 (MQ=255)
                                                                                  tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCa        >  1:2375816/1‑93 (MQ=255)
                                                                                  tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCa        >  1:3382356/1‑93 (MQ=255)
                                                                                   aTTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAg       >  1:2584122/1‑93 (MQ=255)
                                                                                   aTTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAg       >  1:2863223/1‑93 (MQ=255)
                                                                                   aTTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAg       <  1:1787417/93‑1 (MQ=255)
                                                                                   aTTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAg       <  2:1563254/93‑1 (MQ=255)
                                                                                     tctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTc     >  1:4488491/1‑93 (MQ=255)
                                                                                     tctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTc     <  1:5603334/93‑1 (MQ=255)
                                                                                      ctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTcc    <  2:3020921/93‑1 (MQ=255)
                                                                                      ctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTcc    <  1:4999144/93‑1 (MQ=255)
                                                                                       tcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCt   >  2:1830365/1‑93 (MQ=255)
                                                                                              |                                                                                       
GCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTG  >  NC_000913/4013679‑4013857

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