Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,684,604 G→A 100% Q339* (CAG→TAG)  yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,6040GA97.6% 547.0 / 1.4 166Q339* (CAG→TAG) yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (75/86);  total (80/86)
Fisher's exact test for biased strand distribution p-value = 5.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 166 total aligned reads displayed.

GATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACA  >  NC_000913/3684508‑3684699
                                                                                                |                                                                                               
gaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAcc                                                                                              <  2:2809755/100‑1 (MQ=255)
gaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAcc                                                                                              >  2:1356348/1‑100 (MQ=255)
  tCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTa                                                                                            <  2:2437559/100‑1 (MQ=255)
   cGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTaa                                                                                           <  2:23154/100‑1 (MQ=255)
    gCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAg                                                                                          <  2:1672397/100‑1 (MQ=255)
      cGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAgg                                                                                         >  2:3443587/1‑99 (MQ=255)
       gTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTg                                                                                       <  1:1141909/100‑1 (MQ=255)
       gTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTg                                                                                       <  2:325436/100‑1 (MQ=255)
           aaCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTg                                                                                   <  2:1700655/100‑1 (MQ=255)
            aCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGc                                                                                  <  2:1559080/100‑1 (MQ=255)
                   cacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAc                                                                           <  2:3589930/100‑1 (MQ=255)
                   cacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAc                                                                           <  2:484997/100‑1 (MQ=255)
                    acgccaATGCTGCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAc                                                                           <  1:997303/99‑1 (MQ=255)
                    acgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACg                                                                          <  1:2494233/100‑1 (MQ=255)
                    acgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACg                                                                          <  2:1312965/100‑1 (MQ=255)
                       ccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGtt                                                                       >  1:1639905/1‑100 (MQ=255)
                       ccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGtt                                                                       >  2:669229/1‑100 (MQ=255)
                        caATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTc                                                                      <  2:383484/100‑1 (MQ=255)
                         aaTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTcc                                                                     <  2:3190755/100‑1 (MQ=255)
                         aaTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTcc                                                                     <  2:513803/100‑1 (MQ=255)
                         aaTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTcc                                                                     <  1:420554/100‑1 (MQ=255)
                         aaTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTcc                                                                     <  1:126750/100‑1 (MQ=255)
                           tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  2:1092103/100‑1 (MQ=255)
                           tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  2:1818752/100‑1 (MQ=255)
                           tGCTACAGTGCCGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  2:3689900/100‑1 (MQ=255)
                            gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  <  1:3806167/100‑1 (MQ=255)
                            gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  >  2:1521149/1‑100 (MQ=255)
                            gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  >  1:3095411/1‑100 (MQ=255)
                            gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  <  2:2925802/100‑1 (MQ=255)
                             cTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTGGCTCAGGATAGTGAGCACATGCTAACCTAAGGTGATTGCGTGCCACGGTTCATGTa                                                                 >  2:2783148/1‑100 (MQ=255)
                             cTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTa                                                                 >  2:1821019/1‑100 (MQ=255)
                              tACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTac                                                                >  1:2692773/1‑100 (MQ=255)
                              tACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTac                                                                <  1:1435909/100‑1 (MQ=255)
                                  gTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTCACCTAAAGTGATTGCGTGCCACGGGTCCTGCACACCg                                                            >  1:3091757/1‑100 (MQ=255)
                                  gTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCg                                                            >  2:1506101/1‑100 (MQ=255)
                                  gTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCg                                                            >  1:2702249/1‑100 (MQ=255)
                                    gCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTt                                                          >  1:1940375‑M2/1‑99 (MQ=255)
                                     cggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTtg                                                         >  2:1407817‑M2/1‑98 (MQ=255)
                                     cggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTtg                                                         <  2:2080328‑M2/100‑3 (MQ=255)
                                      ggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTtgg                                                        >  2:2140319‑M2/1‑97 (MQ=255)
                                      ggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTtgg                                                        <  2:1195519‑M2/100‑4 (MQ=255)
                                       gacggaGTTGAATATGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTtggc                                                       >  1:2776869‑M2/1‑96 (MQ=255)
                                        acggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAacggtgtacaggaaccgtggca                                                      <  2:3846857‑M1/100‑23 (MQ=255)
                                        acggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAacggtgtacaggaaccgtggca                                                      <  2:2266082‑M1/100‑23 (MQ=255)
                                         cggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAacggtgtacaggaaccgtggcac                                                     <  1:722291‑M1/100‑24 (MQ=255)
                                          ggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAAt                                                    <  2:939537/100‑1 (MQ=255)
                                          ggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAAt                                                    <  2:3250813/100‑1 (MQ=255)
                                          ggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAAt                                                    <  2:220605/100‑1 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   <  2:3705966/100‑1 (MQ=255)
                                           gaGTTGAATACGTTCAATCGGCAGGCGCGCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   <  1:3184292/100‑1 (MQ=255)
                                            aGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGa                                                  >  1:868885/1‑100 (MQ=255)
                                             gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCAacggtgtacaggaaccgtggcacgcaa                                                 <  1:1870522‑M1/100‑28 (MQ=255)
                                             gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                                 <  2:1652691/100‑1 (MQ=255)
                                              ttGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACa                                                <  1:3786313/100‑1 (MQ=255)
                                              ttGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGATGATTGCGTGCCAacggtgtacaggaaccgtggcacgcaat                                                <  1:1035289‑M1/100‑29 (MQ=255)
                                               tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               <  1:1336297/100‑1 (MQ=255)
                                               tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               >  1:2342520/1‑100 (MQ=255)
                                                gAATACGTTCAATCGGCAGGCGCTCGCTCATGGTAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                              >  1:1918024/1‑100 (MQ=255)
                                                gAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                              >  2:392853/1‑100 (MQ=255)
                                                  aTACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCaa                                            >  2:3659962/1‑100 (MQ=255)
                                                      gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                        >  1:2762495/1‑100 (MQ=255)
                                                      gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                        >  1:2584426/1‑100 (MQ=255)
                                                       ttCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCa                                       <  1:2444343/100‑1 (MQ=255)
                                                        tCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCAt                                      >  1:1830588/1‑100 (MQ=255)
                                                          aaTCGGTAGGGTCGCGTTGATGATAGTGTGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTTCTGTCGACCGTTGTCAATGACAGCAAAGGCAt                                      >  1:926259/1‑98 (MQ=255)
                                                           aTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                                   >  1:1960835/1‑100 (MQ=255)
                                                            tCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAcc                                  <  2:3799992/100‑1 (MQ=255)
                                                            tCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAcc                                  <  2:2712232/100‑1 (MQ=255)
                                                              ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                <  1:2918664/100‑1 (MQ=255)
                                                              ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                >  2:3028219/1‑100 (MQ=255)
                                                               gCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCt                               <  2:732282/100‑1 (MQ=255)
                                                                    cgctcgctCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCt                          >  1:3628039/1‑100 (MQ=255)
                                                                     gctcgctgATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTg                         <  1:326210/100‑1 (MQ=255)
                                                                       tcgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcg                       <  1:1521140/100‑1 (MQ=255)
                                                                        cgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcgc                      <  2:3334633/100‑1 (MQ=255)
                                                                        cgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcg                       >  1:3300523/1‑99 (MQ=255)
                                                                           tcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg                   >  2:700791/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  >  1:2882925/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  >  2:921730/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  >  2:875675/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  >  2:2565938/1‑100 (MQ=255)
                                                                            cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCAACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg                  <  1:2849016/100‑1 (MQ=255)
                                                                             aTGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACCGCAAAGACATAACCGCTAATCTGCGCGAGGa                 >  2:1773837/1‑100 (MQ=255)
                                                                             aTGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGa                 >  2:2104745/1‑100 (MQ=255)
                                                                              tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt                >  2:3087049/1‑100 (MQ=255)
                                                                               gATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATc               <  2:3368430/100‑1 (MQ=255)
                                                                                aTAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCa              <  2:1913740/100‑1 (MQ=255)
                                                                                   gTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATAc           >  1:3875457/1‑100 (MQ=255)
                                                                                   gTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATAc           >  1:3452716/1‑100 (MQ=255)
                                                                                   gTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATAc           >  1:353926/1‑100 (MQ=255)
                                                                                   gTGAGCACTTGCTAACCTAAGGTGAGTGCGTGCCGTGGTTCATGTACTCCGTTGGCAATGACAGCTAAGTCATAACCGCTAATCTGCGCGAGGATCATa            <  1:593545/99‑1 (MQ=255)
                                                                                    tGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg          >  2:1140807/1‑100 (MQ=255)
                                                                                      aGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTg        <  2:2584444/100‑1 (MQ=255)
                                                                                          cTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGa    >  2:1709049/1‑100 (MQ=255)
                                                                                          cTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGa    >  1:2734486/1‑100 (MQ=255)
                                                                                           ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  1:1608867/100‑1 (MQ=255)
                                                                                           ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  1:353327/100‑1 (MQ=255)
                                                                                           ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  2:2128022/100‑1 (MQ=255)
                                                                                            tGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACa  >  2:2020947/1‑100 (MQ=255)
                                                                                            tGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACa  >  1:2128007/1‑100 (MQ=255)
                                                                                                |                                                                                               
GATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACA  >  NC_000913/3684508‑3684699

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: