Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,095 C→G 100% L594F (TTG→TTC rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,318,0950CG100.0% 283.6 / NA 76L594F (TTG→TTCrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (37/39);  total (37/39)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCAT  >  NC_000913/2317949‑2318232
                                                                                                                                                  |                                                                                                                                         
aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGaa                                                                                                                                         <  2:762505/149‑1 (MQ=255)
aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGaa                                                                                                                                         <  3:545745/149‑1 (MQ=255)
 ccTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAAc                                                                                                                                        <  1:724808/149‑1 (MQ=255)
  cTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACa                                                                                                                                       <  3:128584/149‑1 (MQ=255)
   tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAg                                                                                                                                      >  2:557704/1‑149 (MQ=255)
   tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAg                                                                                                                                      >  3:41708/1‑149 (MQ=255)
     gaagaaGGGATCAAACAAGCGCCCCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGt                                                                                                                                    <  3:587218/149‑1 (MQ=255)
        gaagGGATCAACCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTgg                                                                                                                                 <  1:230046/149‑1 (MQ=255)
        gaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTgg                                                                                                                                 >  1:238321/1‑149 (MQ=255)
         aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                <  4:351752/149‑1 (MQ=255)
         aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                >  1:265965/1‑149 (MQ=255)
           gggATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGgaga                                                                                                                              <  2:454237/149‑1 (MQ=255)
            ggATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGAt                                                                                                                             <  2:519198/149‑1 (MQ=255)
            ggATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGAt                                                                                                                             <  3:189983/149‑1 (MQ=255)
             gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATg                                                                                                                            <  3:941547/149‑1 (MQ=255)
               tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGAc                                                                                                                          <  1:921234/149‑1 (MQ=255)
                 aaaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                <  1:389855/140‑1 (MQ=255)
                 aaaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGCCGTTACTGAACAGGTTGga                                                                                                                                >  2:389855/1‑140 (MQ=255)
                    cAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgct                                                                                                                     <  2:558901/149‑1 (MQ=255)
                    cAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgct                                                                                                                     <  3:480773/149‑1 (MQ=255)
                     aaGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgctg                                                                                                                    >  4:594189/1‑149 (MQ=255)
                      aGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATTACCTGCTGt                                                                                                                   <  1:543111/149‑1 (MQ=255)
                         gcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaaa                                                                                                                <  1:551284/149‑1 (MQ=255)
                             cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTg                                                                                                                                  >  4:188356/1‑127 (MQ=255)
                             cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTg                                                                                                                                  <  3:188356/127‑1 (MQ=255)
                             cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCa                                                                                                            >  1:95286/1‑149 (MQ=255)
                                ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCg                                                                                                         >  1:898275/1‑149 (MQ=255)
                                ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCg                                                                                                         >  4:686458/1‑149 (MQ=255)
                                ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCg                                                                                                         >  2:508834/1‑149 (MQ=255)
                                ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCg                                                                                                         >  1:1042317/1‑149 (MQ=255)
                                             aTCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa                                                                                            <  4:879164/149‑1 (MQ=255)
                                              tCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAg                                                                                           <  4:408506/149‑1 (MQ=255)
                                               ccccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAgg                                                                                          >  2:926999/1‑149 (MQ=255)
                                                cccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCGGTAAACGCATCGGGTCGCCATGTAAGGc                                                                                         >  2:4653/1‑149 (MQ=255)
                                                   aCGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCAc                                                                                      >  1:458050/1‑149 (MQ=255)
                                                     gCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTg                                                                                    >  4:989444/1‑149 (MQ=255)
                                                       cGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGc                                                                                  <  4:389498/149‑1 (MQ=255)
                                                            tCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                             >  2:461666/1‑149 (MQ=255)
                                                                gAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGt                                                                         <  1:615839/149‑1 (MQ=255)
                                                                 aaCACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGtt                                                                        <  1:926999/149‑1 (MQ=255)
                                                                    acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                     <  3:989444/149‑1 (MQ=255)
                                                                      ggATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGt                                                                                                                   >  4:269860/1‑101 (MQ=255)
                                                                      ggATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGt                                                                                                                   <  3:269860/101‑1 (MQ=255)
                                                                           gagagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCa                                                              <  2:613699/149‑1 (MQ=255)
                                                                           gagagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCa                                                              <  4:869398/149‑1 (MQ=255)
                                                                            agtgaTAATCGCCATCCGCGCCAAGATGCAAAACTATACAGCCGGTATCGGTGAACTTTATGGCGTTACTGTACAGGTTCGAGATCACCTGCTGTAAACGCATCGTGTCGCCGTTTAAGGCCAGTGGCACATCCGGTTCAATAAAGCAg                                                             >  4:653249/4‑149 (MQ=255)
                                                                              agaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                             >  1:25577/1‑131 (MQ=255)
                                                                              agaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                             <  2:25577/131‑1 (MQ=255)
                                                                              agaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTa                                                           >  2:765699/1‑149 (MQ=255)
                                                                              agaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTa                                                           >  2:103276/1‑149 (MQ=255)
                                                                               gaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTAc                                                          >  2:501483/1‑149 (MQ=255)
                                                                                  aaTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAg                                                       >  4:829868/1‑149 (MQ=255)
                                                                                       ccATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTaa                                                  >  4:178987/1‑149 (MQ=255)
                                                                                        cATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAAc                                                 <  3:186930/149‑1 (MQ=255)
                                                                                                 gAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGt                                        >  3:257068/1‑149 (MQ=255)
                                                                                                         aaaaCTATACAGCCGGTATCGGTGAATTTGATGGCGTTAGTGAACAGGTTGGAGATGACCTGCTCTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCATAACTGCTTGCGTACCACCAg                                >  4:1030915/1‑149 (MQ=255)
                                                                                                           aaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCg                              <  3:829868/149‑1 (MQ=255)
                                                                                                             cTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGt                            <  4:179474/149‑1 (MQ=255)
                                                                                                              tataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTa                           >  4:109314/1‑149 (MQ=255)
                                                                                                               ataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTaa                          >  2:464831/1‑149 (MQ=255)
                                                                                                                taCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTaaa                         >  1:985365/1‑149 (MQ=255)
                                                                                                                taCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTaaa                         <  4:970440/149‑1 (MQ=255)
                                                                                                                taCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCAGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTaaa                         <  2:588889/149‑1 (MQ=255)
                                                                                                                   aGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                      >  4:94013/1‑149 (MQ=255)
                                                                                                                    gCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGt                     >  1:608423/1‑149 (MQ=255)
                                                                                                                        gTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGc                 <  4:41708/149‑1 (MQ=255)
                                                                                                                              gTGAATTTTATGGCGTTACTGATCAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt           >  3:857558/1‑149 (MQ=255)
                                                                                                                              gTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt           >  1:244482/1‑149 (MQ=255)
                                                                                                                              gTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt           >  2:815853/1‑149 (MQ=255)
                                                                                                                                gAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAtgt         <  3:171510/149‑1 (MQ=255)
                                                                                                                                   ttttATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGt      <  2:873889/149‑1 (MQ=255)
                                                                                                                                       aTGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGtt                                                                        <  1:39590/79‑1 (MQ=255)
                                                                                                                                       aTGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGtt                                                                        >  2:39590/1‑79 (MQ=255)
                                                                                                                                       aTGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtcat  <  1:990049/149‑1 (MQ=255)
                                                                                                                                           cGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTgg       <  3:983351/140‑1 (MQ=255)
                                                                                                                                           cGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTgg       >  4:983351/1‑140 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
ACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCAT  >  NC_000913/2317949‑2318232

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