Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,746,379 C→T 6.7% G555G (GGC→GGT paoC → aromatic aldehyde dehydrogenase molybdopterin‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,746,3790CT 6.7% 219.2 / 3.7 74G555G (GGC→GGTpaoCaromatic aldehyde dehydrogenase molybdopterin‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base T (3/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

ACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGTCGACGAG  >  NC_002947/3746233‑3746526
                                                                                                                                                  |                                                                                                                                                   
aCGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCa                                                                                                                                                    <  6:96608/148‑1 (MQ=255)
 cGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCATGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCCCGCGTGCAgg                                                                                                                                                  <  5:185380/149‑2 (MQ=255)
  ggTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAACTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGc                                                                                                                                                  >  4:300519/1‑148 (MQ=255)
      tcgtcgGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTGCCCGCGTGCAGGTAGCCATCg                                                                                                                                             <  4:369085/149‑1 (MQ=255)
        gtcggtcgGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGc                                                                                                                                            <  1:72065/148‑1 (MQ=255)
         tcggtcggTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCt                                                                                                                                          >  5:71771/1‑149 (MQ=255)
          cggtcggtTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTg                                                                                                                                         <  6:112648/149‑1 (MQ=255)
           ggtcggtTCGGCCGTGCGCCATGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTg                                                                                                                                         <  1:254837/148‑1 (MQ=255)
           ggtcggtTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTgg                                                                                                                                        >  1:337694/1‑149 (MQ=255)
                gtTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTg                                                                                                                                   >  5:55661/1‑149 (MQ=255)
                gtTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTg                                                                                                                                   <  5:26209/149‑1 (MQ=255)
                     gCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGa                                                                                                                              >  4:224931/1‑149 (MQ=255)
                     gCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACATGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGa                                                                                                                              <  4:247603/149‑1 (MQ=255)
                       cGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGAca                                                                                                                            <  3:308909/149‑1 (MQ=255)
                         tGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGAcatc                                                                                                                          <  6:49010/149‑1 (MQ=255)
                            gcCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACatcatc                                                                                                                       <  5:310432/149‑1 (MQ=255)
                            gcCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACatcatc                                                                                                                       <  8:199211/149‑1 (MQ=255)
                              cAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGc                                                                                                                     <  1:12234/149‑1 (MQ=255)
                               aGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGc                                                                                                                     >  4:5265/1‑148 (MQ=255)
                                ggCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCa                                                                                                                   <  3:64299/149‑1 (MQ=255)
                                 gCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAg                                                                                                                  >  5:283081/1‑149 (MQ=255)
                                   cTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAgcg                                                                                                                <  7:226821/149‑1 (MQ=255)
                                    tGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAgcgc                                                                                                               >  1:125092/1‑149 (MQ=255)
                                      ccAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGccc                                                                                                             >  5:102923/1‑149 (MQ=255)
                                       cAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGTCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGCCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACAACAACGCCAGCGcccc                                                                                                            >  5:284322/1‑149 (MQ=255)
                                                cGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGt                                                                                                     >  7:196460/1‑147 (MQ=255)
                                                cGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGt                                                                                                     <  2:334013/147‑1 (MQ=255)
                                                    ccAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGtt                                                                                                <  6:306817/148‑1 (MQ=255)
                                                     cAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTc                                                                                               >  3:145199/1‑148 (MQ=255)
                                                     cAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTc                                                                                               >  1:263047/1‑148 (MQ=255)
                                                       aaGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGaa                                                                                            <  6:71771/149‑1 (MQ=255)
                                                             tGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCa                                                                                       >  3:261563/1‑148 (MQ=255)
                                                                  cATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGa                                                                                 <  5:372802/149‑1 (MQ=255)
                                                                        cgcAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGa                                                                           <  7:326516/149‑1 (MQ=255)
                                                                         gcAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAg                                                                          >  6:68633/1‑149 (MQ=255)
                                                                         gcAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAg                                                                          >  2:43249/1‑149 (MQ=255)
                                                                               tGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGgc                                                                    >  3:350639/1‑149 (MQ=255)
                                                                                 tACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGgcgc                                                                  >  6:250329/1‑149 (MQ=255)
                                                                                 tACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGgcg                                                                   <  1:43249/148‑1 (MQ=255)
                                                                                    ccccAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGgcgc                                                                  >  5:37619/1‑146 (MQ=255)
                                                                                     cccAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGa                                                              >  2:211094/1‑149 (MQ=255)
                                                                                      ccAAGTGGCCAGCACGGCGCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGAc                                                             <  6:284322/149‑1 (MQ=255)
                                                                                      ccAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGAc                                                             >  8:338047/1‑149 (MQ=255)
                                                                                             gCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCGCGGCGCCAGACGCCAAcc                                                      <  2:263047/149‑1 (MQ=255)
                                                                                             gCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCCCCCCCGAAGGTGAGTAGGAAGTGCAAATCGACGCCGAGCCCGCCGCCAGACGCCAAcc                                                      >  8:105212/1‑149 (MQ=255)
                                                                                                      gCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCAc                                             <  2:337694/149‑1 (MQ=255)
                                                                                                      gCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCAc                                             <  4:350639/149‑1 (MQ=255)
                                                                                                      gCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCAc                                             >  3:11802/1‑149 (MQ=255)
                                                                                                              gCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCGGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGaa                                                                                            >  2:196932/1‑94 (MQ=255)
                                                                                                                 gctCAGTTTCGCCGCTGGCTCCCTGCACGCGTGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGACGGGGAGTTCGCAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTGCTCCACCGCTACCCAtt                                  <  8:251090/147‑1 (MQ=255)
                                                                                                                 gATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAAc                                                       <  2:218282/128‑1 (MQ=255)
                                                                                                                 gATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAAc                                                       >  1:218282/1‑128 (MQ=255)
                                                                                                                  aTCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGACAGCGCCCCCGAAGGGGAGTTCGAAGTGCAAATCGACCCCGCGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTc                                 >  7:132983/1‑149 (MQ=255)
                                                                                                                    cATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCgg                               >  4:282053/1‑149 (MQ=255)
                                                                                                                     aTTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCgg                               <  3:300519/148‑1 (MQ=255)
                                                                                                                       tttCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAg                            >  7:204189/1‑149 (MQ=255)
                                                                                                                            ccGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTc                       <  3:282053/149‑1 (MQ=255)
                                                                                                                              gATGGCTGCCCGCGTGCAGGTAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGa                                                              >  2:6261/1‑108 (MQ=38)
                                                                                                                              gATGGCTGCCCGCGTGCAGGTAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGa                                                              <  1:6261/108‑1 (MQ=38)
                                                                                                                              gATGGCTGCCCGCGTGCAGGTAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTcgc                     >  1:357160/1‑149 (MQ=255)
                                                                                                                               aTGGCTGCCCGCGTGCAGGTAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGACCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTcgcg                    >  6:338103/1‑149 (MQ=255)
                                                                                                                               aTGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTcgcg                    <  2:125092/149‑1 (MQ=255)
                                                                                                                                tGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCgg                   <  5:68633/149‑1 (MQ=255)
                                                                                                                                  gCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGaa                 <  7:257547/149‑1 (MQ=255)
                                                                                                                                  gCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGaa                 <  4:353658/149‑1 (MQ=255)
                                                                                                                                  gCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGaa                 <  6:37619/149‑1 (MQ=255)
                                                                                                                                   cTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCg                                <  1:75728/133‑1 (MQ=255)
                                                                                                                                   cTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCg                                >  2:75728/1‑133 (MQ=255)
                                                                                                                                   cTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAg                <  7:338047/149‑1 (MQ=255)
                                                                                                                                        tGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCAc           >  7:218345/1‑149 (MQ=255)
                                                                                                                                          cACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGt         >  7:269922/1‑149 (MQ=255)
                                                                                                                                             gcgGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGTcga      <  8:204189/149‑1 (MQ=255)
                                                                                                                                                ggCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGTcgacga   <  6:283081/149‑1 (MQ=255)
                                                                                                                                                 gCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGTCGACGAg  <  5:250329/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGTCGACGAG  >  NC_002947/3746233‑3746526

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: