Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,540 Δ1 bp 6.5% intergenic (+56/+38) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,5400C. 6.5% 398.7 / 24.5 138intergenic (+56/+38)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (63/60);  new base . (5/4);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 146 total aligned reads displayed.

GGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCC  >  NC_002947/5212392‑5212687
                                                                                                                                                     |                                                                                                                                                   
   aCCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                                                 <  2:3003742/150‑1 (MQ=255)
    ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAgg                                                                                                                                                >  2:232279/1‑150 (MQ=255)
     cAACGACGAAGCCGTCAACCGCCGGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGa                                                                                                                                               >  2:152816/1‑150 (MQ=255)
         gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                                           <  2:2884673/150‑1 (MQ=255)
            gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCggg                                                                                                                                          >  1:1940509/1‑148 (MQ=255)
            gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCggg                                                                                                                                          <  2:1940509/148‑1 (MQ=255)
            gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGtt                                                                                                                                        <  2:2918338/150‑1 (MQ=255)
               gCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAc                                                                                                                                     >  1:2774796/1‑150 (MQ=255)
                 cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                                                   <  1:2885309/150‑1 (MQ=255)
                 cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                                                   <  1:2099693/150‑1 (MQ=255)
                  gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGCAGCGGGTTTAcccc                                                                                                                                  >  1:2432042/1‑149 (MQ=255)
                       ccGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGaaa                                                                                                                             <  1:1032777/150‑1 (MQ=255)
                        cGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                                            >  2:2923296/1‑150 (MQ=255)
                           cTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                                         >  2:310620/1‑150 (MQ=255)
                                 aaGGAGAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAGTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                                                   <  1:2930143/150‑1 (MQ=255)
                                 aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGACCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                                                   >  2:1292750/1‑150 (MQ=255)
                                  aGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGc                                                                                                                  <  2:668681/150‑1 (MQ=255)
                                   ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                                                 >  2:1377213/1‑150 (MQ=255)
                                   ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                                                 >  1:1385117/1‑150 (MQ=255)
                                    gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGccc                                                                                                                <  1:2659499/150‑1 (MQ=255)
                                     aaaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                               <  2:2790300/150‑1 (MQ=255)
                                     aaaGGTTGCTGAGCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                               <  2:1085467/150‑1 (MQ=255)
                                         gTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcaca                                                                                                           <  1:232279/150‑1 (MQ=255)
                                            gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAg                                                                                                        <  1:2157476/150‑1 (MQ=255)
                                             cTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtc                                                                                                                                                     <  1:770805/104‑1 (MQ=255)
                                             cTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtc                                                                                                                                                     >  2:770805/1‑104 (MQ=255)
                                              tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCg                                                                                                      >  1:2026596/1‑150 (MQ=255)
                                               gATCAAGGACAAGATCGAACGGATCATTCATCCGCCTAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGa                                                                                                     <  1:2892102/150‑1 (MQ=255)
                                               gATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGa                                                                                                     >  1:717357/1‑150 (MQ=255)
                                                    aGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                                >  1:1357943/1‑150 (MQ=255)
                                                    aGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                                >  1:1284862/1‑150 (MQ=255)
                                                     ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTACACAGCGACCAgcgc                                                                                              >  2:2769859/1‑150 (MQ=255)
                                                      gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcgc                                                                                              >  1:2549265/1‑150 (MQ=255)
                                                         aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                                           <  2:688127/102‑1 (MQ=255)
                                                          aGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGTCATCTATTGAATCCACCGGCCTCTTCGCGGGCACGCCCAATCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCCCACACCGCGACCAGCGCTTTc                                                                                          >  1:1023340/1‑150 (MQ=39)
                                                           gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCa                                                                                         >  1:85956/1‑150 (MQ=255)
                                                            aTCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGACGGTGTAGGAGCgggc                                                                                                                                         >  2:1626596/1‑100 (MQ=37)
                                                               gAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                     >  1:957902/1‑150 (MQ=255)
                                                                  cGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcccc                                                                                  <  1:152816/150‑1 (MQ=255)
                                                                   ggATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCa                                                                                 >  1:1082596/1‑150 (MQ=255)
                                                                    gATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAggg                                                                                                                          >  2:1708503/1‑108 (MQ=255)
                                                                    gATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCAt                                                                                >  2:117863/1‑150 (MQ=255)
                                                                     aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                               >  1:2244748/1‑150 (MQ=255)
                                                                     aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCACGCGAAAGCCCGCTCATAGATCTTCGGGGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                               >  1:372909/1‑150 (MQ=255)
                                                                      tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcg                                                                                                                                                    >  1:3020468/1‑80 (MQ=255)
                                                                      tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcg                                                                                                                                                    <  2:3020468/80‑1 (MQ=255)
                                                                        attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCtt                                                                                            >  2:15458/1‑134 (MQ=255)
                                                                          tcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGt                                                                          <  2:1465679/150‑1 (MQ=255)
                                                                            atCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTcc                                                                        <  1:1684529/150‑1 (MQ=255)
                                                                                gCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                                >  1:437219/1‑122 (MQ=255)
                                                                                gCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGa                                                                    >  1:2682536/1‑150 (MQ=255)
                                                                                 cccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                                   <  2:1284862/150‑1 (MQ=255)
                                                                                  ccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAcc                                                                  >  1:409912/1‑150 (MQ=255)
                                                                                   cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGcttgc                                                                                                                                     <  1:2185836/82‑3 (MQ=17)
                                                                                    aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                                         >  2:47731/1‑93 (MQ=255)
                                                                                    aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                                         <  1:47731/93‑1 (MQ=255)
                                                                                      aaccTGCCAGCTATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTATACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACCCAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGt                                                              <  2:487316/148‑1 (MQ=38)
                                                                                      aTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGt                                                              <  1:1763713/150‑1 (MQ=255)
                                                                                       tCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTc                                                             >  2:2081198/1‑150 (MQ=255)
                                                                                       gccTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTc                                                             >  2:2103862/2‑150 (MQ=255)
                                                                                        ccTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCa                                                            <  1:117863/150‑1 (MQ=255)
                                                                                           gACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                         <  1:2081198/150‑1 (MQ=255)
                                                                                           gACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                         <  2:1357943/150‑1 (MQ=255)
                                                                                             cAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTaa                                                      <  1:2769859/150‑1 (MQ=255)
                                                                                                gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                                >  2:2213499/1‑106 (MQ=255)
                                                                                                gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                                <  1:2213499/106‑1 (MQ=255)
                                                                                                 cATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATcc                                                   <  2:957902/150‑1 (MQ=255)
                                                                                                       ttGAATCCACCGGCCTCTTCGCGGGCAAGCCAGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATc                                             >  1:992174/1‑150 (MQ=255)
                                                                                                               aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGTGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                                   <  1:634254/120‑1 (MQ=255)
                                                                                                                 cGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                                                   >  2:2748885/1‑70 (MQ=255)
                                                                                                                  ggCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACaa                                  >  1:211164/1‑150 (MQ=255)
                                                                                                                     ctctTCGCGGGCAAGCCCGCTCCTACACCGTCTTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCg                               >  2:553845/1‑150 (MQ=255)
                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTCCACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                               <  2:2256162/101‑1 (MQ=255)
                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                                       >  2:316837/1‑77 (MQ=255)
                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                               >  1:2256162/1‑101 (MQ=255)
                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCTTCGATCCCTACCGACAACCgc                              >  1:2619987/1‑150 (MQ=255)
                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgc                              <  1:1705413/150‑1 (MQ=255)
                                                                                                                        ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcag                            <  1:310620/150‑1 (MQ=255)
                                                                                                                           gcgGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGcagcag                         <  2:2426410/150‑1 (MQ=255)
                                                                                                                                cAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTCCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtct                    <  1:1527790/150‑1 (MQ=255)
                                                                                                                                 aaacccGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtctc                   >  1:2114323/4‑150 (MQ=38)
                                                                                                                                 aaGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTtctc                   >  1:654954/1‑150 (MQ=255)
                                                                                                                                  aacccGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGtt                      <  1:796690/144‑1 (MQ=255)
                                                                                                                                    cccGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCt                >  1:1772011/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                                      <  2:550134/95‑1 (MQ=38)
                                                                                                                                     ccGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                                      >  1:550134/1‑95 (MQ=38)
                                                                                                                                     ccGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCa                                                                 >  2:2604640/1‑100 (MQ=38)
                                                                                                                                     ccGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                                          >  1:1548069/1‑107 (MQ=38)
                                                                                                                                        cTCCTACACCGGCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATc            <  2:372909/150‑1 (MQ=255)
                                                                                                                                           cTACACCGTCGTAGGAGCGGGTTTACTCGCGAAAGGGCCTGCGGCCCTGACACAGCGACCAGCGCATTCAACGCCCCATCCAGAACCTCGACCAGGTCACGGTAATCCACGATCACTACCGACAACAGCAGCAGGTTCTCACTGATtccc         >  2:2178082/1‑150 (MQ=255)
                                                                                                                                           cTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATcccc         >  1:71923/1‑150 (MQ=255)
                                                                                                                                             acacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTca                                                                                                             <  2:1613819/48‑1 (MQ=255)
                                                                                                                                             acacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTca                                                                                                             >  1:1613819/1‑48 (MQ=255)
                                                                                                                                             acacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAt       >  1:1671762/1‑150 (MQ=255)
                                                                                                                                              cacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATcc                                                   <  1:272797/105‑1 (MQ=255)
                                                                                                                                                acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGTGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCa     >  1:51192/1‑150 (MQ=255)
                                                                                                                                                acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc                 >  2:260153/1‑138 (MQ=255)
                                                                                                                                                acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc                 <  1:260153/138‑1 (MQ=255)
                                                                                                                                                  cgtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAgg                                                                           <  2:776033/78‑1 (MQ=255)
                                                                                                                                                   aacgtACGAGCGGGTTTACCCGCGAAATGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCAccc  >  2:820678/3‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCC  >  NC_002947/5212392‑5212687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: