Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,762 A→C 6.6% L394R (CTG→CGG)  PP_2413 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,759,7620AC 6.6% 185.4 / 4.9 61L394R (CTG→CGG) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base C (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

AGTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCG  >  NC_002947/2759615‑2759899
                                                                                                                                                     |                                                                                                                                          
aGTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAg                                                                                                                                           >  5:100615/1‑149 (MQ=255)
 gTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGc                                                                                                                                          <  3:313150/149‑1 (MQ=255)
        gTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCgg                                                                                                                                   >  4:174871/1‑149 (MQ=255)
         tCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGt                                                                                                                                  <  1:247655/149‑1 (MQ=255)
             ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              <  1:207172/149‑1 (MQ=255)
             ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              <  2:134840/149‑1 (MQ=255)
             ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              >  2:207172/1‑149 (MQ=255)
             ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                              <  6:100615/149‑1 (MQ=255)
             ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGACCGATCCGCTGACCGGCCTgg                                                                                                                              <  7:342949/149‑1 (MQ=255)
              cACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTggg                                                                                                                             >  2:329303/1‑149 (MQ=255)
                      gtCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGa                                                                                                                     >  6:136231/1‑149 (MQ=255)
                       tCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAt                                                                                                                    >  6:14338/1‑149 (MQ=255)
                          aTGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGa                                                                                                                 <  7:118703/149‑1 (MQ=255)
                              cTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccg                                                                                                             >  4:338985/1‑149 (MQ=255)
                                aTCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggc                                                                                                           >  2:286090/1‑149 (MQ=255)
                                 tCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                          <  8:336312/149‑1 (MQ=255)
                                            cTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                               >  4:334421/1‑149 (MQ=255)
                                               gAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGc                                                                                                                                          <  4:385823/103‑1 (MQ=255)
                                               gAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGc                                                                                                                                          >  3:385823/1‑103 (MQ=255)
                                                  gTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTgcagca                                                                                         <  3:338985/149‑1 (MQ=255)
                                                     gTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAg                                                                                       <  5:124505/148‑1 (MQ=255)
                                                       cGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGtt                                                                                     >  3:391063/1‑148 (MQ=255)
                                                             ccAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggccg                                                                                                         <  3:276248/122‑1 (MQ=255)
                                                             ccAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggccg                                                                                                         >  4:276248/1‑122 (MQ=255)
                                                              cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAtt                                                                                                                   >  5:275212/1‑111 (MQ=255)
                                                              cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAtt                                                                                                                   <  6:275212/111‑1 (MQ=255)
                                                              cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                          <  3:348731/120‑1 (MQ=255)
                                                              cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                          >  4:348731/1‑120 (MQ=255)
                                                              cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTg                                                                             >  5:15726/1‑148 (MQ=255)
                                                                    aaTTGCTGCAAACTCCCGCTTTTCACCCTCCAGTAACGACATTACCATATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAg                                                                        <  1:136555/149‑1 (MQ=255)
                                                                    aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                               <  6:192481/125‑1 (MQ=255)
                                                                    aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                               >  5:192481/1‑125 (MQ=255)
                                                                    aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGc                                                                       >  3:382116/1‑148 (MQ=255)
                                                                          tGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGc                                                                       <  8:352858/142‑1 (MQ=255)
                                                                          tGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAAcgcgc                                                                >  6:245763/1‑149 (MQ=255)
                                                                          tGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAAcgcgc                                                                >  5:18978/1‑149 (MQ=255)
                                                                          tGCAAACTCCCGCTATTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGc                                                                       >  7:352858/1‑142 (MQ=255)
                                                                                 tCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAg                                                         >  1:271252/1‑149 (MQ=255)
                                                                                 tCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAg                                                         >  2:372269/1‑149 (MQ=255)
                                                                                  cccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCtt                                                            <  6:40915/145‑1 (MQ=255)
                                                                                  cccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCtt                                                            >  5:40915/1‑145 (MQ=255)
                                                                                  cccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGc                                                        >  6:384422/1‑149 (MQ=255)
                                                                                   ccGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCt                                                       >  6:280431/1‑149 (MQ=255)
                                                                                   aggattttCATCCTCAAGTTACGAGA‑‑GGCTTATTTCCAGGCTTCGGCGTTTGCCCAGGCCTGTCAGCGGATCGTTCTGGGCGTCCAGATTGAGGCGACCGATACGCTGTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTTAg                                                         <  8:130876/143‑1 (MQ=255)
                                                                                    cGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTTCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAAAAGCAACGCGCGCTTGAGCTc                                                      >  5:387912/1‑149 (MQ=255)
                                                                                    cGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTc                                                      >  4:282246/1‑149 (MQ=255)
                                                                                     gCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCg                                                     <  5:14338/149‑1 (MQ=255)
                                                                                         ttCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGgc                                                   <  3:174871/147‑1 (MQ=255)
                                                                                             gCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTc                                             <  6:18978/149‑1 (MQ=255)
                                                                                                  cAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGaa                                          <  5:280431/147‑1 (MQ=255)
                                                                                                  cAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCGCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAg                                         <  2:69013/148‑1 (MQ=255)
                                                                                                     tAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGCGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTAc                                      <  1:286090/149‑1 (MQ=255)
                                                                                                         gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGctc                                                                                                 >  4:64793/1‑86 (MQ=255)
                                                                                                         gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGctc                                                                                                 <  3:64793/86‑1 (MQ=255)
                                                                                                         gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGc                                 >  1:335701/1‑149 (MQ=255)
                                                                                                                          aGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCaa                  <  6:15726/149‑1 (MQ=255)
                                                                                                                                   cggTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTcggc         <  2:335701/149‑1 (MQ=255)
                                                                                                                                        gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAg                                   >  7:325621/1‑118 (MQ=38)
                                                                                                                                        gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAg                                   <  8:325621/118‑1 (MQ=38)
                                                                                                                                        gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTAc    >  4:283715/1‑149 (MQ=255)
                                                                                                                                          ccAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCg  >  1:123060/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
AGTACCTGGTCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCG  >  NC_002947/2759615‑2759899

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: