Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,103,886 A→G 32.3% Y364H (TAT→CAT)  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,8860AG32.3% 138.3 / 91.4 114Y364H (TAT→CAT) cpxAsensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base A (42/35);  new base G (19/18);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

ATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACG  >  NC_000913/4103791‑4103984
                                                                                               |                                                                                                  
aTCTTCCGTGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGc                                                                                                <  2:210151/100‑1 (MQ=255)
 tCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCa                                                                                               <  1:1617175/100‑1 (MQ=255)
  cTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTTATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCag                                                                                              <  2:1428210/100‑1 (MQ=255)
  cTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCag                                                                                              >  2:921925/1‑100 (MQ=255)
  cTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCag                                                                                              >  1:616701/1‑100 (MQ=255)
  cTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCag                                                                                              >  2:346434/1‑100 (MQ=255)
  cTTCCGGGCTAACGCCAAGACCATCGTCGTCCAACGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCGCTTCAACATTCGTACGGGAATGACGc                                                                                                >  2:1480268/1‑98 (MQ=255)
   ttCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGAAGCaga                                                                                             >  1:621241/1‑100 (MQ=255)
     ccGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGc                                                                                           <  1:970905/100‑1 (MQ=255)
      cGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTCTCTACCGCAAAGCCCACTTCAATCTTCGTATTGGAATAACGCAGAGCa                                                                                          <  1:1413304/100‑1 (MQ=255)
       gggCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCAt                                                                                         >  1:514741/1‑100 (MQ=255)
         gCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTa                                                                                       >  1:1569283/1‑100 (MQ=255)
          cTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTAc                                                                                      >  1:1872985/1‑100 (MQ=255)
            aaCGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGTAATAACGCAGAGCATTACGa                                                                                    >  2:687156/1‑100 (MQ=255)
            aaCGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGa                                                                                    >  2:363761/1‑100 (MQ=255)
            aaCGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGa                                                                                    >  2:1491310/1‑100 (MQ=255)
            aaCGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGa                                                                                    >  2:2147531/1‑100 (MQ=255)
            aaCGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGa                                                                                    >  2:332608/1‑100 (MQ=255)
            aaCGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGa                                                                                    >  1:1137875/1‑100 (MQ=255)
              cGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAAc                                                                                  <  1:1447561/100‑1 (MQ=255)
               gCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGTGAATAACGCAGAGCATTACGAACa                                                                                 >  2:58357/1‑100 (MQ=255)
                      ccATCGTCGTCCACGGTAATTTTTATACCGGCTTTAGCTACCGTAAAGCCCACTGTAATCGTCGTATGGGAATAACGCAGAGCATTACGAACAATGtttt                                                                          >  2:1305263/1‑100 (MQ=255)
                      ccATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGtttt                                                                          >  1:1423034/1‑100 (MQ=255)
                      ccATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGtttt                                                                          >  2:1612674/1‑100 (MQ=255)
                        aTCGTCGTCCACCGTGATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTcc                                                                        <  2:2081164/100‑1 (MQ=255)
                        aTCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTcc                                                                        <  2:454020/100‑1 (MQ=255)
                           gtcgtcCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAgc                                                                     >  1:8155/1‑100 (MQ=255)
                           gtcgtcCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAgc                                                                     >  1:1613649/1‑100 (MQ=255)
                            tcgtcCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAgcg                                                                    <  2:1518592/100‑1 (MQ=255)
                             cgtcCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAgcgc                                                                   >  2:1066931/1‑100 (MQ=255)
                              gtcCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCa                                                                  >  1:1166010/1‑100 (MQ=255)
                                 cACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACtt                                                               >  2:1763601/1‑100 (MQ=255)
                                 cACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACtt                                                               >  1:1132997/1‑100 (MQ=255)
                                  aCCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACATTGTTTTCCAGCGCACttt                                                              >  2:1401671/1‑100 (MQ=255)
                                  aCCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACttt                                                              >  1:1340449/1‑100 (MQ=255)
                                   ccGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTc                                                             >  2:1542656/1‑100 (MQ=255)
                                     gTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCa                                                           >  1:1413470/1‑100 (MQ=255)
                                      tAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCa                                                           >  1:452433/1‑99 (MQ=255)
                                      tAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAg                                                          <  1:921926/100‑1 (MQ=255)
                                      tAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAg                                                          <  1:1168164/100‑1 (MQ=255)
                                      tAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAg                                                          <  1:1568656/100‑1 (MQ=255)
                                      tAATGGTGATACCGTCTATATATACAGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCAAACGCACTTTCCAg                                                          >  2:1445343/1‑100 (MQ=255)
                                        aTGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAgg                                                         <  2:511804/99‑1 (MQ=255)
                                        aTGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAggg                                                        <  2:785364/100‑1 (MQ=255)
                                         tGGTGATACCGTCTTTATCTACCGCAAAGCCGACTTCTATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGc                                                       <  2:1509376/100‑1 (MQ=255)
                                         tGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGc                                                       >  2:1389013/1‑100 (MQ=255)
                                         tGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGc                                                       <  1:1595864/100‑1 (MQ=255)
                                          ggTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCg                                                      <  2:1019992/100‑1 (MQ=255)
                                          ggTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCg                                                      >  2:1481326/1‑100 (MQ=255)
                                          ggTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCg                                                      <  2:1219021/100‑1 (MQ=255)
                                             gctACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCATCTCACTTTCCAGGGCGTTc                                                   <  2:1872973/98‑1 (MQ=255)
                                             gATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTc                                                   >  2:1172101/1‑100 (MQ=255)
                                               tACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCgg                                                 <  2:1263682/100‑1 (MQ=255)
                                                 ccGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGAt                                               <  2:2037109/100‑1 (MQ=255)
                                                  cGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGAtt                                              <  1:759514/100‑1 (MQ=255)
                                                   gTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTg                                             <  2:570750/100‑1 (MQ=255)
                                                    tCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGc                                            <  2:1153435/100‑1 (MQ=255)
                                                     cTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGcc                                           >  2:1014866/1‑100 (MQ=255)
                                                     cTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGcc                                           >  1:1310449/1‑100 (MQ=255)
                                                      tttATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGTATGACGCAGAGCATTACGAACAATGTTTTCCAGCTCACTTTCCAGGGCGTTCGGATTGCCg                                          <  2:730401/100‑1 (MQ=255)
                                                      tttATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCg                                          <  1:161642/100‑1 (MQ=255)
                                                          tCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACa                                      >  1:1293064/1‑100 (MQ=255)
                                                          tCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACa                                      <  2:517181/100‑1 (MQ=255)
                                                             aCCGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCg                                   <  1:1918890/100‑1 (MQ=255)
                                                              ccGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCgg                                  >  1:1052257/1‑100 (MQ=255)
                                                               cGCAAAGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGc                                 >  2:1995142/1‑100 (MQ=255)
                                                                gCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGcc                                <  2:1384112/100‑1 (MQ=255)
                                                                  aaaGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGTGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCAc                              <  1:797405/100‑1 (MQ=255)
                                                                   aaGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACg                             >  2:1250500/1‑100 (MQ=255)
                                                                    aGCCCACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACgg                            <  1:9160/100‑1 (MQ=255)
                                                                     gCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACgg                            >  2:1748899/1‑99 (MQ=255)
                                                                     gCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGc                           <  2:1646941/100‑1 (MQ=255)
                                                                     gCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGc                           >  1:843249/1‑100 (MQ=255)
                                                                      cccACTTCAATCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCagatcg                          >  1:2164991/1‑95 (MQ=255)
                                                                       ccACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGccc                         >  1:2223852/1‑100 (MQ=255)
                                                                       ccACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGccc                         >  1:623902/1‑100 (MQ=255)
                                                                       ccACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGccc                         >  1:1957344/1‑100 (MQ=255)
                                                                           ttCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCa                        >  1:501368/1‑97 (MQ=255)
                                                                            tCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCg                    <  2:2158895/100‑1 (MQ=255)
                                                                            tCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCg                    <  1:775209/100‑1 (MQ=255)
                                                                             cAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCgg                   <  2:6345/100‑1 (MQ=255)
                                                                              aaTCTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCggg                  <  1:1389070/100‑1 (MQ=255)
                                                                              aaTCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCggg                  <  1:710477/100‑1 (MQ=255)
                                                                                 cTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAg               <  1:363760/100‑1 (MQ=255)
                                                                                 cTTCGTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAg               <  1:1491320/100‑1 (MQ=255)
                                                                                  ttCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGt              <  2:1005825/100‑1 (MQ=255)
                                                                                  ttCGTATGGGAATAACGCAGACCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCCGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAg               <  2:1915132/99‑1 (MQ=255)
                                                                                    cGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTa            <  2:1464436/100‑1 (MQ=255)
                                                                                     gTATGGGAATGACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTaa           >  1:395347/1‑100 (MQ=255)
                                                                                      tATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAAc          <  2:314545/100‑1 (MQ=255)
                                                                                      tATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAAc          >  2:1739767/1‑100 (MQ=255)
                                                                                        tGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTg        >  2:1415186/1‑100 (MQ=255)
                                                                                        tGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTg        >  2:1691056/1‑100 (MQ=255)
                                                                                         gggAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGt       <  1:1415196/100‑1 (MQ=255)
                                                                                         gggAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGt       >  1:313699/1‑100 (MQ=255)
                                                                                          ggAATGACGCAGAGTATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTc      <  2:2081949/100‑1 (MQ=255)
                                                                                          ggAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTc      >  1:1646247/1‑100 (MQ=255)
                                                                                             aTAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAAc   <  1:965240/100‑1 (MQ=255)
                                                                                              tAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACg  >  1:1986495/1‑100 (MQ=255)
                                                                                              tAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACg  <  1:992766/100‑1 (MQ=255)
                                                                                               |                                                                                                  
ATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATACCGTCTTTATCTACCGCAAAGCCCACTTCAATCTTCGTATGGGAATAACGCAGAGCATTACGAACAATGTTTTCCAGCGCACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACG  >  NC_000913/4103791‑4103984

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: