Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,981,889 A→C 6.4% Q845P (CAG→CCG)  PP_5242 → GAF domain‑containing sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,981,8890AC 6.4% 200.2 / 4.6 63Q845P (CAG→CCG) PP_5242GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base C (2/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

CGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGA  >  NC_002947/5981744‑5982036
                                                                                                                                                 |                                                                                                                                                     
cgctcgGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGt                                                                                                                                                    >  8:58287/1‑149 (MQ=255)
 gctcgGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTc                                                                                                                                                   <  8:49249/149‑1 (MQ=255)
 gctcgGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTc                                                                                                                                                   >  8:209552/1‑149 (MQ=255)
     ggtggtGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACCTGCTGCCAGTCGAgg                                                                                                                                               <  6:327133/149‑1 (MQ=255)
     ggTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAgg                                                                                                                                               >  1:121107/1‑149 (MQ=255)
      gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGc                                                                                                                                              >  4:218653/1‑149 (MQ=255)
      gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACCTGCTGCCAGTCGAGGc                                                                                                                                              <  2:184209/149‑1 (MQ=255)
      gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGCTAGCAGCCCGGTGCGGGGCCTCGACTCGCAGCAGGTCGAGGc                                                                                                                                              <  8:119891/149‑1 (MQ=255)
      gtgtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGAACAACTTGCGCGGTGACTATCAGACGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGCGCCTCGACTGGCACCAGGTCGAGGc                                                                                                                                              >  3:152180/1‑149 (MQ=255)
        gtTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGc                                                                                                                                              >  5:146339/1‑147 (MQ=255)
            gACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTTTCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAg                                                                                                                                                >  6:65310/1‑141 (MQ=255)
            gACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTTTCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAg                                                                                                                                                <  5:65310/141‑1 (MQ=255)
            gACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCGTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCa                                                                                                                                        >  1:61611/1‑149 (MQ=255)
              ccAGGCACGGGCGTGCCGGCTGCAGGCGATGCTGTACACCTTGCGCGGTGACTATCAGACGGGCATCGACAGCGCGCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATg                                                                                                                                      <  4:152180/149‑1 (MQ=255)
               cAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGc                                                                                                                                     <  1:137144/149‑1 (MQ=255)
                 ggCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGcc                                                                                                                                    >  3:147974/1‑148 (MQ=255)
                         cGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAg                                                                                                                         <  4:131074/149‑1 (MQ=255)
                                 cTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGAt                                                                                                                 <  3:218653/149‑1 (MQ=255)
                                  agCAGGCCATGCTGTACACCATGCGCCTAGACTATCAGGCGGCCATCGACAGCGCCCTGTCCGGCCTCGACTTACTGGGTATCAGCCCGGTGCGGGCCCTCGACTGCCACCAGGTCTAGGCCAGTCATGCACATGTTCAGGCGTTGCTc                                                                                                                >  8:154997/2‑149 (MQ=255)
                                    cAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATc                                                                                                                <  6:146339/147‑1 (MQ=255)
                                     aGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGaa                                                                                                             <  4:188217/149‑1 (MQ=255)
                                      ggCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAAc                                                                                                            <  2:264078/149‑1 (MQ=255)
                                       gcgcTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACa                                                                                                           <  4:147974/149‑1 (MQ=255)
                                              tgTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGAc                                                                                                    >  5:58909/1‑149 (MQ=255)
                                               gTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGAcc                                                                                                   >  8:255868/1‑149 (MQ=255)
                                                tACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCg                                                                                                  >  5:18976/1‑149 (MQ=255)
                                                tACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCg                                                                                                  <  6:58909/149‑1 (MQ=255)
                                                  cacCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCgg                                                                                                 <  7:209552/148‑1 (MQ=255)
                                                   acCTTGCGCGGTGACTTTCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGc                                                                                                                       <  1:307795/125‑1 (MQ=255)
                                                   acCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGc                                                                                                                       >  2:307795/1‑125 (MQ=255)
                                                     cTTGCGCGGTGACTATCAGGCTGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCAtt                                                                                             <  7:154997/149‑1 (MQ=255)
                                                     cTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCAtt                                                                                             >  4:93006/1‑149 (MQ=255)
                                                          gcgGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTg                                                                                        >  8:10529/1‑149 (MQ=255)
                                                          gcgGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTg                                                                                        <  7:197/149‑1 (MQ=255)
                                                          gcgGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTg                                                                                        >  8:195573/1‑149 (MQ=255)
                                                            ggTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                                      <  7:255868/149‑1 (MQ=255)
                                                            ggTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                                      <  6:18976/149‑1 (MQ=255)
                                                            ggTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                                      <  8:278237/149‑1 (MQ=255)
                                                                            ccATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                                      >  1:47584/1‑133 (MQ=255)
                                                                            ccATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGa                                                                                      <  2:47584/133‑1 (MQ=255)
                                                                               tCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGTGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGa                                                                    >  6:83596/1‑148 (MQ=255)
                                                                               tCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGAc                                                                   <  7:195573/149‑1 (MQ=255)
                                                                                  aCAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCg                                                                >  2:66684/1‑149 (MQ=255)
                                                                                    aGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCgga                                                              <  2:121107/149‑1 (MQ=255)
                                                                                    aGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCgga                                                              <  4:207956/149‑1 (MQ=255)
                                                                                     gcgcCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCggag                                                             <  5:43494/149‑1 (MQ=255)
                                                                                               cggccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACgg                                                   >  7:8066/1‑149 (MQ=255)
                                                                                                 gccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGtt                                                 <  2:61611/149‑1 (MQ=255)
                                                                                                   cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGc                                               >  1:119112/1‑149 (MQ=255)
                                                                                                   cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGc                                               <  2:119112/149‑1 (MQ=255)
                                                                                                       aCTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGTGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATc                                           <  5:83596/149‑1 (MQ=255)
                                                                                                           aCTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGa                                                                                                                  >  6:260980/1‑74 (MQ=255)
                                                                                                           aCTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGa                                                                                                                  <  5:260980/74‑1 (MQ=255)
                                                                                                             tGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCct                                      >  1:253297/1‑149 (MQ=255)
                                                                                                               ggTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCGCATGTGAAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGACACTGCCGCGCCCCGAAACGGAGGAAGTGACGTTTGCCATCACCctgc                                    >  1:82700/1‑149 (MQ=255)
                                                                                                                       cccGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCAc                             <  5:162041/147‑1 (MQ=255)
                                                                                                                       cccGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCAc                             >  6:162041/1‑147 (MQ=255)
                                                                                                                           gTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGtc                       <  3:93006/149‑1 (MQ=255)
                                                                                                                           gTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGt                        >  1:247381/1‑148 (MQ=255)
                                                                                                                               ggggCCTCGACTGGGAGGAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGATGGTCAGGCATTGCCTGGATTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGtcatca                   <  2:111032/149‑1 (MQ=255)
                                                                                                                                   ccTCGACCTGCTGCCAGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCct               >  7:100915/1‑149 (MQ=255)
                                                                                                                                   ccTCGACCTGCTGCCAGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCct               >  8:43287/1‑149 (MQ=255)
                                                                                                                                        aCTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCttcttc          <  4:306872/149‑1 (MQ=255)
                                                                                                                                            gcagcaGGTCGAGCCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCAATGCCTGGAGTCACTGCCGCGCACCGACTCGGACGAAGTGAAGGTTGCCATCAGCCTGCTGGACACTCTCTCATCATCCTTCATCGTg       >  8:265966/1‑148 (MQ=255)
                                                                                                                                                cagGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGga  <  1:66684/149‑1 (MQ=255)
                                                                                                                                                 agGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGc                                <  5:175496/118‑1 (MQ=255)
                                                                                                                                                 agGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGc                                >  6:175496/1‑118 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGCTCGGTGTTGGACCAGGCACGGGCGTGCCGGCTGCAGGCGCTGCTGTACACCTTGCGCGGTGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCATTGCCTGGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGA  >  NC_002947/5981744‑5982036

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