Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,772 T→C 7.0% T391A (ACC→GCC)  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,7720TC 7.0% 161.9 / 2.9 57T391A (ACC→GCC) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCA  >  NC_002947/2759624‑2759920
                                                                                                                                                      |                                                                                                                                                     
tCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGt                                                                                                                                                       <  1:247655/149‑1 (MQ=255)
    ccACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTgg                                                                                                                                                   <  6:100615/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‑‑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)
                                         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                                                                                                                                        <  6:275212/111‑1 (MQ=255)
                                                     cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGAtt                                                                                                                                        >  5:275212/1‑111 (MQ=255)
                                                     cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                                               >  4:348731/1‑120 (MQ=255)
                                                     cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAgccggcc                                                                                                                               <  3:348731/120‑1 (MQ=255)
                                                     cAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTg                                                                                                  >  5:15726/1‑148 (MQ=255)
                                                           aaTTGCTGCAAACTCCCGCTTTTCACCCTCCAGTAACGACATTACCATATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAg                                                                                             <  1:136555/149‑1 (MQ=255)
                                                           aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                                                    >  5:192481/1‑125 (MQ=255)
                                                           aaTGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCtt                                                                                                                    <  6:192481/125‑1 (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                                                                              >  2:372269/1‑149 (MQ=255)
                                                                        tCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAg                                                                              >  1:271252/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                                                                                                                      <  3:64793/86‑1 (MQ=255)
                                                                                                gACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGctc                                                                                                                      >  4:64793/1‑86 (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                                                        <  8:325621/118‑1 (MQ=38)
                                                                                                                               gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAg                                                        >  7:325621/1‑118 (MQ=38)
                                                                                                                               gCCCAGACCGATCCGCTGACCGGCCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTAc                         >  4:283715/1‑149 (MQ=255)
                                                                                                                                 ccAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCg                       >  1:123060/1‑149 (MQ=255)
                                                                                                                                                  tCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCg      <  2:123060/149‑1 (MQ=255)
                                                                                                                                                    ggTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCgg     <  3:282246/148‑1 (MQ=255)
                                                                                                                                                      tCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCa  >  5:253259/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCGCCCACGTCGTGTCCATGGCTATCATTTATCCGCTTGAAGTGGTCGATATCCAGCGCAATGGCTGCAAACTCCCGCTTTTCAGCCTCCAGTAACGACA‑‑GGCTGATTTCCAGGCTGCGGCGGTTGCCCAGGCCGGTCAGCGGATCGGTCTGGGCGTCCCGATTGAGCCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCA  >  NC_002947/2759624‑2759920

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: