Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,333,820 T→A 5.3% L251Q (CTG→CAG)  fadA → 3‑ketoacyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,333,8200TA 5.3% 236.1 / 2.1 75L251Q (CTG→CAG) fadA3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (2/2);  total (39/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 = 9.08e-01

ACAAGCTGACAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCGGTGCCCAAG  >  NC_002947/2333676‑2333968
                                                                                                                                                |                                                                                                                                                    
aCAAGCTGACAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTcc                                                                                                                                                  <  5:51721/149‑1 (MQ=255)
aCAAGCTGACAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGCCCGCCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTAGAAGGCGGTTGAGTCACGGCCCGCAATGCCAGTCAGCTGTcc                                                                                                                                                  >  8:101405/1‑149 (MQ=255)
    gCTGACAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATg                                                                                                                                              >  1:193604/1‑149 (MQ=255)
    gCTGACAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGACCGGCAATGCCAGTCAGCTGTCCGATg                                                                                                                                              <  1:251080/149‑1 (MQ=255)
        aCAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGgc                                                                                                                                            <  6:323721/147‑1 (MQ=255)
              gtcgtcAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGc                                                                                                                                    <  3:275709/149‑1 (MQ=255)
                  tcAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGc                                                                                                                                <  6:202341/149‑1 (MQ=255)
                    aGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGt                                                                                                                              >  1:288303/1‑149 (MQ=255)
                    aGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGt                                                                                                                              <  7:163155/149‑1 (MQ=255)
                     gCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTg                                                                                                                             <  5:247284/149‑1 (MQ=255)
                     gCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTg                                                                                                                             >  7:62280/1‑149 (MQ=255)
                      ccATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCAGTCCGATGGCGCCAGTGCCTGCGTGc                                                                                                                            >  8:289840/1‑149 (MQ=255)
                      ccATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCAGTCCGATGGCGCCAGTGCCTGCGTGc                                                                                                                            >  8:285391/1‑149 (MQ=255)
                         tGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGa                                                                                                                         <  1:85645/149‑1 (MQ=255)
                                gTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAgg                                                                                                                  >  4:43330/1‑149 (MQ=255)
                                  gcgcCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGt                                                                                                                              <  1:328834/135‑1 (MQ=255)
                                  gcgcCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGt                                                                                                                              >  2:328834/1‑135 (MQ=255)
                                   cgcCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTgc                                                                                                               >  6:155286/1‑149 (MQ=255)
                                   cgcCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTg                                                                                                                <  3:103960/148‑1 (MQ=255)
                                     ccTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAg                                                                                                                                       <  4:235468/123‑1 (MQ=255)
                                     ccTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAg                                                                                                                                       >  3:235468/1‑123 (MQ=255)
                                     ccTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTgcgc                                                                                                             >  1:169791/1‑149 (MQ=255)
                                      cTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCt                                                                                                            >  4:179190/1‑149 (MQ=255)
                                          cccAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCGCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTggcg                                                                                                        <  7:101405/149‑1 (MQ=255)
                                           ccAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTggcgg                                                                                                       >  4:7421/1‑149 (MQ=255)
                                            cAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTggcggc                                                                                                      >  4:46664/1‑149 (MQ=255)
                                             aGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGgcggcg                                                                                                     >  6:26649/1‑149 (MQ=255)
                                                 cGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCt                                                                                                 >  7:167479/1‑149 (MQ=255)
                                                  gAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTc                                                                                                >  3:109408/1‑149 (MQ=255)
                                                  gAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTc                                                                                                >  3:21249/1‑149 (MQ=255)
                                                     gggAACCGGCCGGCGACGACCCTGGCCGATCTGCAGCGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCggg                                                                                             >  2:188448/1‑149 (MQ=255)
                                                          ccggccggCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATcg                                                                                        <  7:63604/149‑1 (MQ=255)
                                                           cggccggcGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGAt                                                                                                                        >  6:250617/1‑116 (MQ=255)
                                                           cggccggcGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGAt                                                                                                                        <  5:250617/116‑1 (MQ=255)
                                                               cggcGACGACCCTGGCCGATCTGCAGGGGCTAAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGcc                                                                                   >  5:40899/1‑149 (MQ=255)
                                                                ggcGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCAGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGccc                                                                                  <  7:285391/149‑1 (MQ=255)
                                                                ggcGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCAGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGccc                                                                                  <  7:289840/149‑1 (MQ=255)
                                                                      gacCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGc                                                                            <  3:7421/149‑1 (MQ=255)
                                                                         ccTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCtt                                                                          >  1:210095/1‑148 (MQ=255)
                                                                                   cTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTAt                                                               >  4:70488/1‑149 (MQ=255)
                                                                                   cTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGG‑AGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATc                                                              >  1:88039/1‑149 (MQ=255)
                                                                                    tGCAGGGGCTCAAGCCAGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTAta                                                              >  5:275783/1‑148 (MQ=255)
                                                                                      cAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTa                                                                <  8:62280/145‑1 (MQ=255)
                                                                                           gCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGtgg                                                       >  1:336000/1‑149 (MQ=255)
                                                                                             tCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGtggct                                                     >  5:117420/1‑149 (MQ=255)
                                                                                               aaGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGtggctgg                                                   >  7:12683/1‑149 (MQ=255)
                                                                                                aGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGtggctggc                                                  <  2:169791/149‑1 (MQ=255)
                                                                                                 gCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGctggcc                                                 >  1:286114/1‑149 (MQ=255)
                                                                                                  cccGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGctggcct                                                >  5:46462/1‑149 (MQ=255)
                                                                                                    cGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGctggcctgg                                              >  4:172548/1‑149 (MQ=255)
                                                                                                            aaGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAg                                      <  6:117420/149‑1 (MQ=255)
                                                                                                            aaGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCAGGCCCCCGAg                                      <  2:288303/149‑1 (MQ=255)
                                                                                                                cGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGaaa                                  <  2:193604/149‑1 (MQ=255)
                                                                                                                   ttGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATgg                               <  2:286114/149‑1 (MQ=255)
                                                                                                                   ttGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATgg                               <  3:43330/149‑1 (MQ=255)
                                                                                                                     gtgtCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGG‑AGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGc                             <  2:88039/148‑1 (MQ=255)
                                                                                                                      tgtCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGc                             <  3:70488/148‑1 (MQ=255)
                                                                                                                      tgtCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCa                            >  7:4997/1‑149 (MQ=255)
                                                                                                                       gtCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCAt                           >  3:9493/1‑149 (MQ=255)
                                                                                                                        tCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCAt                           <  4:109408/148‑1 (MQ=255)
                                                                                                                        tCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATc                          >  2:152570/1‑149 (MQ=255)
                                                                                                                            ggccggcAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGtggctgg                                                   >  7:31424/1‑120 (MQ=255)
                                                                                                                            ggccggcAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGtggctgg                                                   <  8:31424/120‑1 (MQ=255)
                                                                                                                             gccggcAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGccc                     <  3:172548/149‑1 (MQ=255)
                                                                                                                              ccggcAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCg                    <  6:46462/149‑1 (MQ=255)
                                                                                                                                ggcAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGgt                  <  6:40899/149‑1 (MQ=255)
                                                                                                                                   aaTGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGtt               >  6:330439/1‑149 (MQ=255)
                                                                                                                                    aTGCCAGTCAGGTGTCCGATGGCGCCTGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTg              <  1:146860/149‑1 (MQ=255)
                                                                                                                                    aTGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTg              <  4:21249/149‑1 (MQ=255)
                                                                                                                                     tGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTgg             >  1:171566/1‑149 (MQ=255)
                                                                                                                                        cAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCgg          >  7:123760/1‑149 (MQ=255)
                                                                                                                                        cAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCgg          <  3:179190/149‑1 (MQ=255)
                                                                                                                                         aGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCGGt         >  2:161027/1‑149 (MQ=255)
                                                                                                                                          gTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCGGTg        <  5:155286/149‑1 (MQ=255)
                                                                                                                                           tCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCTCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGCCGGTGc       <  2:336000/149‑1 (MQ=255)
                                                                                                                                                tGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCGGTGCCCAAg  >  3:257991/1‑149 (MQ=255)
                                                                                                                                                tGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCGGTGCCCAAg  >  6:129276/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACAAGCTGACAGGTGTCGTCAGCCATGAAGAGGTGCGCCTTACCCAGGACGAAGGGAACCGGCCGGCGACGACCCTGGCCGATCTGCAGGGGCTCAAGCCCGTGGTCGAAGGCGGTTGTGTCACGGCCGGCAATGCCAGTCAGCTGTCCGATGGCGCCAGTGCCTGCGTGCTGATGGAAGGTGCGCTGGCGGCGCGCTCGGGCATCGCGGCCCTGGGGCTTTACCGCGGTATCGCCGTGGCTGGCCTGGCCCCCGAGGAAATGGGCATCGGCCCGGTGTTGGCGGTGCCCAAG  >  NC_002947/2333676‑2333968

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