Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,049 C→G 16.7% G90A (GGC→GCC)  yijO ← AraC family putative transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,0490CG16.7% 471.3 / 88.4 181G90A (GGC→GCC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (77/72);  new base G (13/17);  total (91/90)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 187 total aligned reads displayed.

GAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTC  >  NC_000913/4147950‑4148148
                                                                                                   |                                                                                                   
 aaCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGcc                                                                                                    <  1:2626007/100‑1 (MQ=255)
 aaCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGcc                                                                                                    <  1:2030329/100‑1 (MQ=255)
  aCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCa                                                                                                   >  1:744622/1‑100 (MQ=255)
   ccTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCCGAAGCCGAGTTGCTGTTTGCCaa                                                                                                  <  1:195338/100‑1 (MQ=255)
    cTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCaaa                                                                                                 <  2:374121/100‑1 (MQ=255)
     tATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAAc                                                                                                >  1:2961447/1‑100 (MQ=255)
       tGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAg                                                                                              >  2:237971/1‑100 (MQ=255)
        gCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGc                                                                                             <  1:2518998/100‑1 (MQ=255)
        gCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGc                                                                                             <  1:237971/100‑1 (MQ=255)
         cGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGc                                                                                             >  2:2467419/1‑99 (MQ=255)
          gTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCAc                                                                                           >  2:2023824/1‑100 (MQ=255)
           tGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaac                                                                                          <  2:242608/100‑3 (MQ=255)
           tGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaac                                                                                          <  1:1439499/100‑3 (MQ=255)
           tGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaac                                                                                          <  1:3296956/100‑3 (MQ=255)
             ggCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactc                                                                                        <  1:3589864/100‑5 (MQ=255)
             ggCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCt                                                                                        >  1:312970/1‑100 (MQ=255)
              gCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactcg                                                                                       <  1:148716/100‑6 (MQ=255)
               ccGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactcgg                                                                                      <  1:2537991/100‑7 (MQ=255)
               ccGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactcgg                                                                                      <  1:1275805/100‑7 (MQ=255)
               ccGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactcg                                                                                       <  2:3261067/99‑6 (MQ=255)
                 gcgTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                    <  1:409097/100‑1 (MQ=255)
                  cgTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAgg                                                                                   <  1:1155283/100‑1 (MQ=255)
                   gTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgt                                                                                  <  1:412182/100‑1 (MQ=255)
                     cGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtag                                                                                >  1:402996/1‑100 (MQ=255)
                      gggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagt                                                                               >  1:981849/1‑100 (MQ=255)
                      gggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagt                                                                               <  2:3040589/100‑1 (MQ=255)
                      gggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagt                                                                               <  1:2422194/100‑1 (MQ=255)
                           aCGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgg                                                                          >  1:102819/1‑100 (MQ=255)
                              tgtTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggggg                                                                       >  2:3378386/1‑100 (MQ=255)
                                ttGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCt                                                                     <  1:2932058/100‑1 (MQ=255)
                                ttGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCt                                                                     <  1:212832/100‑1 (MQ=255)
                                 tGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCtt                                                                    >  2:268112/1‑100 (MQ=255)
                                 tGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCtt                                                                    >  2:479178/1‑100 (MQ=255)
                                  gCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTg                                                                   <  1:1837249/100‑1 (MQ=255)
                                  gCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTg                                                                   >  1:2064980/1‑100 (MQ=255)
                                   cTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGc                                                                  <  2:2095366/100‑1 (MQ=255)
                                     tCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCa                                                                >  2:2947845/1‑100 (MQ=255)
                                     tCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCa                                                                <  2:561602/100‑1 (MQ=255)
                                       gCCAGATTTTGCTATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAt                                                               >  2:1391974/1‑99 (MQ=255)
                                       gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTGGCGGGGGCTTGCCAt                                                               >  2:784748/1‑99 (MQ=255)
                                       gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                              >  1:764186/1‑100 (MQ=255)
                                       gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                              >  1:509470/1‑100 (MQ=255)
                                       gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGATAAAGGTAGTAGCGGGGGCTTGCCAtt                                                              >  1:698005/1‑100 (MQ=255)
                                        ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtg                                                             <  1:3162622/100‑1 (MQ=255)
                                        ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAACGGGGGCTTGCCATtg                                                             >  2:2123938/1‑100 (MQ=255)
                                          aGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgtg                                                           <  2:3532609/100‑1 (MQ=255)
                                             ttttGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaa                                                        >  1:371890/1‑100 (MQ=255)
                                                tGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAAt                                                     <  1:2646207/100‑1 (MQ=255)
                                                 gATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAAtt                                                    >  1:1611333/1‑100 (MQ=255)
                                                  atatTGGTTGCCATCCCATCGAACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTCGTAGCGGGGGCTTGACATTGTGGaa                                                        >  2:229983/1‑95 (MQ=255)
                                                   tatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTcc                                                  >  1:1549511/1‑100 (MQ=255)
                                                    atTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCa                                                 <  1:2284433/100‑1 (MQ=255)
                                                     ttGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAg                                                <  2:2540043/100‑1 (MQ=255)
                                                     ttGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAg                                                <  1:3616626/100‑1 (MQ=255)
                                                     ttGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAg                                                >  2:844870/1‑100 (MQ=255)
                                                       ntttGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGcc                                              >  1:101685/2‑100 (MQ=255)
                                                       gTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGcc                                              >  1:2759108/1‑100 (MQ=255)
                                                       gTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGcc                                              >  1:2184239/1‑100 (MQ=255)
                                                         ttGCCCTTCCATTACACGACGCTGAAGCAGAGTTGCTGTTTGACAAATAGAACGCTAAAGGCAGTAGCGAGGGCTTGCCATTGTGGAAAATTCCAGc                                               <  2:1741934/97‑1 (MQ=255)
                                                          tGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccacc                                           <  1:2114270/100‑1 (MQ=255)
                                                          tGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccacc                                           >  1:348175/1‑100 (MQ=255)
                                                          tGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccacc                                           <  2:1994701/100‑1 (MQ=255)
                                                           gCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccacca                                          >  1:1620080/1‑100 (MQ=255)
                                                            ccATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcag                                         >  1:1950827/1‑100 (MQ=255)
                                                             catcgcgctGCACGACACTGAAGCAGAGTCGCCGTTTGCCAAACAGAACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcag                                         <  2:440170/91‑1 (MQ=255)
                                                             cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagc                                        <  1:3090498/100‑1 (MQ=255)
                                                              aTCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagca                                       <  1:2944338/100‑1 (MQ=255)
                                                                cccATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAgg                                     <  2:2766826/100‑1 (MQ=255)
                                                                  cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAACCACCAGCAGGaa                                   >  2:2784598/1‑100 (MQ=255)
                                                                  cATTGCACGACGCTGAAGCCGAGTTGATGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCGTGCCATTGTGGAAAATTCCAGCCACCAGCAGGaa                                   >  2:1368937/1‑100 (MQ=255)
                                                                   aTTGCACGACGCTGAAGCCGAGTTGCTGTTTGACAAACAGCACGCTAAATGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAAc                                  >  2:3426925/1‑100 (MQ=255)
                                                                      gCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTa                               >  1:59421/1‑100 (MQ=255)
                                                                       cACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTAc                              <  2:1834546/100‑1 (MQ=255)
                                                                         cgacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAg                            >  1:1105316/1‑100 (MQ=255)
                                                                         cgacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAg                            >  1:135598/1‑100 (MQ=255)
                                                                          gacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGc                           <  1:2262896/100‑1 (MQ=255)
                                                                          gacgCTGAAGCCGAGTTGCTGTGTGCCAAACAGCACGCTAAAGGTAGTACCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCATGAACGTACAGc                           <  2:2826157/100‑1 (MQ=255)
                                                                                gAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCg                     >  2:201101/1‑100 (MQ=255)
                                                                                 aaGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGc                    >  2:1974937/1‑100 (MQ=255)
                                                                                    ccGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCgg                 <  1:3027545/100‑1 (MQ=255)
                                                                                     cGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCggg                >  2:442377/1‑100 (MQ=255)
                                                                                      gAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCGCCAGCAGGAACGTACAGCACATCGCCgg                 <  1:2410192/98‑1 (MQ=255)
                                                                                      gAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCggg                <  1:1578100/99‑1 (MQ=255)
                                                                                       gcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGaa              >  1:1888115/4‑100 (MQ=255)
                                                                                        gTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc             <  2:981842/100‑1 (MQ=255)
                                                                                        cgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCAGATCggaagagcg             >  2:465796/3‑92 (MQ=255)
                                                                                        cgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc             >  2:177355/3‑100 (MQ=255)
                                                                                        cgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc             >  1:2046168/3‑100 (MQ=255)
                                                                                         gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt            >  1:2917877/2‑100 (MQ=255)
                                                                                         gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt            >  1:289270/2‑100 (MQ=255)
                                                                                            cTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAAc         >  1:1726050/1‑100 (MQ=255)
                                                                                            cTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAAc         >  1:866959/1‑100 (MQ=255)
                                                                                             tGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCATCAGCAGGAACGTACAGCACATCGCCGGGAACTAACg        <  1:2564567/100‑1 (MQ=255)
                                                                                             tGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACg        >  2:1828970/1‑100 (MQ=255)
                                                                                              gTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGt       >  1:222264/1‑100 (MQ=255)
                                                                                               tttGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGtt      <  2:1531978/100‑1 (MQ=255)
                                                                                                 tGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc    >  2:2073126/1‑100 (MQ=255)
                                                                                                  gCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCt   <  2:3202695/100‑1 (MQ=255)
                                                                                                   ccAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTc  >  2:274439/1‑100 (MQ=255)
                                                                                                   ccAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTc  <  2:1150846/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTC  >  NC_000913/4147950‑4148148

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