Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,268,065 A→T 5.5% Y179F (TAT→TTT) ‡ azoR1 → FMN‑dependent NADH‑azoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,268,0650AT 5.5% 227.8 / 2.7 73Y179F (TAT→TTT) ‡azoR1FMN‑dependent NADH‑azoreductase
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base T (2/2);  total (37/36)
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.64e-01

GGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTT  >  NC_002947/3267917‑3268210
                                                                                                                                                    |                                                                                                                                                 
ggCAACACGGAACTGGTGCTGGTTCCCACTGCGGGGGGCCTGCACGCTCGCCAGCCGACTGCCGCCGGGCATGAACACTTCCTGAAACTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTa                                                                                                                                                   <  8:294168/149‑1 (MQ=255)
      aaGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGccc                                                                                                                                             >  8:239059/1‑149 (MQ=255)
      aaGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGccc                                                                                                                                             >  6:339503/1‑149 (MQ=255)
      aaGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGccc                                                                                                                                             >  3:20832/1‑149 (MQ=255)
      aaGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGccc                                                                                                                                             >  6:145282/1‑149 (MQ=255)
       cgAAAGGGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGTCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCg                                                                                                                                            <  6:76365/148‑1 (MQ=255)
         aaaGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCATAGGCCAGGCCGTGGGCCCGGa                                                                                                                                          <  2:280781/149‑1 (MQ=255)
          aaGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGa                                                                                                                                          <  6:129673/148‑1 (MQ=255)
            gtggtgCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGagca                                                                                                                                       <  4:20832/149‑1 (MQ=255)
               gtgCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAgcagcg                                                                                                                                    >  3:136457/1‑149 (MQ=255)
                  cTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAgcagcgca                                                                                                                                  <  7:239059/148‑1 (MQ=255)
                       ttCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAgg                                                                                                                            >  6:219353/1‑149 (MQ=255)
                        tCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAgg                                                                                                                            <  4:231203/148‑1 (MQ=255)
                           aCTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcgat                                                                                                                        >  7:64174/1‑149 (MQ=255)
                           aCTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGcgat                                                                                                                        >  6:308010/1‑149 (MQ=255)
                                gggTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGANTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGAtg                                                                                                                   >  1:128414/1‑149 (MQ=255)
                                   tGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGAtgctg                                                                                                                <  8:64174/149‑1 (MQ=255)
                                      ccTGCACGCTGGCCAGCCGACTAGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGca                                                                                                              >  5:335447/1‑148 (MQ=255)
                                       cTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGAtgctgta                                                                                                              >  7:178436/1‑145 (MQ=255)
                                       cTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGca                                                                                                              >  7:387082/1‑147 (MQ=255)
                                       cTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCATAGGCCAGGCCGTGGGCCCGGa                                                                                                                                          >  8:301477/1‑119 (MQ=38)
                                          cACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgc                                                                                                         >  8:162782/1‑149 (MQ=255)
                                          cACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgc                                                                                                         <  5:145282/149‑1 (MQ=255)
                                             gCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGAGTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGgcg                                                                                                        >  7:206931/1‑147 (MQ=255)
                                                 gCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGTACAGGCGCAGAt                                                                                                   <  8:178436/148‑1 (MQ=255)
                                                 gCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGAtt                                                                                                  <  1:305519/149‑1 (MQ=255)
                                                    aGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGc                                                                                                >  7:97518/1‑148 (MQ=255)
                                                     gCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGc                                                                                                >  2:376991/1‑147 (MQ=255)
                                                         aCTGGCGCCGGGCATGAAGTCTTCCTGAAAGTGTTCCTGGGCTTCATCGGGTTTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGa                                                                                          <  3:248419/149‑1 (MQ=255)
                                                         aCTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGa                                                                                          <  3:9260/149‑1 (MQ=255)
                                                             gcgcCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCt                                                                                       <  4:136457/148‑1 (MQ=255)
                                                                     cATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGcaca                                                                                                            <  1:360143/119‑1 (MQ=255)
                                                                     cATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGcaca                                                                                                            >  2:360143/1‑119 (MQ=255)
                                                                        gAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGc                                                                           >  2:270515/1‑149 (MQ=255)
                                                                         aaGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGc                                                                           <  6:335447/148‑1 (MQ=255)
                                                                         aaGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGc                                                                           <  2:22767/148‑1 (MQ=255)
                                                                            aCTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGACCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGa                                                                       >  7:282620/1‑149 (MQ=255)
                                                                              ttCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGTCCCACGTCCTGGCCTATTGCAAGGAGCAGCGCAGCCAGGCGATCGATTCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGTCTGCATGAt                                                                      >  3:342669/1‑148 (MQ=255)
                                                                              ttCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGAt                                                                      >  8:253413/1‑148 (MQ=255)
                                                                                 cTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGcc                                                                                                                                              <  7:103731/73‑1 (MQ=255)
                                                                                 cTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGcc                                                                                                                                              >  8:103731/1‑73 (MQ=255)
                                                                                 cTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCaa                                                                   <  5:35905/148‑1 (MQ=255)
                                                                                 cTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAg                                                                  <  5:219353/149‑1 (MQ=255)
                                                                                  tGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAgg                                                                 >  8:359793/1‑149 (MQ=255)
                                                                                  tGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCATAGGCCAGGCCGTGGGCCCgg                                                                                                                                           >  5:33037/1‑75 (MQ=21)
                                                                                  tGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCATAGGCCAGGCCGTGGGCCCgg                                                                                                                                           <  6:33037/75‑1 (MQ=21)
                                                                                   gAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAgg                                                                 >  2:341854/1‑148 (MQ=255)
                                                                                          ttCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGtt                                                          <  7:162782/148‑1 (MQ=255)
                                                                                            ccTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTg                                                        <  5:308010/148‑1 (MQ=255)
                                                                                             cTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCg                                                                                <  8:26821/123‑1 (MQ=255)
                                                                                             cTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCg                                                                                >  7:26821/1‑123 (MQ=255)
                                                                                             cTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGt                                                      <  5:339503/149‑1 (MQ=255)
                                                                                                ggCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGc                                                   <  8:387082/149‑1 (MQ=255)
                                                                                                      aTCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGt                                             >  4:85117/1‑149 (MQ=255)
                                                                                                            aTTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCt                                       >  5:373573/1‑149 (MQ=255)
                                                                                                               aCTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGt                                     <  7:359793/148‑1 (MQ=255)
                                                                                                                cTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTc                                    <  8:97518/148‑1 (MQ=255)
                                                                                                                cTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTc                                    <  3:353252/148‑1 (MQ=255)
                                                                                                                cTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCg                                   <  6:373573/149‑1 (MQ=255)
                                                                                                                   aCCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGa                                                                                          <  6:322723/91‑1 (MQ=255)
                                                                                                                   aCCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGa                                                                                          >  5:322723/1‑91 (MQ=255)
                                                                                                                   aCCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGaaa                                >  3:382874/1‑149 (MQ=255)
                                                                                                                   aCCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGaaa                                <  5:184478/149‑1 (MQ=255)
                                                                                                                        gAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAacac                             <  3:85117/147‑1 (MQ=255)
                                                                                                                                  gggCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTcc                  >  6:200022/1‑148 (MQ=255)
                                                                                                                                    gCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAAc               <  1:341854/149‑1 (MQ=255)
                                                                                                                                      ccACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCt             >  2:347562/1‑149 (MQ=255)
                                                                                                                                        aCGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCACGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGt           <  2:128414/149‑1 (MQ=255)
                                                                                                                                           gcctggcctATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGa        <  1:270515/149‑1 (MQ=255)
                                                                                                                                            cctggcctATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGaa       >  3:65088/1‑149 (MQ=255)
                                                                                                                                            cctggcctATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGaa       <  1:204892/149‑1 (MQ=255)
                                                                                                                                             ctggcctATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAg      >  8:304825/1‑149 (MQ=255)
                                                                                                                                                   tATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGtt  >  5:127289/1‑147 (MQ=255)
                                                                                                                                                   tATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGtt  >  2:265349/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCAACAAGAAAGTGGTGCTGGTTTCCACTGCGGGTGGCCTGCACGCTGGCCAGCCGACTGGCGCCGGGCATGAAGACTTCCTGAAAGTGTTCCTGGGCTTCATCGGGATTACTGACCTGGAAATCGTCCGGGCCCACGGCCTGGCCTATGGCCCGGAGCAGCGCAGCCAGGCGATCGATGCTGCACAGGCGCAGATTGCCAGCGAGCTGTTTGCGGCTGCCTGATCAAGGGAGCGTTTGGTTGCAAGCGTGGCCCTGTCGAAACACAAGGCCACTCCAACCTGTAGAAGCGTT  >  NC_002947/3267917‑3268210

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