Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,103,769 A→G F403L (TTC→CTC)  cpxA ← sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,103,7690AG97.4% 406.6 / ‑0.2 116F403L (TTC→CTC) cpxAsensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (61/52);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 124 total aligned reads displayed.

CCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCAC  >  NC_000913/4103673‑4103865
                                                                                                |                                                                                                
cccAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTCCGATAGAgcgg                                                                                               >  1:1224639/1‑100 (MQ=255)
cccAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCAGTACGATAGAgcgg                                                                                               >  2:2314969/1‑100 (MQ=255)
  cAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAgcggac                                                                                             <  1:99157/100‑1 (MQ=255)
  cAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAgcggac                                                                                             >  2:2069905/1‑100 (MQ=255)
   aGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAgcggacg                                                                                            >  2:1681693/1‑100 (MQ=255)
    gCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAgcggacgg                                                                                           <  1:1690221/100‑1 (MQ=255)
        cGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGaaaa                                                                                       <  2:309913/100‑1 (MQ=255)
        cGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGaaaa                                                                                       >  2:3163843/1‑100 (MQ=255)
         gATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAAt                                                                                      <  1:3180636/100‑1 (MQ=255)
          aTGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATc                                                                                     <  1:3320048/100‑1 (MQ=255)
          aTGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATc                                                                                     <  2:138337/100‑1 (MQ=255)
           tgctgcTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCt                                                                                    >  1:1726297/1‑100 (MQ=255)
           tgctgcTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACAGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCt                                                                                    >  1:1413139/1‑100 (MQ=255)
             ctgctgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGt                                                                                  >  1:2232548/1‑100 (MQ=255)
             ctgctgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGt                                                                                  <  1:2069909/100‑1 (MQ=255)
              tgctgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGtt                                                                                 <  1:2314975/100‑1 (MQ=255)
               gctgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTc                                                                                >  1:129811/1‑100 (MQ=255)
               gctgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTc                                                                                >  1:41531/1‑100 (MQ=255)
                 tgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTcgc                                                                              <  2:73482/100‑1 (MQ=255)
                 tgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTcgc                                                                              >  1:2806326/1‑100 (MQ=255)
                 tgAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTcgc                                                                              >  2:2978056/1‑100 (MQ=255)
                  gAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTcgcg                                                                             >  1:713963/1‑100 (MQ=255)
                  gAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTcgcg                                                                             >  2:673747/1‑100 (MQ=255)
                   aaTGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGa                                                                            >  1:1665305/1‑100 (MQ=255)
                   aaTGGCGGTTTCAACAATCACCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGa                                                                            <  2:2872971/100‑1 (MQ=255)
                    aTGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGAt                                                                           >  1:1882645/1‑100 (MQ=255)
                     tGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATc                                                                          <  1:227032/100‑1 (MQ=255)
                      ggCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCt                                                                         >  1:2572819/1‑100 (MQ=255)
                       gCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCtt                                                                        >  1:831135/1‑100 (MQ=255)
                       gCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCtt                                                                        >  2:672115/1‑100 (MQ=255)
                         ggTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTcc                                                                      >  2:3284586/1‑100 (MQ=255)
                          gTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCg                                                                     >  2:2457890/1‑100 (MQ=255)
                            ttCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCggg                                                                   <  1:1766329/100‑1 (MQ=255)
                            ttCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCggg                                                                   <  2:3355567/100‑1 (MQ=255)
                             tCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGc                                                                  <  1:1711689/100‑1 (MQ=255)
                             tCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGc                                                                  <  2:1251143/100‑1 (MQ=255)
                              caacaaTCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCt                                                                 >  2:120867/1‑100 (MQ=255)
                              caacaaTCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCt                                                                 >  2:866940/1‑100 (MQ=255)
                              caacaaTCACCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCt                                                                 >  2:3406683/1‑100 (MQ=255)
                                acaaTCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTaa                                                               >  1:217293/1‑100 (MQ=255)
                                 caaTCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAAc                                                              >  1:3148482/1‑100 (MQ=255)
                                 atctCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAAc                                                              <  1:2493306/97‑1 (MQ=255)
                                   aTCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACg                                                             <  1:3029026/99‑1 (MQ=255)
                                     ctccAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCa                                                          <  2:1763365/98‑1 (MQ=255)
                                     cGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAAcagat                                                          >  2:2493298/1‑96 (MQ=255)
                                      gCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAg                                                         <  1:2085622/100‑1 (MQ=255)
                                      gCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAg                                                         <  1:171450/100‑1 (MQ=255)
                                      gCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAg                                                         >  2:1230245/1‑100 (MQ=255)
                                      gCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGGGCTGCATCGGGACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAg                                                         >  2:1988403/1‑100 (MQ=255)
                                       ccAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCaag                                                        >  1:1763369/1‑98 (MQ=255)
                                         aGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGAc                                                      <  2:1922724/100‑1 (MQ=255)
                                         aGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGAc                                                      >  1:2272303/1‑100 (MQ=255)
                                          gCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGCTCTTCCGGGCTAACGCCAGGAc                                                      >  2:2379542/1‑99 (MQ=255)
                                          gCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGAcc                                                     >  2:3169482/1‑100 (MQ=255)
                                           ccccAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCa                                                    <  1:866941/100‑1 (MQ=255)
                                           ccccAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCa                                                    <  1:1971656/100‑1 (MQ=255)
                                             ccAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCAtc                                                  >  2:115415/1‑100 (MQ=255)
                                               aaaCCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCAtcgt                                                <  2:756502/100‑1 (MQ=255)
                                                    tGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTcc                                           >  1:1589356/1‑100 (MQ=255)
                                                    tGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTcc                                           <  2:3084548/100‑1 (MQ=255)
                                                    tGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTcc                                           <  2:1589353/100‑1 (MQ=255)
                                                        ccGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGAACCATCGGCGTCCACCg                                       <  2:299669/100‑1 (MQ=255)
                                                         cGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGt                                      <  2:2817260/100‑1 (MQ=255)
                                                          gCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTa                                     >  1:785123/1‑100 (MQ=255)
                                                             aGATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATg                                  >  1:167270/1‑100 (MQ=255)
                                                              gATTCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATgg                                 <  1:424888/100‑1 (MQ=255)
                                                                ttCACGATCGCGTGCTTCATCGGTACGATAGAGCGGCCGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTg                               <  2:1882641/100‑1 (MQ=255)
                                                                ttCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTg                               >  1:1899281/1‑100 (MQ=255)
                                                                ttCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTg                               <  2:217292/100‑1 (MQ=255)
                                                                ttCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTGCGGGCGAACGCCAGGACCATCGTCGTCCACCGTAATGGTg                               >  2:3164356/1‑100 (MQ=255)
                                                                 tCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGa                              >  1:64757/1‑100 (MQ=255)
                                                                 tCACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGa                              >  1:569479/1‑100 (MQ=255)
                                                                  cACGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGAt                             <  1:2429140/100‑1 (MQ=255)
                                                                    cGATCGCGTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATAc                           <  2:2910460/100‑1 (MQ=255)
                                                                        cgcgTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTc                       <  2:3116465/100‑1 (MQ=255)
                                                                        cgcgTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTc                       <  2:3388982/100‑1 (MQ=255)
                                                                         gcgTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCt                      <  1:3113235/100‑1 (MQ=255)
                                                                         gcgTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCt                      <  1:317703/100‑1 (MQ=255)
                                                                           gTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCttt                    >  2:48747/1‑100 (MQ=255)
                                                                           gTGCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCttt                    <  1:1276157/100‑1 (MQ=255)
                                                                            tGCTTCATCGCTCCGCTAGAGCGCACGGACAATCTGTTCGCGATCTTCCGGGCTAACCCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCt                      >  2:1686947/1‑97 (MQ=255)
                                                                             gCTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTAt                  >  1:1341645/1‑100 (MQ=255)
                                                                              cTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATc                 <  2:1376953/100‑1 (MQ=255)
                                                                              cTTCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATc                 >  1:264210/1‑100 (MQ=255)
                                                                               ttCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCt                >  1:2690623/1‑100 (MQ=255)
                                                                               ttCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCt                >  1:889035/1‑100 (MQ=255)
                                                                                tCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTa               >  1:2984805/1‑100 (MQ=255)
                                                                                tCATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTa               >  2:2896701/1‑100 (MQ=255)
                                                                                 cATCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTAc              <  2:287471/100‑1 (MQ=255)
                                                                                  aTCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTAcc             <  1:264444/100‑1 (MQ=255)
                                                                                   tCGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCg            <  1:2564728/100‑1 (MQ=255)
                                                                                    cGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGc           >  2:1202573/1‑100 (MQ=255)
                                                                                    cGGTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGc           <  1:2978068/100‑1 (MQ=255)
                                                                                     ggTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACcaca          >  1:963691/1‑100 (MQ=255)
                                                                                     ggTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCa          >  2:889932/1‑100 (MQ=255)
                                                                                      gTACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCaa         >  2:1538464/1‑100 (MQ=255)
                                                                                       tACGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCaaa        <  2:569478/100‑1 (MQ=255)
                                                                                        aCGATAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAg       <  1:3163855/100‑1 (MQ=255)
                                                                                         cGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGc      >  1:1069387/1‑100 (MQ=255)
                                                                                           aTAGAGCGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGccc    <  1:2471713/100‑1 (MQ=255)
                                                                                                |                                                                                                
CCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCAC  >  NC_000913/4103673‑4103865

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