Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,918,259 C→G 6.4% intergenic (‑35/+83) PP_1718 ← / ← prc EAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,918,2590CG 6.4% 202.6 / 3.8 62intergenic (‑35/+83)PP_1718/prcEAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base G (2/2);  total (28/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 = 3.75e-01

CTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGAG  >  NC_002947/1918111‑1918407
                                                                                                                                                    |                                                                                                                                                    
cTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGc                                                                                                                                                      <  6:303109/149‑1 (MQ=255)
 tCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCa                                                                                                                                                     >  5:311032/1‑149 (MQ=255)
  cATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCAt                                                                                                                                                    >  6:7993/1‑149 (MQ=255)
    tGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGc                                                                                                                                                  >  4:236754/1‑149 (MQ=255)
        cGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCgg                                                                                                                                              <  7:185212/149‑1 (MQ=255)
         gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGc                                                                                                                                             <  4:107779/149‑1 (MQ=255)
               ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCa                                                                                                                                        >  6:6745/1‑148 (MQ=255)
                cTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                      >  5:322146/1‑149 (MQ=255)
                    gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCGCATATTTc                                                                                                                                  >  8:345775/1‑149 (MQ=255)
                     ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGTTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCTTACTTTTCa                                                                                                                                 <  1:107946/149‑1 (MQ=38)
                     ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTc                                                                                                                                  <  3:236754/148‑1 (MQ=255)
                     ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCa                                                                                                                                 <  1:91396/149‑1 (MQ=255)
                     ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCa                                                                                                                                 <  2:32487/149‑1 (MQ=255)
                        aCTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGa                                                                                                                              >  8:15390/1‑149 (MQ=255)
                                 caTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTTAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                     <  2:266857/149‑1 (MQ=255)
                                 caTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                     <  2:325907/149‑1 (MQ=255)
                                      ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACt                                                                                                                 <  5:7993/148‑1 (MQ=255)
                                          cTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtg                                                                                                            <  7:186612/149‑1 (MQ=255)
                                              aaCAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                        <  1:301064/149‑1 (MQ=255)
                                                 aGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                        >  3:274040/1‑146 (MQ=255)
                                                 aGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                        >  1:242957/1‑146 (MQ=255)
                                                   aCTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                        >  2:290286/1‑144 (MQ=255)
                                                       tgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                               <  7:15390/149‑1 (MQ=255)
                                                       tgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                               >  7:252943/1‑149 (MQ=255)
                                                         aaTACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTTTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACat                                                                                              <  5:191140/148‑1 (MQ=255)
                                                          aTACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACatta                                                                                            <  3:175857/149‑1 (MQ=255)
                                                                  cTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCca                                                                                    <  8:81583/149‑1 (MQ=255)
                                                                     aGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCTGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCcatc                                                                                  <  2:117759/148‑1 (MQ=255)
                                                                                  cTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCACATAAATCACGACAGCACGATCACTGTATGATGACAGATCACATTATTCTGCCATCATAGCTGAGGTGcc                                                                    >  4:153901/1‑149 (MQ=255)
                                                                                  cTCAAGGCACTGAGCTGTTCGGCAAAGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGcc                                                                    >  7:194604/1‑149 (MQ=255)
                                                                                    ctccgctcTGCGCTGGCCGGCAACGCGCAGGCTGGCTTCCTGAAAAGTAAGCGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGATGTGcc                                                                    <  4:197203/140‑1 (MQ=255)
                                                                                     aaGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGTCAGATCACATTATTCTGCCATCATTGCTGAGGTGCCt                                                                   <  2:76951/147‑1 (MQ=255)
                                                                                     aaGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCt                                                                   <  2:242957/147‑1 (MQ=255)
                                                                                            tGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                          <  3:295441/149‑1 (MQ=255)
                                                                                            tGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                          >  7:210351/1‑149 (MQ=255)
                                                                                              aGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                          <  5:6745/147‑1 (MQ=255)
                                                                                              aGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACACCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                          <  6:67040/147‑1 (MQ=255)
                                                                                              aGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAACGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATAACAATATTCTGCCATCATTGCTGAGGTgac                                                                    >  3:197203/1‑135 (MQ=255)
                                                                                                   ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCAc                                                     <  3:258073/147‑1 (MQ=255)
                                                                                                      ggCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTg                                                 <  1:290286/148‑1 (MQ=255)
                                                                                                       gCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGt                                               <  8:252943/149‑1 (MQ=255)
                                                                                                       gCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGt                                               >  7:31218/1‑149 (MQ=255)
                                                                                                         aaCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGc                                             >  2:34670/1‑149 (MQ=255)
                                                                                                          aCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCt                                            <  8:194604/149‑1 (MQ=255)
                                                                                                            ggTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTgg                                          >  4:356302/1‑149 (MQ=255)
                                                                                                                 tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACt                                      >  3:33565/1‑148 (MQ=255)
                                                                                                                  gatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACt                                      >  3:163331/1‑147 (MQ=255)
                                                                                                                   atgGATTCCTTAAAAGTAAGGGCCGGGCATCTGCATGCCCGTCCCTCATATTTCACGACAGCACGATGCCTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                   <  2:85221/149‑1 (MQ=255)
                                                                                                                   atgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                    <  1:34670/148‑1 (MQ=255)
                                                                                                                   atgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                    <  8:31218/148‑1 (MQ=255)
                                                                                                                   atgGATTCCTGAAAAGTAAGGGCCCGGCATCTGCATGCCCGGCCCACATATTTAACCACACCAAGATGACTCTTTGATGACAGAACACATTATTATGCCATAATTGCTGAGGTTCCTTACTTCTTGGACACCTTGGTGCTGGCACTGAg                                   >  2:170430/1‑149 (MQ=255)
                                                                                                                    tgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                   <  4:345052/148‑1 (MQ=255)
                                                                                                                            tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                          >  7:39891/1‑149 (MQ=255)
                                                                                                                            tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATAACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGACACCTTGGTGCTGGCACTGAGGTAGTCGa                           >  7:142355/1‑148 (MQ=255)
                                                                                                                              aaaaGTAAGGGCCGGGCAGCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                          <  6:311032/147‑1 (MQ=255)
                                                                                                                               aaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCa                        <  3:356302/148‑1 (MQ=255)
                                                                                                                                  gTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCGTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                   >  3:100012/1‑134 (MQ=255)
                                                                                                                                   tAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                           >  8:272011/1‑125 (MQ=255)
                                                                                                                                      gggCCTGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGa                <  4:274040/149‑1 (MQ=255)
                                                                                                                                        gCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc              >  2:107132/1‑149 (MQ=255)
                                                                                                                                        gCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc              >  8:152608/1‑149 (MQ=255)
                                                                                                                                          cGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCgg            >  5:199212/1‑149 (MQ=255)
                                                                                                                                                 cTGCATGCCGGGCCCTCATATTTCACGACAGCACGATGTCTGTTTGATGACAGATCTCATTATTCTGCCATGATTGCTGAGGTGGCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTGcg      <  3:153901/148‑1 (MQ=255)
                                                                                                                                                    cATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGATGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGa   <  6:109192/148‑1 (MQ=255)
                                                                                                                                                    cATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGAg  <  8:210351/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGAG  >  NC_002947/1918111‑1918407

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