Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,776 T→A 18.1% intergenic (+24/‑15) crp → / → yhfK cAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,7760TA18.1% 468.9 / 97.8 205intergenic (+24/‑15)crp/yhfKcAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (80/86);  new base A (15/22);  total (95/110)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-02

Only 100 of 212 total aligned reads displayed.

CATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCG  >  NC_000913/3486677‑3486875
                                                                                                   |                                                                                                   
cATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                     <  1:3841452/100‑1 (MQ=255)
cATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                     >  1:2852439/1‑100 (MQ=255)
cATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                     >  2:3669825/1‑100 (MQ=255)
 aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                                                                                    <  1:1252683/100‑1 (MQ=255)
 aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                                                                                    >  2:1082504/1‑100 (MQ=255)
   tCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGt                                                                                                  <  2:14301/100‑1 (MQ=255)
     tGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAg                                                                                                <  1:2999817/100‑1 (MQ=255)
     tGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAg                                                                                                >  2:2999833/1‑100 (MQ=255)
        aGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgc                                                                                             <  2:3211834/100‑1 (MQ=255)
          aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                           >  2:1710535/1‑100 (MQ=255)
          aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                           >  2:658239/1‑100 (MQ=255)
          aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcgc                                                                                           <  2:3700442/100‑1 (MQ=255)
          aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcgc                                                                                           <  2:246516/100‑1 (MQ=255)
          aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAATGGCGCGCTACCAGGTAacgcgc                                                                                           <  2:3283270/100‑1 (MQ=255)
           tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAACACCATCGTCGTTTACGGCACTCGTTAAACCCGTCGGAGTGGCGCGTTACATGGTAgcccg                                                                                            >  1:3372059/1‑98 (MQ=255)
           tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                          <  2:3573666/100‑1 (MQ=255)
           tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                          <  2:3122739/100‑1 (MQ=255)
           tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                          <  2:2595258/100‑1 (MQ=255)
            gCTGGAAGATCAGCACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACGTGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                          <  1:830258/99‑1 (MQ=255)
            gCTGGAAGATCAGAACCTGGTCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCa                                                                                         <  2:1164078/100‑1 (MQ=255)
            gCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCa                                                                                         <  2:1998979/100‑1 (MQ=255)
             cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCac                                                                                        <  2:1896184/100‑2 (MQ=39)
             cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCac                                                                                        <  2:1769362/100‑2 (MQ=39)
             cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCac                                                                                        <  2:1320352/100‑2 (MQ=39)
               ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                                                                                      >  2:2459329/1‑100 (MQ=255)
               ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                                                                                      >  1:42720/1‑100 (MQ=255)
               ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                                                                                      >  1:886871/1‑100 (MQ=255)
                  aGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGt                                                                                   >  2:460885/1‑100 (MQ=255)
                  aGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGt                                                                                   <  1:2811152/100‑1 (MQ=255)
                   gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                                  <  2:2509414/100‑1 (MQ=255)
                   gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                                  <  1:1739750/100‑1 (MQ=255)
                   gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTAGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                                  >  2:685718/1‑100 (MQ=255)
                    aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGttt                                                                                 >  2:3516540/1‑100 (MQ=255)
                    aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGttt                                                                                 >  2:4076450/1‑100 (MQ=255)
                    aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGttt                                                                                 >  2:3043917/1‑100 (MQ=255)
                    aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                                  <  2:3974243/99‑1 (MQ=255)
                       aGAACGGGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCACGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTctc                                                                              >  1:4025459/1‑98 (MQ=255)
                       aGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTccc                                                                              <  1:2157819/100‑1 (MQ=255)
                       aGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTcc                                                                               <  2:271096/99‑1 (MQ=255)
                             tGATCTCCGCCCACGGTAAAACCATCGTCGTTTCCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTcc                                                                               >  1:3264739/1‑93 (MQ=255)
                             tGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAtg                                                                        <  2:2178128/100‑1 (MQ=255)
                              gATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTGTGTTTCCCCCGAtgt                                                                       >  2:4340789/1‑100 (MQ=255)
                               aTCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAtgtg                                                                      >  2:902732/1‑100 (MQ=255)
                                  tacGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGgcg                                                                   <  2:2397074/98‑1 (MQ=255)
                                        cacGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACt                                                             >  1:3385766/1‑100 (MQ=255)
                                           ggTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                          <  2:55086/100‑1 (MQ=255)
                                           ggTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                          <  2:3332855/100‑1 (MQ=255)
                                           ggTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                          <  1:2568237/100‑1 (MQ=255)
                                              aaaaCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTa                                                       >  1:1241096/1‑100 (MQ=255)
                                                  ccATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGCTTCCCCCGATGTGGCGCAGACTGATTTATCAc                                                   >  1:235845/1‑100 (MQ=255)
                                                  ccATCGACGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTCGTCGCGCTCCCTTTCGTCTCCCCCGATGTCGCGCAGACTGATTTAt                                                      >  1:2743618/1‑97 (MQ=255)
                                                      cgtcgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCg                                               <  1:797431/100‑1 (MQ=255)
                                                         cgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGata                                            <  2:3771058/100‑1 (MQ=255)
                                                            ttACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTCGCGTGCCATTTTGTTTCACACGATGTGGCGAAGACTGATTTATCACCCCGATATCa                                         >  1:3344721/1‑100 (MQ=255)
                                                             tACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCTCCCGATGTTGCGCAGACTGATTTATCACCCCGatatga                                         >  1:56823/1‑97 (MQ=255)
                                                              aCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAc                                       <  1:1668329/100‑1 (MQ=255)
                                                               cGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACt                                      <  2:4239612/100‑1 (MQ=255)
                                                                ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCaa                                        <  1:3200788/97‑1 (MQ=255)
                                                                ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTa                                     <  1:115884/100‑1 (MQ=255)
                                                                ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTa                                     >  2:860777/1‑100 (MQ=255)
                                                                 gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                                    >  1:1487870/1‑100 (MQ=255)
                                                                 gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCACCGATGTGGCGCAGACTGATTTATCATCCCGATATCAACTAt                                    >  1:3083233/1‑100 (MQ=255)
                                                                 gCAATCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                                    <  1:1976133/100‑1 (MQ=255)
                                                                  cACTCGTTCCTCCCGTCGGAGTGGCGCGTTACCTGATAGCGCGTCATTTTGTTTCCCCCGATGTGTCGCAGACTGATTTATCACCCCGATATCAACTa                                     <  1:1040146/98‑1 (MQ=255)
                                                                    cTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCa                                 >  2:2519096/1‑100 (MQ=255)
                                                                     tCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCAc                                >  2:1541163/1‑100 (MQ=255)
                                                                     tCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCAc                                <  2:247672/100‑1 (MQ=255)
                                                                      cGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACt                               >  2:3444364/1‑100 (MQ=255)
                                                                         tAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCg                            >  2:829769/1‑100 (MQ=255)
                                                                         tAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCg                            <  1:3940847/100‑1 (MQ=255)
                                                                          aaTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGa                           <  1:2188374/100‑1 (MQ=255)
                                                                          aaTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGa                           <  2:2362241/100‑1 (MQ=255)
                                                                           aTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGAc                          >  2:359643/1‑100 (MQ=255)
                                                                              gaaacaaaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTTGACtaa                       >  1:3562564/10‑100 (MQ=25)
                                                                              ccGTCTGAGTGGCGCGGTACCTGCGAGCACGCCAGTGTGGTACGCCCGGTGTGGCGCAGACTGATTAATGACCCCGATATCAACTATGCACTTCGACaa                        <  1:3730814/99‑1 (MQ=255)
                                                                               cGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAAc                      <  1:3669803/100‑1 (MQ=255)
                                                                                gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACg                     >  2:3562499/1‑100 (MQ=255)
                                                                                  caaaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCt                   >  2:2791293/6‑100 (MQ=38)
                                                                                    aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACTCCGATATCAACTATGCACTTCGACAAACGCTg                  >  1:1881776/4‑99 (MQ=38)
                                                                                    aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg                 >  2:913514/4‑100 (MQ=38)
                                                                                    aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg                 >  2:298131/4‑100 (MQ=38)
                                                                                    aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg                 >  1:256115/4‑100 (MQ=38)
                                                                                    aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg                  >  2:1512057/4‑99 (MQ=38)
                                                                                      gTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTACCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTg               <  1:527741/100‑1 (MQ=255)
                                                                                       tGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATCTCAACTATGCACTTCGACAAACGCTGGTGc              >  1:4201290/1‑100 (MQ=255)
                                                                                       tGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGc              <  2:1051173/100‑1 (MQ=255)
                                                                                       tGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGc              >  2:4162378/1‑100 (MQ=255)
                                                                                       tGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGc              <  1:1780034/100‑1 (MQ=255)
                                                                                        ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGca             >  1:1051169/1‑99 (MQ=255)
                                                                                        ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt             >  2:119633/1‑100 (MQ=255)
                                                                                        ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt             >  2:538052/1‑100 (MQ=255)
                                                                                         gcgcGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTa            <  1:880238/100‑1 (MQ=255)
                                                                                           gcgTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATg          <  1:3757320/100‑1 (MQ=255)
                                                                                            cgctACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGt         >  2:4037195/4‑100 (MQ=255)
                                                                                               tACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTg      >  2:2726248/1‑100 (MQ=255)
                                                                                                aCCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGc     <  1:261697/100‑1 (MQ=255)
                                                                                                 ccTGGTAGCGGGCCATTTTGTTTCCACCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGcc    <  2:42720/100‑1 (MQ=255)
                                                                                                  cTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGccc   >  1:676772/1‑100 (MQ=255)
                                                                                                   tGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCg  <  1:2488545/100‑1 (MQ=255)
                                                                                                   tGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCg  >  2:1339381/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCG  >  NC_000913/3486677‑3486875

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: