Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,918,259 C→G 6.3% 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.3% 251.9 / 6.8 79intergenic (‑35/+83)PP_1718/prcEAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base G (2/3);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

GCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGC  >  NC_002947/1918117‑1918404
                                                                                                                                              |                                                                                                                                                 
gCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGccc                                                                                                                                             >  8:61091/1‑149 (MQ=255)
  cGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGGGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCgg                                                                                                                                           >  3:310394/1‑149 (MQ=255)
  cGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCgg                                                                                                                                           >  8:261436/1‑149 (MQ=255)
  cGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCgg                                                                                                                                           <  6:222789/149‑1 (MQ=255)
   gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGc                                                                                                                                          <  7:261436/149‑1 (MQ=255)
   gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGc                                                                                                                                          <  3:93869/149‑1 (MQ=255)
      tACCCTATGTTAACTCATGCACATGATTGGATGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCt                                                                                                                                       <  8:159464/149‑1 (MQ=255)
         ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCa                                                                                                                                     >  3:126383/1‑148 (MQ=255)
         ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAAAGGTCATGATGGATTACTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGACCTCa                                                                                                                                     >  3:17196/1‑148 (MQ=255)
                        gCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCAc                                                                                                                     >  2:123761/1‑149 (MQ=255)
                           caTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                         <  7:182834/126‑1 (MQ=255)
                           caTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                         >  8:182834/1‑126 (MQ=255)
                           caTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                  >  1:322709/1‑149 (MQ=255)
                           caTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                  <  1:123761/149‑1 (MQ=255)
                            aTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                  >  4:251629/1‑148 (MQ=255)
                            aTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                  >  7:125892/1‑148 (MQ=255)
                               aTTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGAc                                                                                                               <  1:100677/148‑1 (MQ=255)
                               aTTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACt                                                                                                              <  4:139455/149‑1 (MQ=255)
                               aTTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATAATTCACGACAGCACGATGAc                                                                                                               <  3:275926/148‑1 (MQ=255)
                                ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCCAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGa                                                                                                                   <  4:105494/143‑1 (MQ=255)
                                ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGa                                                                                                                   >  3:105494/1‑143 (MQ=255)
                                ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACt                                                                                                              <  3:295449/148‑1 (MQ=255)
                                    cTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCTGAAACGGTCATGATGGATTCCTGAAAAGTAAGGGACGGGCATCTGCATGCCCGGCCGGCATATTTCCCGCCAGCATGATGACTGatt                                                                                                          >  5:278561/1‑148 (MQ=255)
                                     tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtga                                                                                                        >  8:20009/1‑149 (MQ=255)
                                     tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtga                                                                                                        >  4:54579/1‑149 (MQ=255)
                                            gACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGa                                                                                                 <  3:302849/149‑1 (MQ=255)
                                              cTGTGAATAATGCGCTGAGCCAGGATGGCACTAAAGGCAGTGAGCTGTTCGGCAACTGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGAt                                                                                                <  4:17196/148‑1 (MQ=255)
                                                      aCTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACattat                                                                                        >  4:299773/1‑148 (MQ=255)
                                                        tGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCt                                                                                     >  6:326847/1‑149 (MQ=255)
                                                        tGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCt                                                                                     <  8:125892/149‑1 (MQ=255)
                                                        tGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCt                                                                                     >  6:331692/1‑149 (MQ=255)
                                                           gcTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGcc                                                                                  <  7:61091/149‑1 (MQ=255)
                                                            cTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCca                                                                                 <  4:126383/149‑1 (MQ=255)
                                                                 ccAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCAtt                                                                            <  6:215727/149‑1 (MQ=255)
                                                                       tGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAg                                                                      <  3:54579/149‑1 (MQ=255)
                                                                         gCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGt                                                                    >  7:77626/1‑149 (MQ=255)
                                                                               aaGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTa                                                              >  8:297092/1‑149 (MQ=255)
                                                                               aaGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTACCGCCAGCACGATGACTGTTTCATGACACAGCACATAATTCTGACATCATTGCTTAGGTGCCtt                                                               >  7:59850/1‑148 (MQ=255)
                                                                                aGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTAc                                                             >  5:65838/1‑149 (MQ=255)
                                                                                        ggCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                       <  4:310394/146‑1 (MQ=255)
                                                                                        aGCTGTTCTGCAACGGTCATTATGGATTCCTGAAAATTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                       <  8:201854/147‑1 (MQ=255)
                                                                                        aGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCAAATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGc                                                                  <  3:321021/136‑1 (MQ=255)
                                                                                        aGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCAAATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGc                                                                  >  4:321021/1‑136 (MQ=255)
                                                                                          cTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCa                                                   <  3:226341/149‑1 (MQ=255)
                                                                                             ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCt                                                                <  2:239745/133‑1 (MQ=255)
                                                                                             ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCt                                                                >  1:239745/1‑133 (MQ=255)
                                                                                                  cAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTg                                           <  8:77626/149‑1 (MQ=255)
                                                                                                    aCGGTCATTATGTATTCCTTAAATTTATGGGCCGGGCATGTGCATGCCGGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCt                                         <  2:322709/149‑1 (MQ=255)
                                                                                                     cGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                        <  7:50745/149‑1 (MQ=255)
                                                                                                     cGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                        <  6:65838/149‑1 (MQ=255)
                                                                                                     cGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                        <  7:20009/149‑1 (MQ=255)
                                                                                                       gTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTgg                                       >  6:253923/1‑148 (MQ=255)
                                                                                                       gTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGc                                      <  5:326847/149‑1 (MQ=255)
                                                                                                       gTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGc                                      <  5:331692/149‑1 (MQ=255)
                                                                                                              tgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                <  3:299773/148‑1 (MQ=255)
                                                                                                                 aTTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAg                            >  2:132792/1‑149 (MQ=255)
                                                                                                                     cTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                        <  3:269536/133‑1 (MQ=255)
                                                                                                                     cTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                        >  4:269536/1‑133 (MQ=255)
                                                                                                                      tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGGCTTACTTCTTGGCCACGTTGGTGCTGGCACTGAGGTAGTCGAt                       >  4:185459/1‑149 (MQ=255)
                                                                                                                      tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCGCTGAGGTAGTCGAt                       >  3:59563/1‑149 (MQ=255)
                                                                                                                      tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGa                        >  1:41916/1‑148 (MQ=255)
                                                                                                                      tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                       >  6:295758/1‑149 (MQ=255)
                                                                                                                      tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                       >  3:149992/1‑149 (MQ=255)
                                                                                                                         aaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAg                                                                      <  4:259183/99‑1 (MQ=255)
                                                                                                                         aaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAg                                                                      >  3:259183/1‑99 (MQ=255)
                                                                                                                         aaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGt                                                                    <  1:212410/101‑1 (MQ=255)
                                                                                                                         aaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATAATTGCTGAGGt                                                                    >  2:212410/1‑101 (MQ=255)
                                                                                                                             tAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATg                >  2:80445/1‑149 (MQ=255)
                                                                                                                              aaGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGc               >  7:179954/1‑149 (MQ=255)
                                                                                                                                gggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAgg                   >  4:260828/1‑143 (MQ=255)
                                                                                                                                gggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAgg                   <  3:260828/143‑1 (MQ=255)
                                                                                                                                 ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGc               >  6:315300/1‑146 (MQ=255)
                                                                                                                                 ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGc               <  5:315300/146‑1 (MQ=255)
                                                                                                                                  gCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTAATTATTGGCCACCTTGGTGCTGGCTCTGAGGTAGTCGATCAGGATGCGccc           >  4:176722/1‑149 (MQ=255)
                                                                                                                                  gCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                       >  6:341690/1‑137 (MQ=255)
                                                                                                                                  gCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                       <  5:341690/137‑1 (MQ=255)
                                                                                                                                    cGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCCACCgg         >  4:116731/1‑149 (MQ=255)
                                                                                                                                    cGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCAACATTGCTGAGGTGCCTTAATTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCgg         >  1:109780/1‑149 (MQ=255)
                                                                                                                                         aTCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTcc    >  3:22513/1‑149 (MQ=255)
                                                                                                                                              cATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCcgc  <  2:41916/146‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGC  >  NC_002947/1918117‑1918404

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: