Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,746,372 A→G 6.5% H553R (CAC→CGC) ‡ 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,3720AG 6.5% 223.0 / 4.5 78H553R (CAC→CGC) ‡paoCaromatic aldehyde dehydrogenase molybdopterin‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (34/38);  new base G (3/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

GGATCGTTTACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGT  >  NC_002947/3746224‑3746519
                                                                                                                                                    |                                                                                                                                                   
ggATCGTTTACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGTCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCa                                                                                                                                                     >  1:334013/1‑149 (MQ=255)
     gTTTACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCAcgcg                                                                                                                                                 <  6:189367/148‑1 (MQ=255)
     gTTTACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCgg                                                                                                                                                <  5:115989/149‑1 (MQ=255)
     gTTTACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCgg                                                                                                                                                <  2:107105/149‑1 (MQ=255)
      tttACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCggg                                                                                                                                               >  3:353658/1‑149 (MQ=255)
        tACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCggg                                                                                                                                               >  2:72065/1‑147 (MQ=255)
         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                                                                                                                <  8:199211/149‑1 (MQ=255)
                                     gcCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACatcatc                                                                                                                <  5:310432/149‑1 (MQ=255)
                                      ccAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCAcgc                                                                                                                                                  >  2:374332/1‑114 (MQ=255)
                                      ccAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCAcgc                                                                                                                                                  <  1:374332/114‑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                                                                   >  2:43249/1‑149 (MQ=255)
                                                                                  gcAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAg                                                                   >  6:68633/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                                      <  4:350639/149‑1 (MQ=255)
                                                                                                               gCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCAc                                      <  2:337694/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                                                >  1:218282/1‑128 (MQ=255)
                                                                                                                          gATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAAc                                                <  2:218282/128‑1 (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                                                       <  1:6261/108‑1 (MQ=38)
                                                                                                                                       gATGGCTGCCCGCGTGCAGGTAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGa                                                       >  2:6261/1‑108 (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          <  6:37619/149‑1 (MQ=255)
                                                                                                                                           gCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGaa          <  4:353658/149‑1 (MQ=255)
                                                                                                                                           gCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGaa          <  7:257547/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)
                                                                                                                                                    |                                                                                                                                                   
GGATCGTTTACGGTGTCGTCGGTCGGTTCGGCCGTGCGCCAGGCCTGCCAGGTACTGCGTACCAAAGTGCTGGAACATGCGCGCAGGATGTACCCCCAAGTGGCCAGCACGGCCCGGGAGCAGATCATTTTCGCCGATGGCTACCTGCACGCGGGCAGCCATCGCCTGGCCTTGGCGGACATCATCGCCAGCGCCCCCGAAGGTGAGTTCGAAGTGCAAATCGACGCCGAGCCCGGCGCCAGACGCCAACCGTTCTCCACCGCTACCCATTCGGCAGTATTCGCGGAAGTCCACGT  >  NC_002947/3746224‑3746519

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