Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,418,264 G→C 5.4% P64A (CCG→GCG) ‡ ulaR ← transcriptional repressor for the L‑ascorbate utilization divergent operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,418,2640GC 5.4% 321.3 / 8.1 92P64A (CCG→GCG) ‡ulaRtranscriptional repressor for the L‑ascorbate utilization divergent operon
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base C (3/2);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CACATTTCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGC  >  NC_000913/4418116‑4418411
                                                                                                                                                    |                                                                                                                                                   
cacaTTTCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCgg                                                                                                                                                     <  1:159424/149‑1 (MQ=255)
  caTTTCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCggg                                                                                                                                                    >  4:444814/1‑148 (MQ=255)
   aTTTCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGgcg                                                                                                                                                  >  1:579732/1‑149 (MQ=255)
    nttCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGgcg                                                                                                                                                  >  3:108134/2‑148 (MQ=255)
     ttCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGgcgc                                                                                                                                                 >  4:470303/1‑148 (MQ=255)
        ccGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGtt                                                                                                                                             <  2:581816/149‑1 (MQ=255)
        ccGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGt                                                                                                                                              <  4:678648/148‑1 (MQ=255)
        ccGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCCCGCGCTgga                                                                                                                                             <  1:655684/149‑3 (MQ=255)
        ccGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCCCGCGCTgg                                                                                                                                              <  3:178851/148‑2 (MQ=255)
          gCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTgg                                                                                                                                           <  1:39816/149‑1 (MQ=255)
          gCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTgg                                                                                                                                           <  1:372314/149‑1 (MQ=255)
          gCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTgg                                                                                                                                           <  2:579732/149‑1 (MQ=255)
          gCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTg                                                                                                                                            <  3:470303/148‑1 (MQ=255)
          gCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTg                                                                                                                                            <  4:482160/148‑1 (MQ=255)
            ccAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGt                                                                                                                                         <  2:188427/149‑1 (MQ=255)
             cAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCCAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCGGTTGGGt                                                                                                                                         >  3:485975/1‑148 (MQ=255)
                acAGGAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTCATa                                                                                                                                     >  2:66667/2‑149 (MQ=255)
                 cAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATa                                                                                                                                     >  4:867199/1‑148 (MQ=255)
                  aGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAg                                                                                                                                    <  4:874458/148‑1 (MQ=255)
                         cgGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCtt                                                                                                                                 >  2:162118/1‑144 (MQ=255)
                         cgGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCtt                                                                                                                                 <  1:162118/144‑1 (MQ=255)
                              gAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCAtt                                                                                                                       <  2:534277/149‑1 (MQ=255)
                               aGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCAtt                                                                                                                       <  4:657643/148‑1 (MQ=255)
                               aGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTg                                                                                                                      <  1:89736/149‑1 (MQ=255)
                                    cAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCa                                                                                                                  >  3:406245/1‑148 (MQ=255)
                                    cAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCa                                                                                                                  <  4:1055291/148‑1 (MQ=255)
                                      gTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCAATGGCGCaga                                                                                                                <  4:673415/148‑3 (MQ=255)
                                          atgaCTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACttttt                                                                                                            <  3:444814/148‑1 (MQ=255)
                                          atgaCTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACttttt                                                                                                            >  4:25448/1‑148 (MQ=255)
                                            gaCTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCa                                                                                                         <  1:358159/149‑1 (MQ=255)
                                              cTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGt                                                                                                       >  2:206255/1‑149 (MQ=255)
                                                aCGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAg                                                                                                        >  4:883876/1‑146 (MQ=255)
                                                aCGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAg                                                                                                        <  3:883876/146‑1 (MQ=255)
                                                aCGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGttt                                                                                                     >  1:777686/1‑149 (MQ=255)
                                                    tttCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGt                                                                                                       >  2:416899/1‑143 (MQ=255)
                                                    tttCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGt                                                                                                       <  1:416899/143‑1 (MQ=255)
                                                         cccGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTc                                                                                            <  2:494041/149‑1 (MQ=255)
                                                         cccGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTGTCGTCGTGATTCTGTGCCTTTTGCAGAGTCATCGGCGGCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTATTTCAGTTTGCCGCttg                                                                                             <  4:110443/148‑2 (MQ=255)
                                                              aTTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCa                                                                                        >  4:353609/1‑148 (MQ=255)
                                                              aTTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCa                                                                                        <  4:92135/148‑1 (MQ=255)
                                                               ttAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCaa                                                                                       <  3:267572/148‑1 (MQ=255)
                                                                 aaCCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGt                                                                                     >  3:867249/1‑148 (MQ=255)
                                                                  aCCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGtt                                                                                    >  4:546829/1‑148 (MQ=255)
                                                                  aCCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATACCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGtt                                                                                    >  3:713554/1‑148 (MQ=255)
                                                                   ccAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGttt                                                                                   >  3:822387/1‑148 (MQ=255)
                                                                       cTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTg                                                                               >  3:442134/1‑148 (MQ=255)
                                                                             ggCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATAtcgc                                                                        <  2:118885/149‑1 (MQ=255)
                                                                             ggCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATAtcgc                                                                        <  2:826853/149‑1 (MQ=255)
                                                                                cGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATAtctcgt                                                                      >  4:913551/1‑148 (MQ=255)
                                                                                cGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATAtcgcgt                                                                      <  4:108134/148‑1 (MQ=255)
                                                                                        cGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGt                                                              >  4:417706/1‑148 (MQ=255)
                                                                                        cGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTg                                                             >  1:568568/1‑149 (MQ=255)
                                                                                          aTACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTgg                                                            >  4:428673/1‑148 (MQ=255)
                                                                                           tACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGGACTTTTTTCAGTTTGCCGCTTGCGTCAAGTTTATTGATATCGCGTCGCGCAGTAGc                                                           >  4:652866/1‑148 (MQ=255)
                                                                                            aCGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCa                                                                                                                  >  4:1005210/1‑92 (MQ=255)
                                                                                            aCGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCa                                                                                                                  <  3:1005210/92‑1 (MQ=255)
                                                                                            aCGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTACGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAg                                                         <  1:336713/149‑1 (MQ=255)
                                                                                              gTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGc                                                       >  2:713732/1‑149 (MQ=255)
                                                                                               tACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGc                                                       >  3:142083/1‑148 (MQ=255)
                                                                                               tACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGc                                                       >  3:762257/1‑148 (MQ=255)
                                                                                                aCTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGa                                                     >  1:496484/1‑149 (MQ=255)
                                                                                                  tttttCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGaaa                                                   >  1:682345/1‑149 (MQ=255)
                                                                                                   ttttCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAAt                                                  >  2:10045/1‑149 (MQ=255)
                                                                                                        aTCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTccc                                              <  3:353609/148‑1 (MQ=255)
                                                                                                        aTCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTccc                                              <  4:142083/148‑1 (MQ=255)
                                                                                                         tCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTa                                                                                  >  1:539441/1‑111 (MQ=255)
                                                                                                         tCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTa                                                                                  <  2:539441/111‑1 (MQ=255)
                                                                                                         tCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAg                                            >  1:485349/1‑149 (MQ=255)
                                                                                                          cGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAg                                            >  3:849413/1‑148 (MQ=255)
                                                                                                           gTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGAc                                          >  2:464800/1‑149 (MQ=255)
                                                                                                               ttCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCt                                       <  4:510432/148‑1 (MQ=255)
                                                                                                               ttCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTc                                      >  1:568951/1‑149 (MQ=255)
                                                                                                                  tGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAAc                                   <  1:206255/149‑1 (MQ=255)
                                                                                                                     gcCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGa                                 <  3:417706/148‑1 (MQ=255)
                                                                                                                     gcCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGAc                                >  2:683877/1‑149 (MQ=255)
                                                                                                                      ccTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGAc                                >  3:245737/1‑148 (MQ=255)
                                                                                                                      ccTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACt                               <  2:747814/149‑1 (MQ=255)
                                                                                                                        tGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACttt                             >  1:466038/1‑149 (MQ=255)
                                                                                                                            gCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCa                         >  2:1020746/1‑149 (MQ=255)
                                                                                                                             cAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCaa                        <  1:683877/149‑1 (MQ=255)
                                                                                                                             cAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCa                         <  4:245737/148‑1 (MQ=255)
                                                                                                                               gATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACg                      <  2:799318/149‑1 (MQ=255)
                                                                                                                                aTTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACgg                     >  2:435868/1‑149 (MQ=255)
                                                                                                                                  tCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTc                   <  1:1020746/149‑1 (MQ=255)
                                                                                                                                    aTCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTcac                 >  2:586163/1‑149 (MQ=255)
                                                                                                                                    aTCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTcac                 >  2:44885/1‑149 (MQ=255)
                                                                                                                                    aTCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTcac                 <  1:694639/149‑1 (MQ=255)
                                                                                                                                         ccaacAGCGCCCGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGc            >  1:1069224/5‑149 (MQ=255)
                                                                                                                                          caacAGCGCCCGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGcc           >  1:465834/4‑149 (MQ=255)
                                                                                                                                           tCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGcc           >  4:944745/1‑148 (MQ=255)
                                                                                                                                           aacAGCGCCCGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGcc           >  3:558948/3‑148 (MQ=255)
                                                                                                                                              aGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAAt       >  2:922718/1‑149 (MQ=255)
                                                                                                                                                 gcgGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATtg     <  3:867199/148‑1 (MQ=255)
                                                                                                                                                 gcgGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATtg     <  3:363802/148‑1 (MQ=255)
                                                                                                                                                   gggCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATtgtg   <  4:713554/148‑1 (MQ=255)
                                                                                                                                                   gggCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATtgtg   <  4:849413/148‑1 (MQ=255)
                                                                                                                                                   gggCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGc  >  1:614797/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACATTTCCCGCCCAAGCAGAAACGCGGTGGAGCCGCAGTTGATGACTACGCTTTCGCCCGGATTAACCAGCTGCGAGGCCGCTTTAGCGATACGTACTTTTTCATCGTGATTCTGCGCCTGATGCAGATTCATCGGCGTCCAGCGCGGGCGCTGTTGGGTAATAGCTTCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGC  >  NC_000913/4418116‑4418411

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