Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 908,893 C→T 5.2% E106K (GAG→AAG)  cra ← catabolite DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947908,8930CT 5.2% 236.0 / 2.2 77E106K (GAG→AAG) cracatabolite DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base T (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

TGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGCCGGCTGC  >  NC_002947/908745‑909039
                                                                                                                                                           |                                                                                                                                                  
tgtgCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGcgctc                                                                                                                                                    >  1:280567/1‑149 (MQ=255)
    cGGGGTCCAGACGGCGGTCGATGGCAATGACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTg                                                                                                                                                >  1:212933/1‑149 (MQ=255)
       ggTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTct                                                                                                                                             <  7:79509/149‑1 (MQ=255)
          ccAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCTACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggc                                                                                                                                          >  4:120923/1‑149 (MQ=255)
          ccAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggc                                                                                                                                          >  2:355364/1‑149 (MQ=255)
          ccAGACGGCGGTCGATGGCAAGCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggc                                                                                                                                          <  7:218165/149‑1 (MQ=255)
           cAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggct                                                                                                                                         <  7:318315/149‑1 (MQ=255)
            aGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggctg                                                                                                                                        <  5:191027/149‑1 (MQ=255)
            aGACGGCGGTCGAAGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggctg                                                                                                                                        <  5:307785/149‑1 (MQ=255)
             gACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggctgg                                                                                                                                       <  2:212933/149‑1 (MQ=255)
             gACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctgg                                                                                                                                           <  4:80433/145‑1 (MQ=255)
             gACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctgg                                                                                                                                           >  3:80433/1‑145 (MQ=255)
                  cggTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATc                                                                                                                                  >  4:330925/1‑149 (MQ=255)
                   ggTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCg                                                                                                                                 >  4:121540/1‑149 (MQ=255)
                    gTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCgc                                                                                                                                >  7:2782/1‑149 (MQ=255)
                     tCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCgct                                                                                                                               <  6:106920/149‑1 (MQ=255)
                       gATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCgctgc                                                                                                                             >  7:81895/1‑149 (MQ=255)
                       gATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCgctgc                                                                                                                             <  6:87293/149‑1 (MQ=255)
                           gCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCgctgct                                                                                                                            >  8:172910/1‑146 (MQ=255)
                           gCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCgctgct                                                                                                                            <  7:172910/146‑1 (MQ=255)
                               tCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATc                                                                                                                     <  3:330925/149‑1 (MQ=255)
                                   cGGCTGGGCTTTGTACTTCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTTGCAACGATCAGCGCAGCGCATCGGCGCGCACTGTACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATc                                                                                                                     <  5:177530/145‑1 (MQ=255)
                                         gCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTTGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTgg                                                                                                            >  7:247463/1‑148 (MQ=255)
                                         gCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTTGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTgg                                                                                                            >  7:252369/1‑148 (MQ=255)
                                         gCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGt                                                                                                           <  2:252327/149‑1 (MQ=255)
                                            tttGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggc                                                                                                                                          >  3:38347/1‑115 (MQ=255)
                                            tttGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTctggc                                                                                                                                          <  4:38347/115‑1 (MQ=255)
                                            tttGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGc                                                                                                        >  3:203040/1‑149 (MQ=255)
                                                 ccTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCgg                                                                                                   >  4:326911/1‑149 (MQ=255)
                                                  cTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCggg                                                                                                  >  1:182600/1‑149 (MQ=255)
                                                  cTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCggg                                                                                                  >  8:263719/1‑149 (MQ=255)
                                                   tGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGgc                                                                                                 >  6:301464/1‑149 (MQ=255)
                                                   tGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGgc                                                                                                 >  7:247866/1‑149 (MQ=255)
                                                   tGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGgc                                                                                                 >  1:55773/1‑149 (MQ=255)
                                                           gcgGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTgcgc                                                                                         <  3:120923/149‑1 (MQ=255)
                                                           gcgGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTgcg                                                                                          <  2:182600/148‑1 (MQ=255)
                                                           gcgGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTgc                                                                                           <  2:194017/147‑1 (MQ=255)
                                                             ggTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCt                                                                                       <  8:81895/149‑1 (MQ=255)
                                                               tAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTg                                                                                     <  1:277501/149‑1 (MQ=255)
                                                                     tCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCa                                                                                <  1:355364/148‑1 (MQ=255)
                                                                     tCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCa                                                                                <  8:238703/148‑1 (MQ=255)
                                                                       aTCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGtt                                                                             >  7:105990/1‑149 (MQ=255)
                                                                          ttctGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCt                                                                          >  1:92423/1‑149 (MQ=255)
                                                                          ttctGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTACt                                                                          >  8:28238/1‑149 (MQ=255)
                                                                           tctGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAACAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCtt                                                                         >  6:40707/1‑149 (MQ=255)
                                                                               gcggcCGGCAACTGGCAACGAACAGCGGCTCGCAGCGGCGCGCACGTATCAGCCCTGACCTCTTCCACACGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAg                                                                               <  1:123078/146‑1 (MQ=255)
                                                                               gcggcAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTg                                                                        <  4:161200/146‑1 (MQ=255)
                                                                                 ggcaggcaACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAAt                                                                   <  8:2782/149‑1 (MQ=255)
                                                                                    aggcaACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCg                                                                >  2:93262/1‑149 (MQ=255)
                                                                                     ggcaACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTTGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCgg                                                               <  8:252369/149‑1 (MQ=255)
                                                                                     ggcaACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTTGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCgg                                                               <  8:247463/149‑1 (MQ=255)
                                                                                     ggcaACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCgg                                                               <  8:247866/149‑1 (MQ=255)
                                                                                              cAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTg                                                      >  7:136553/1‑149 (MQ=255)
                                                                                               aaCGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTgg                                                     >  4:312303/1‑149 (MQ=255)
                                                                                               aaCGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTgg                                                     >  2:244651/1‑149 (MQ=255)
                                                                                                aCGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTggg                                                    <  6:174244/149‑1 (MQ=255)
                                                                                                aCGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTggg                                                    <  7:158141/149‑1 (MQ=255)
                                                                                                 cGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCGGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTgggg                                                   <  7:158/149‑1 (MQ=255)
                                                                                                       gcgcATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTtctc                                              >  2:213202/1‑148 (MQ=255)
                                                                                                       gcgcATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTcc                                             <  3:326911/149‑1 (MQ=255)
                                                                                                                     cgcgcACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTTCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAAt                               <  1:93262/149‑1 (MQ=255)
                                                                                                                      gcgcACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAAt                               >  4:46017/1‑148 (MQ=255)
                                                                                                                        gcACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGaa                            <  3:121540/149‑1 (MQ=255)
                                                                                                                        gcACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTACTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGaa                            <  7:28238/149‑1 (MQ=255)
                                                                                                                        gcACGGAACAG‑‑TTG‑‑‑‑C‑GGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGaa                            <  4:203040/149‑1 (MQ=255)
                                                                                                                          aCGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGTGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGc                          >  8:345444/1‑149 (MQ=255)
                                                                                                                            ggAACCG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCGTGCGCCCGGTGCTTGGCAATGCGGGCGTAGCTGGGGGTCGCCAGCTCCGCGAGAATGAAcccc                        >  8:342292/1‑149 (MQ=255)
                                                                                                                              aaCAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGc                                                             <  6:66838/110‑1 (MQ=255)
                                                                                                                              aaCAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGc                                                             >  5:66838/1‑110 (MQ=255)
                                                                                                                               aCAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAgg                     <  2:280567/149‑1 (MQ=255)
                                                                                                                               aCAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAgg                     <  2:92423/149‑1 (MQ=255)
                                                                                                                               aCAG‑‑GTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAgg                     <  5:301464/149‑1 (MQ=255)
                                                                                                                                cAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAggg                    <  3:46017/149‑1 (MQ=255)
                                                                                                                                  g‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTg                  >  2:236185/1‑149 (MQ=255)
                                                                                                                                               gcAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCggg                                    >  7:286150/1‑125 (MQ=255)
                                                                                                                                                   ttGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGcc        >  2:41640/1‑149 (MQ=255)
                                                                                                                                                   ttGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGcc        >  1:171517/1‑149 (MQ=255)
                                                                                                                                                      gcgcTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGCCGGc     >  2:321813/1‑149 (MQ=255)
                                                                                                                                                        gctcgctGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCATGGTGCGGGTATGCCGGCTg   >  3:232697/1‑149 (MQ=255)
                                                                                                                                                        gctcgctGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGACGGCTg   >  3:237442/1‑149 (MQ=255)
                                                                                                                                                         ctcgctGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGCCGGCTgc  <  2:55773/149‑1 (MQ=255)
                                                                                                                                                           |                                                                                                                                                  
TGTGCGGGGTCCAGACGGCGGTCGATGGCAATCACCGGCAGGCCTTTGTCCTGCAGCTCGCGGTAGCTGTCATCTTCTGGCGGCAGGCAACTGGCAACGAACAGCGCATCGCATCGGCGCGCACGGAACAG‑‑TTG‑‑‑‑C‑TGCAGTTGGCGCTCGCTGTCTGGCTGGTCATCGCTGCTGGCGATCAGCAGCTGGTAGCCGCGGGCGCGTGCGCCTTGCTCCAGTTGCTTGGCAATGCGGGCGTAGCTGGGGTTCTCCAGATCCGGGAGAATGAAGCCCAGGGTGCGGGTATGCCGGCTGC  >  NC_002947/908745‑909039

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: ‑